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This  is  a report  of  an  in-house  effort  conducted  under  Project 
1368,  "Advanced  Structures  for  Military  Aerospace  Vehicles"/  Task 
136802,  "Structural  Integration  for  Advanced  Military  Aerospace 
Vehicles"/Work  Unit  13680229,  "Bird-Proof  Windshields  for  Military 
Aircraft"  in  the  Structures  Development  and  Integration  Group, 
Structural  Development  Branch,  Structures  Division  AF  Flight 
Dynamics  Laboratory.  The  principal  investigator  for  this  work  was 
Mr.  Werner  R.  Jansen,  Aerospace  Engineer. 

The  author  gratefully  acknowledges  the  interests  and  advices  of 
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SECTION  I 
INTRODUCTION 

Efficient  bird-proof  design  is  a requirement  for  military  air- 
craft to  reduce  losses  due  to  bird  strike  and  to  reduce  acquisition  cost 
by  reducing  aircraft  buys  required  to  maintain  the  fleet  size. 

A stringent  Air  Force  requirement  is  to  operate  ai. craft  in  take-off, 
landing,  and  low  level,  high  speed  penetration  modes  where  maximum  crew 
protection  is  afforded  for  mission  accomplishment.  Windshields  and 
canopies  of  aircraft  are  continuously  vulnerable.  In  1971,  a F-lli 
airplane  was  lost  as  a result  of  a bird  strike  (Reference  1,2,  and  Table  1) . 

Aircraft  bird  strike  hazard  programs  may  fit  into  two  categories, 

(1)  concerned  with  the  overall  bird  strike  problem  (collision  avoidance), 
and  (2)  with  the  design  of  the  aircraft  itself  to  withstand  the  damaging 
effects  of  bird  strikes.  Bird  impact  resistance  of  aircraft  transparencies 
under  dynamic  conditions  is  the  theme  of  this  report. 

Development  of  conceptionally  new  bird  proof  windshield  systems  is 
costly.  Design  modifications  of  existing  windshield  systems  that  otherwise 
perform  perfectly  may  be  elaborate  too.  Therefore,  one  should  consider 
other  means  of  bird  strike  alleviation.  The  objective  of  Reference  3 was 
to  determine  the  general  feasibility  of  using  the  deflected  air  stream 
or  pressure  field  caused  by  vanes,  fuselage  shaping,  or  shock  waves  to 
prevent  bird  strikes  on  aircraft  windshields.  This  study  concluded  that 
bird  strike  prevention  by  means  of  deflected  air  flow  using  vanes,  fuselage 
shaping,  or  shock  waves,  was  not  a feasible  concept. 
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In  this  report,  techniques  for  reducing  bird  damage  to  USAF  aircraft 
and  air  crew  compartments  are  examined.  Design  concepts  that  would  result 
in  cockpit  enclosures  which  can  withstand  bird  impact  are  analyzed.  Consid- 
erations include  force-time  loading  factors,  material  response,  attachment 
design,  and  structural  transition  zones. 

1.  BIRD  STRIKE  ENVIRONMENT 

An  international  conference  on  bird  strike  on  aircraft  held  in  Sept. 
1969  revealed  that  bird  incidents  involving  aircraft  in  the  vicinity  of  air- 
ports had  significantly  decreased  (References  4 and  5) . ■ The  only  measures 
offering  any  real  assurance  of  reducing  the  enroute  hazard  are  those  of  bird- 
proofing the  aircraft.  This  approach  consists  of  increasing  the  strength 
and  the  energy-absorbing  capabilities  of  the  vulnerable  parts  of  the  aircraft 
that  are  likely  to  be  subjected  to  bird  impacts.  Reference  2 Indicates  that 
the  other  approach  to  the  problem,  that  of  trying  to  avoid  the  birds,  also 
appears  to  be  a possible  solution.  Although  a great  deal  of  effort  has  been 
devoted  to  this  concept,  success  is  not  in  sight. 

In  response  to  AF  Systems  Command  Document  "Coordinating  Committee 
on  Bird  Strike"  of  1972  and  Technical  Need  TN4-10  "Integrated  Design  Concepts 
to  Provide  Aircraft  Windscreen  Bird  Strike  Protection",  a design  concept  for 
a bird-proof,  high-visibility  windshield  by  means  of  controlled  shock  isola- 
tion lias  been  proposed  ("Proposal  of  Bird  Proof  Windshield",  AFFDL/FBS, 

March  1972).  The  AFSC  document  suggested  a two  phase  effort  attempting  to 
treat  the  windscreen  bird  strike  problem  as  an  integrated  design  approach 
using  and  optimizing  various  specific  techniques,  such  as  deflectors,  wind- 
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shield  shaping,  specific  operations,  materials,  structural  strength,  etc., 
in  order  to  achieve  the  total  objective  of  adequate  windscreen  bird  proof 
designs. 

This  concept  was  pursued  to  determine  feasibility.  A detailed  DDC 
literature  search  (Reference  6)  showed  no  analysis  methods  and  data  with 
respect  :o  shock-absorbing  windows. 

2.  GENERAL  WINDSHIELD  CRITERIA 

Design  criteria  for  a tandem  side  opening  canopy,  within  the  scope  of 
structural  finite  element  analysis  conducted  prior  to  1972,  lacked  the 
requirement  of  proofing  against  bird  strike.  The  magnitude  of  the  design 
operating  stresses  in  a canopy  structure  must  be  reliably  known  to  insure 
both  minimum  weight  and  structural  integrity.  Minimum  weight  is  required 
since  the  canopy  ejection  thruster  size  increases  as  the  canopy  weight 
increases.  And,  when  pilot  ejection  through  the  canopy  glass  is  required, 
the  glass  must  be  as  thin  as  possible  (Reference  7) . The  windshield 
must  also  satisfy  strength  requirements  in  new  load  carrying  concepts 
as  part  of  the  external  shell  of  the  fuselage,  subjected  to  the  pressure 
differential  due  to  aerodynamic  and  pressurization  loads  totaling  several 
tons  (Reference  8) . 

Visibility  thru  the  windshield  at  all  times  during  inclement  weather 
is  an  obvious  requirement  for  aircraft  (Table  II) . In  addition  to  the 
strength  and  optical  requirements,  electrical  and  thermal  properties  must  be 
considered.  A satisfactory  electrical  conducting  film  must  have  a high  level 
of  light  transmission  and  good  electrical  properties  for  heating  and  demisting 


*N 


4 


AFFDL-TR-74-155 


TABLE  II 

WINDSHIELD  DESIGN  CRITERIA,  CONSIDERATIONS  AND  REQUIREMENTS 


Life  Cycle  General  Design  Criteria 


Structural  Static  and  Fatigue  Strength 
Environmental  Fracture  Toughness 
Impact  Resistance  Against  Bird,  Hail,  Bullet 
Cyclic  Pressure  and  Temperature  Service  Capability 
Weathering  Resistance 


Optics  and  Field  of  Vision  Qualification 


Design  Considerations* 
Weight 

Abrasion  Resistance 
Bird  Impact  Resistance 
Cabin  Pressure  Retention 
Ice  Protection 
Fail-Safe  Characteristics 
Fog  Protection 
Size 
Contour 

(*  From  Reference  9) 


Design  Requirements* 
Temperature  Requirements 
Heat/Power  Requirements 
Thermal  Design 
Electrical  Design 
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(Reference  10).  Its  adherence  and  lamination  strengths  are  certainly  adjunc- 
tive requirements. 

Even  bullet  impact  has  to  be  considered  as  a plausible  risk,  either 
from  outside  the  aircraft  or  fired  from  within  (Reference  8) . Hail  stone 
impact,  with  masses  up  to  bird  size  masses,  also  needs  to  be  considered 
(Reference  11) , Reference  12  notes  that  there  is  no  test  evidence  to  com- 
pare the  damage  to  windscreens  caused  by  hail  and  bird  strike  directly, 
and  that  the  bird  impact  requirement  is  the  most  severe  from  an  overall 
strength  point  of  view.  However,  hail  can  cause  extensive  damage  to  the 
outer  plies  of  a laminated  windscreen. 

An  indicator  used  for  qualifying  windscreen  panels  against  panel 
failure  is  the  mean  unit  time  between  failures  (MTBF).  In  the  early 
sixties  this  mean  unit  time  for  the  Boeing  707  was  6,700  hours  (Reference  13). 
Reference  14  gives  failure  rates  for  the  747.  The  cumulative  results  of 
all  the  design  changes  as  of  November  1972  were  that  the  mean  time  before 
repair  (MTBR)  had  been  raised  from  a low  of  about  1,400  hrs.  in  early  1971 
for  the  original  design  to  a level  of  over  6,000  hrs.  for  the  curren- 
"balanced"  laminated  design.  Over  25  percent  of  the  removals  were  for 
causes  not  attributable  to  windshield  performance  per  se.  This  data  does 
not  separately  identify  the  stiff  main  ply  interlayer  or  the  improved  seal, 
and  the  MTBF  of  the  latest  configuration  is  undoubtedly  much  higher. 

Design  criteria  for  curved  and  flat  transport  and  helicopter  wind- 
shields for  environmental  conditions  may  best  be  illustrated  by  the  technical 
data  of  the  Sierracin  Windshield  Flight  Environment  Simulator  (WFES) . This 


AFFDL-TR-74-155 


mirror-image  wind  tunnel  set-up  (Reference  15)  permits  simultaneous  repre- 
sentative accelerated  life  testing,  within  a time  span  of  less  than  20 
minutes,  of  two  windshields  under  entirely  different  test  conditions.  These 
conditions  include  continuous  pressure  and  temperature  cycling  for  selected 
flight  profiles,  outside  temperature  range  from  -100°F  using  liquid  nitrogen 
to  150°F,  pressures  up  to  15  psi,  air  velocity  (flow  pattern)  over  windshield 
up  to  200  kt  (far  surpassing  the  maximum  convective  heat  transfer  coefficient 
of  sea  level  air  that  most  transport  windshields  can  encounter  in  flight; 
external  heat  transfer  rates  equivalent  to  800  kt  at  35,000  ft),  high  speed 
rain  simulation,  humidity  testing,  ultraviolet  testing,  continuous  windshield 
and  frame  deflection  measurements,  and  visual  and  photographic  observation. 
The  WFES  makes  it  possible  to  evaluate  the  reliability  of  current  and  newly 
upgraded  designs  in  a short  time  span,  possibly  prior  to  actual  production 
of  the  design.  Representative  flight  cycles  can  be  programmed  at  rates  of 
150  cycles  per  day. 

3.  BIRD  PROOF  DESIGN  CONSIDERATIONS 

The  loads  generated  due  to  bird  impact  at  usual  flight  speeds  are 
impressive.  A 4 pound  bird  at  an  impact  velocity  of  500  ft/sec  causes 
1500  psi  compressive  stress  over  an  assumed  area  of  20  in2  which  is 
equivalent  to  30,000  pounds  (Reference  16). 

Conventional  design  practice  has  to  consider,  and  usually  results  in 
a trade-off  between  several  influencing  factors  (Reference  17) . In  bird- 
proof  designs  for  high  speed  flight,  however,  impact  strength  along  with 
shock  absorbing  capacity  is'  predominant  and  several  of  the  conventional 
factors  may  be  neglected. 
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Therefore,  a reasonable  basic  specification  is  that  the  transparency 
should  withstand  a bird  of  a specified  weight  at  a specified  speed  (Reference 
18) . A more  detailed  requirement  in  the  civil  aviation  field  is  to  resist  a 
4 pound  bird  at  maximum  true  airspeed  in  operational  service  up  to  an  alti- 
tude of  8,000  feet  considering  the  most  adverse  temperature  conditions.  The 
Federal  Aviation  Agency  (FAA)  has,  for  several  years,  required  a specific 
level  of  bird  impact  resistance  for  crew  station  transparencies  on  Transport 
Category  aircraft.  FAA  regulation  FAR  25  requires  a windshield,  with  all 
its  plies  broken,  not  to  be  penetrated  nor  to  shatter  into  the  pilot's  face 
(Reference  19) . 

The  frequency  of  bird  strikes  on  707  and  720  jet  liners  has  been  less 
than  one  per  50,000  flight  hours.  None  of  these  produced  penetrations  af- 
fecting flight  safety.  Bird  penetration  had  not  occurred  once  in  over  500,000 
flight  hrs  in  1961  (Reference  20).  Reference  20  does  not  report  how  often 
cockpit  windows  were  impacted  by  birds;  however,  it  does  reflect,  that  bird 
strike  is  no  problem  for  the  707  and  720  windshields,  and  that  there  is  indeed 
a bird-proof  category  of  airplanes  flying  today.  Reference  17  states  that 
windshields  can  be  made  bird  resistant  by  adding  vinyl  interlayers  to  the 
windshield  assembly  and  heating  the  assembly  to  reduce  its  brittleness. 

Also,  windshield  panes  and  their  supporting;  structures  directly  in  front 
of  the  pilots  should  be  designed  to  withstand  bird  impact  requirements 
generally  specified  in  MIL-A-008860A  (USAF) . 

Results  of  investigations  lend  themselves  directly  to  Specification 
MIL-A-008865A  (USAF)  for  windshield/canopies  and  support  structure  and 
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to  AFSC  DH2-.1  Design  Handbook  revisions.  A follow-on  effort  ("Windshield 
Bird  Strike  Structure  Criteria  - Phase  II",  AFFDL-FB,  Sept  1973)  has  been 
proposed  to  review  the  criteria  developed  under  Phase  I (Reference  21)  to 
recommend  new  bird  strike  data  collection,  windshield  structural  proof 
tests,  compliance  tests,  and  flight  limitation  restrictions. 

4.  USAF  REQUIREMENTS 

The  USAF  probability  of  bird  strikes  is  much  lower  than  that  for 
commercial  aircraft  (1965-1972  ICAO  data),  being  .055/1000  flight  hrs. 
Extrapolation  of  data  shows  356  bird  strikes  annually  on  the  average  with 
up  to  1/3  of  these  to  be  expected  on  crew  enclosures. 

A four  pound  bird  as  a criteria  for  impact  resistant  rating  has  beer, 
found  to  be  statistically  sound.  It  is  reported  that  less  than  5%  of  the 
strikes  on  all  airplanes  involve  birds  heavier  than  four  pounds  (Reference 
20) , Most  USAF  bases  are  located  inland  while  many  large  commercial  air- 
ports are  near  large  bodies  of  water  and  major  cities.  This  may  account 
for  the  great  number  of  sea  gulls  and  pigeons  struck  by  commercial  air- 
liners . The  average  USAF  flight  may  also  be  of  longer  duration  than  the 
average  airline  flight  of  approximately  75  min.  and  above  the  main  bird 
flight  altitudes. 

For  an  update  of  criteria,  Reference  21  was  prepared  considering  the 
possibility  that  future  USAF  aircraft  may  require  windshields  capable  of 
resisting  the  impact  of  a bird  while  flying  at  velocities  up  to  Mach  1.2 
at  sea  level. 
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The  minimum  requirement  for  the  F-lll  Improved  Windshield  Program  is 
a 4 pound  bird  impact  at  500  kt  true  air  speed  (KTAS).  The  F-lll  Trans- 
parencies Development  Plan  of  Project  1926  requires  proof  testing  for  four 
(4)  aircraft  life  times  of  6000  flight  hours  utilizing  cyclic  thermal, 

pressure,  and  mechanical  loading  envelopes.  j 

The  mission  envelope  for  the  B-l  compares  with  the  one  of  the  F-lll 
i for  low  level  penetration  and  dictates  the  need  for  a similar  proven  bird 

resistant  windshield  concept,  except  that  shielding  against  a 4 pound 

i 

bird  at  650  inph  (565  KTAS)  is  required.  The  design  of  the  windshield  on 

the  B-l  is  based  on  transparency  state-of-the-art  as  it  existed  in  1970, 

- . ! 

when  bird  impact  technology  was  not  nearly  as  well  developed  as  it  is  now,  j 

and  consequently  was  not  incorporated  in  the  B-l  Airframe  System.  1 


l 

\ 

| 
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SECTION  II 
APPROACH  TO  PROBLEM 

The  methodology  has  to  meet  the  following  objectives:  (1)  prove 

the  feasibility  of  bird-proof  design  with  existing  high  strength  glazing 
materials  of  laminated  construction  in  combination  with  novel  window  frame 
concepts  capable  of  absorbing  part  of  the  shock  energy;  (2)  prove  the 
validity  of  the  mathematical  model  developed  for  determination  of  fail-safe 
configurations  by  comparing  estimated  results  and  actual  test  data  of  rela- 
ted efforts  to  optimize  for  essential  parameters;  and  (3)  determine  bene- 
ficial applications  of  conceived  design  concepts  and  existing  material 
configurations  co  suit  adverse  bird  strike  conditions  and  other  environ*- 
mental  influences. 

As  the  manufacturer  of  a bird  proof  windshield/canopy  is  usually 
required  to  dimension  it  for  structural  limit  conditions,  additional  vari- 
ables, such  as  location  of  impact,  should  be  considered.  The  safety  margin 
should  be  based  on  dynamic  maximum  stresses  in  the  mounted  window  pane  due 
to  the  design  allowable  impact  relative  to  the  dynamic  stresses  resulting 
from  an  ultimate  impulse  loading  causing  catastrophic  failure  of  the  window/ 
mounting  frame  system.  The  problem  is  to  understand  the  functional  relation- 
ship between  stresses  and  impulse  for  several  particular  design  details. 

A first  approximation  of  the  required  pane  thickness  can  be  calculated 
from  the  Grashof  equation  for  rectangular  plates  (Reference  10)  with  deflec- 
tion which  is  small  as  compared  to  the  dimensions  of  sides,  This  equation 
relates  the  dimensional  parameters,  pressure  differential,  and  allowable 
stress  for  a plate  with  unrestrained  edges. 
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The  structural  load  experienced  by  an  impacted  windshield  may  be 
estimated  by  assuming  conditions  of  1.2  Mach  at  500  ft  "ltitude  (1336  ft/ 
sec)  and  a bird  dimension  of  8 inches.  Impact  time  is  .0005  sec. , accel- 
eration is  1,336,000  ft/sec^,  and  one  half  of  the  flight  speed  is  taken 
as  the  final  speed  of  the  bird.  The  force  of  impact  is  166,000  pouncs  and 
the  mass  is  0.124  pound-sec2/ft.  The  normal  force  acting  on  the  windshield, 
inclined  at  24°,  is  67,500  pounds.  Assuming  this  load  equally  distributed 
over  a window  area  of  600  in2,  the  instant  transverse  pressure,  p,  is  113  psl. 

"Quasi-static"  solution  for  a plate  of  dimensions,  s,  and  thickness,  t, 
with  all  edges  fixed  and  uniform  load  over  the  entire  surface  (Reference  22) 
gives  the  following  maximum  stress , O m3y , at  mid-edge  and  maximum  elastic 
deflection,  Wmax  , at  the  center  of  the  plate: 


O 

max 


V 

max 


.014 


(1) 

(2) 


For  one  case  calculated,  deflection  was  wmax=,77  in,  and  the  material 
must  be  stronger  than  angxm9.35  ksi.  The  deflections  and  stresses  would  be 
different  if  mass  inertia  effects  were  taken  into  account. 

Another  approach  to  the  windshield  bird  strike  is  the  classical  vibration 
problem.  The  plate  is  hit  for  an  instance  of  time,  then  unloaded,  and  vi- 
brates at  an  undampened  frequency.  The  amplitude  of  vibration  changes  due 
to  the  severity  of  impact  impulse.  An  estimation  for  the  frequency,  f,  is 
in  Reference  23. 
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The  first  mode  of  free  vibration  of  a mounted  square  plate  of  specific 
density, p , is: 


where, 


35.99 


D = 


F.*  t 

12(1-^) 


(3) 

(4) 


f = 


(0 

2TT 


(5) 


A sample  calculation  within  the  parametrical  range  of  interest  (see 
Section  III-4)  shows  that  the  time  for  one  cycle  of  vibration  varies  between 
.COS  and  .015  sec.  This  indicates  that  a high  Speed  impulse,  typically  las- 


ting 0.0005  sec,  is  much  shorter  than  the  vibrational  time. 

The  shock  absorbing  window  pane/mounting  frame  (Figure  1)  represents 
another  aspect  of  bird  strike  which  accounts  for  the  total  energy  received 
due  to  impact.  At  the  end  of  impact,  this  energy  is  almost  entirely  kinetic 
energy  of  the  system,  and  practically  no  energy  has  been  converted  into  po- 
tential or  absorbed  energy.  The  sum  of  these  three  energies  (kinetic,  po- 
tential and  absorbed)  is  constant,  that  is 


Eimp  = Ekin  + Epot  + Eabs  = const*  (*>) 


The  kinetic  energy  of  the  one-degree-of-freedom  transverse  movement  of  trans- 
lational mass  elements  dominates  and  can  be  taken  to  represent  the  kinetic 
energy,  E]^  , in  Equation  6.  The  potential  energy,  E_0{.  , consists  of  the 
elastic  bending  energy  in  the  pane  and  frame  in  the  longitudinal  and  lateral 
directions,  and  the  energy  in  the  translational  and  rotational  springs  along 
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the  boundary  of  the  pane  In  the  window  mount.  The  absorbed  energy,  E^g, 
is  composed  of  the  viscous  material  behavior  of  the  pane  and  the  viscous 
effectiveness  of  the  shock  absorbing  frame. 

Reference  24  indicates  that  the  need  to  streamline  the  aircraft  to 
minimize  aerodynamic  dreg  is  in  direct  conflict  with  the  requirement  for 
good  pilot  vision.  The  existing  windshield  for  the  F-lll  8 aircraft  has  a 
21.5°  baseline.  A 30°  windshield  was  proposed  providing  improved  visibility. 
The  modification  would  have  increased  the  susceptibility  to  destructive 
bird  strike  energy  by  a factor  of  (sin  30°/sin  21.5°)^  »=  1.858. 

Designed  shock-absorption  for  the  purpose  of  stress  variation  is  tech- 
nically feasible  (Reference  25) . Time  and  distance  can  be  varied  to  stop 
a moving  mass.  The  resultant  stress  can  be  varied  by  tailoring  the  shock 
absorbers  to  fit  a particular  requirement.  Active  shock  isolation  has  been 
defined  practically  (Reference  26) . This  seems  to  be  more  economical  for 
continuous  requirements.  Substantial  R&D  effort  should  go  into  flight 
hardware  integration  of  a particular  proven  concept. 

The  original  bird  proof  designs  consisted  primarily  of  laminated  tempered 
glass  with  polyvinyl-butyral  (PVB)  interlayers  (Reference  16).  Bird  Impact 
testing  of  laminated  windshields  made  with  PVB  plastic  has  shown  that  impact 
resistance  is  primarily  a function  of  the  interlayer  thickness,  the  percent 
plasticizer  content  in  the  plastic  interlayer,  the  flexibility  of  the  sup- 
porting structure,  and  the  temperature  of  the  plastic  interlayer. 

Two  techniques  exist  for  achieving  bird  proof  designs  at  higher  flight 
velocities.  One  technique  relates  to  the  optimally  designed  mounting  frame 
of  Figure  1 with  capacity  to  accommodate  a sufficient  part  of  the  shock 


15 


"v 


AFFDL-TR-74-1 55 


energy  thereby  reducing  dynamic  peak  bending  stresses  in  the  pane  shortly 
after  impact.  The  second  uses  the  laminated  glazing  concept  (Figure  2) 
with  cushioning,  transparent  rubber-like  interlayer  materials  which  are 
stronger  than  PVB  under  identical  conditions  and  available  at  maximum 
allowable  thickness  for  optical  qualification.  Separate  and  combined  effects 
of  these  concepts  were  investigated. 


1.  PRINCIPLE  OF  SHOCK  ABSORPTION 

The  principle  of  structural  shock-absorption  due  to  Impact  loading  may 
be  explained  as  follows.  The  sudden  application  of  impact  force  causes 
the  window  mounting  to  vibrate  in  a transient  manner  and  in  many  natural 
modes  simultaneously  (Reference  27).  For  a qualitative  understanding,  a 
simplified  analysis  is  given  herein  where  the  masses  are  considered  to  be 
concentrated  and  the  flexible  members  to  be  massless  (Figure  3) , The 
equations  of  motions  of  such  a system  are: 


m_w  xw  + kw  (xw  - xf ) = F(t) 

raj  xf  + cf  xf  + kf  Xf  - 1^  (xw  - xf)  = 0 (7); 

(_t“0:  xw=0,  xw=0,  Xf=0,  Xf=0) 
where  the  effective  mass  of  the  window 


«V=V<4  (8) 

is  defined  by  the  bending  stiffness  and  the  first  mode  of  vibration.  These 
equations  may  be  solved  by  one  of  several  methods.  An  accurate  solution  is 
only  useful  for  a special  case.  For  our  purpose,  a discussion  of  the  nature 
of  the  solution  is  appropriate. 

Bird  impact  is  adequately  represented  by  an  impulse,  Q,  which  may  ;have 
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ATTACHMENT 

BOLT 

LOCATION ' 


r SHATTERED 
GLASS  PLIES 


VINYL 

PLY 


(REFERENCE  16) 


VINYL  PLY  ABSORBS  IMPACT 
ENERGYBY  STRETCHING  ANO 
DEFORMING.  SHATTERED 
GLASS  PLIES  REMAIN 
ADHERING  TO  VINYL. 


NOTE:  DATA  SUPPLIED  BY  LOCKHEED. 


Figure  2.  Manner  in  which  a Bird  Proof  Windshield  Panel 
Resists  a Bird  Strike 
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Figure  3.  Mass-Spring-Damper  System: 
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the  force  varying  as  a half  sine  wave,  i.e., 

Q * 
o 

F(t)  = F0*sin<j  t (0  <ti|)  (9) 

F(t)  = 0 (t  at  jo 

The  duration  of  this  load  is  on  the  order  of  .0005  sec  at  Mach  1.2. 

If  this  is  very  short  compared  to  the  natural  periods  of  motion  of  this  sys- 
tem, the  impulse  will  cause  a displacement,  xw,  and  the  frame  mass  is  res- 
trained by  the  shock  isolation  spring.  A soft  spring  allows  motion  of  the 
frame  mass  and  this  in  turn  acts  to  reduce  the  stresses  in  the  plate. 

In  the  limiting  case  of  an  extremely  soft  isolation  spring,  the  plate 
with  its  frame  acts  as  if  its  edges  were  free.  This  is  the  optimum  condi- 
tion to  resist  breaking  due  to  bird  impact  because  the  plate  is  allowed  to 
move.  However,  a maximum  displacement  is  technically  allowable.  In  the 
practical  case,  the  flexibility  of  support  consists  or  a mixture  of  torsional 
and  transverse  flexibilities, 

Two  types  of  transient  stresses  must  be  considered  during  impact. 

These  are  bending  stresses  due  to  the  relative  displacement  (;cw  - x^) , 
and  dynamic  bending  stresses  due  to  the  inertia  of  the  pane  and  the  difference 
in  acceleration  of  the  pane  center  and  the  window  support,  (xw  - xf).  These 
direct  stresses  are  to  be  superposed  on  initial  stresses  due  to  aerodynamic 
pressure  and  thermal  expansion.  The  initial  stresses  reflect  the  effects  of 
aerodynamic  contour,  pressure  differential,  and  thermal  shock.  It  is  anti- 


u> 


F(t)  dt 
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cipated  that  the  pane,  sized  for  impact  resistance,  will  also  take  the 
stresses  due  to  the  flight  environment. 

The  resonance  conditions  of  the  undamped  system  of  Figure  3 are  given 


by: 


^1,2  = 


1 + -kfJL  + 


V»f 


M 

k V 

U)  = ( nT-^l  2>* 
1,2  mw 


Vw 

kvraf 


a = 

wl  l-x 


1 m£ 

— a a = ( i + _± £_  \ ) a 

fl  ; «2  1“  • >V  *2  f2 


(10) 


This  type  of  dynamic  system  vibrates  in  one  of  the  modes  (60^  or{J2) 

with  deflections  wwl  j2=awl,  2* sin  (wl,2t+^^  a°d  Mf  l,2=afl ,2’ sin^l,2t+^ 
in  phase,  yielding  0=0  (Reference  28).  With  n^,  m£,  and  positive 
and  real,  a three-parameter  analysis  with  n^/m^  ranging  between  10-  and 
10^,  kw/k£=l  thru  10^,  and  1^/11^=  10  thru  10^  was  conducted.  The 
eigenvalues  (^  and  A^) , circular  frequencies  2 (sec”l)  , and  normalized 
amplitudes  in  the  form,  -a^ajrj^  and  +af/aw|2  are  shown  in  Figures  4,  5 
and  6.  Assuming  the  softest  f i anie  support-  shown  in  Figure  6,  ky/k^lODG, 
the  amplitudes  of  the  frame,  a£^  and  a£2,  are  the  largest  with  increasing 
mass  ratios  n^/m^,  which  is  equivalent  to  decreasing  frame  masses,  m£. 

A direct  solution  of  the  mass-spring-damper  system  (Figure  3)  was  obtained 
using  the  analog  computer  for  different  mounting  conditions  of  PPG  Industries 
flat  26"  x 26"  panel  test,  WT-18,  performed  under  AFFDL  Contract  F33615-73-C- 
3099.  This  study  indicated  that  substantial  reductions  of  impact  stresses  in 
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Figure  4.  Eigenvalues  of  Quadratic 
Eigenvalue  Problem 
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Figure  6.  Normalised  Amplitudes  of 


Quadratic  Eigenvalue  Problem 
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the  pane  could  be  achieved.  Optimization  was  performed  by  trial  and  error. 

The  original  test  deflection  data  (Figure  7,  lower  diagram,  a photo- 
mechanical record  of  the  quasi-rigidly  mounted  pane)  were  curve-fitted  for 
the  following  three  parameters: 

(1)  impulse  per  unit  mass  of  window,  (in/sec) 

(2)  natural  frequency,  -v/Ky/My  (rad/sec), 

(3)  damping  factor,  Cy/My  (rad/sec). 

F0t^/Mw  was  simulated  as  a square  wave  for  convenience  in  the  calculations. 

Since  a constant  frequency  model  was  used  to  simulate  the  actual  data 
curve,  which  is  non-linear  shortly  after  impact,  two  calibrations  were  ap- 
plied. Calibration  "A"  (Equations  13)  was  to  get  both  curves  in  agreement 
past  the  second  amplitude,  which  resulted  in  an  initiation  of  the  simulated 
time  history.  Calibration  "B"  (Equations  14)  averaged  the  first  four  ampli- 
tudes such  that  the  origins  of  both  time  histories  and  the  fifth  amplitudes 
were  congruent. 


Calibration  (A) : 
FQt]/Mw  = l,0G0  in/sec 
“ 410  rad/sec 
CW/MW  = 68  rad/sec 

Calibration  (B) ; 
F0tj/Mw  = 1,200  in/sec 
•\/Kw/My  = 467  rad/sec 
Cy/My  * 76  rad/sec 


-3 

(tj~10  sec) 

(13) 


-3 

(tj'vlO  sec) 

(14) 
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Then,  the  original  analog  circuit  for  the  rigidly  mounted  window  pane 
was  modified  to  represent  a dampened,  flexible  mount  for  the  pane.  This 
modification  added  three  additional  parameters: 

(A)  mounting  frame  - window  mass  ratio,  M^/Mp, 

(5)  stiffness  ratio,  Kp/Ky,  and 

(6)  dampening  ratio,  Cp/Cy. 

Results  for  panel  test  WT-18  are  discussed  in  Section  IV-lb. 

2.  SHOCK  ABSORBING  FRAME 

One  configuration  of  mounting  frame  for  damping  the  impact  load  is 
composed  of  a two-part  hollow  sliding  frame  containing  viscous  fluid  that 
is  released  thru  orifices  under  impact.  Steel  springs  in  the  frame  provide 
deceleration  of  the  impacted  window  and  spring-back  to  the  original  position. 
The  frame  contains  distributed  orifices  which  open  to  the. outside.  They 
dissipate  energy  by  controlled  release  of  fluid  (Figure  8)  under  impact 
pressures.  Other  damping  devices,  such  as  an  expendable  low  density  cellular 
crushing  material  developed  for  the  Apollo  crew  couch  impact  strut  assembly 
(Reference  6) , may  be  feasible. 

When  a viscous  fluid  system  is  used,  the  fluid  has  to  remain  functional 
in  the  operational  temperature  regime.  The  system  must  provide  sealing 
characteristics  as  required  for  cabin  pressurization  and  against  atmospheric 
moisture.  A fluid  could  reduce  maintenance  of  the  system.  Replacement  of 
the  fluid  after  expenditure  is  possible  with  auxiliary  pumping  equipment, 
even  uhaer  flight  conditions.  There  should  be  no  adverse  interference  of 
the  fluid  with  the  mechanical  springs. 
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Figure  8 


'MAX 


Principle  of  Viscous  Damping  System 
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Impact  activated  expenditure  of  the  fluid  to  the  outside  should 
eliminate  any  danger  of  spillage  into  the  cabin  interior.  The  fluid  must  be 
non-corrosive-,  may  act  as  a lubricant  for  the  sliding  mechanism,  must  not 
be  a fire  hazard  or  toxic,  and  must  not  block  visibility  after  release. 
Spring-back  of  the  frame  to  its  original  position  after  impact  is  necessary 
for  reconstitution  of  aerodynamic  cleanliness  and  supersonic  flight  perfor- 
mance of  the  vehicle.  The  spring  mechanism  has  to  absorb  and  discharge  that 
part  of  the  impact  energy  required  to  overcome  the  aerodynamic  drag  on  the 
windscreen  ano  frictional  sliding  forces  during  the  spring-back  phase.  The 
frame  must  accommodate  detrimental  changes  due  to  aerodynamic  heating  or  cold 
spells  at  take-off,  as  well  as  static  and  dynamic  flight  loads,  and  different 
impact  loads  including  eccentric  impact. 

Although  a detailed  design  of  the  viscous  flow  damping  system  is  not 
presently  anticipated,  its  principle  of  performance  could  be  explained  (Figure 
8)  and  a preliminary  numerical  evaluation  (Reference  29)  conducted.  The 
pressure  drop,  Ap,  due  to  viscous  friction  in  small  diameter  capillary  tubes 
reacts  the  load  due  to  impact,  F (Equation  15).  A number  of  pipes,  N,  spread 
circumferentially  around  the  frame  should  be  capable  of  providing  necessary 
release  of  the  viscous  fluid  at  a displacement,  u„„„,  during  the  relatively 

uia.-* 

short  response  time  of  impact. 

Ap  = 4FATD2 
uf  =(l/N>(D/d)2-  u 
Ap  = A-(L/d}($f/2)-u| 


*•  (15) 

(16) 
(17) 
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X = 64/Re  (for  Re  £ 2320)  ( 18) 

Ra  •-  U£-d/v£  (19) 

F =(8uL/N>(D/d)4-gf.vf-u  (20) 

F = < u (21) 


In  the  prej iminary  design  cycle,  the  law  of  continuity  for  incompres~ 
sible  fluids  (Equation  16)  is  used.  The  application  of  Bernoulli's  equation 
for  stationary  flow  conditions  (Equation  17)  includes  the  factor  X L/d  to  take 
into  account  the  kinematic  fluid  friction  at  the  walls  of  the  capillary 
tubes.  Equation  18  defines  the  factor  A which  is  laminar  for  Re£2320. 
Combination  of  Equations  15  thru  19  yields  Equations  20  and  21,  defining 
the  viscous  damping  coefficient  (ku)  which  should  assume  the  value  10  pound- 
sec/in  for  an  impact  load  (F)  of  50,000  pounds  and  a frame  velocity  (u)  of 
5000  In/sec.  k^  =10  is  technically  feasible  (Reference  30).  Assuming  a 
diameter  ratio,  D/d=150,  with  the  equivalent  piston  diameter,  D,  and  that  of 
the  capillary,  d,  a length  of  the  capillary,  L=2  in,  a fluid  of  water-like 
density  with  ^-=.00009  pound-sec2/ln4,  and  a kinematic  viscosity  (like  that 
of  mineral  oil  at  R.T.)  Vf“l*  cent! -stokes*1. 0015  in2/sec,  the  number  of 
capillary  tubes  necessary  is  N=344. 

Another  technically  feasible  concept  of  impact  absorption  (the  hollow 
rubber  unit  type)  was  treated  theoretically  and  confirmed  by  experiment 
(Reference  31).  The  force-deflection  characteristics  of  various  shapes  and 
geometrical  dimensions  were  determined. 
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3.  PANE  VISCO- ELASTIC  INTERLAYERS 

Good  interlayers  must  have  shock-absorbing  capability,  optical 
transparency,  and  adhesive  strength  (Reference  16).  Poly-vinyl-butyral 
(PVB)  is  one  of  the  early  interlayer  materials,  and  has  become  the  standard  of 
of  comparison.  Design  details  for  interlayers  have  been  suggested  in 
Reference  32.  One  design  has  metal  reinforced  interlayer  edges  with  the 
extended  sheet  metal  bolted  to  the  windscreen  frame  with  parting  agents  at 
all  bonded  edges  and  individual  laminations  at  the  areas  of  high  local  load 
concentrations.  The  necessity  of  these  design  details  challenges  low  cost 
low  weight  effective  design.  A 1/4-inch  soft  vinyl  cushion  between  main 
plies  and  abrasion  shield  did  not  improve  the  panel  resistance  to  bird 
weight  (Reference  33) . 

The  following  approximate  formulation  for  elastic  flexural  stiffness, 

El,  for  the  laminated  window  pane  was  chosens  for  use  in  the  analysis: 


= b<  Z E4  t 

ei  + 2-  Ei  \ e?> 

(22) 

upper 

lower 

layers 

layers 

2: 

= SZ  Et  t e* 

(23) 

upper 

lower 

layers 

layers 

where  E^  = effective  Young's  modulus  of  layer  i 

t^  = thickness  of  layer  i 
e^  » distance  of  layer  i from  neutral  plane 
The  position  of  the  neutral  plane  is  determined  by  Iteration.  The  elastic 
properties,  E^,  can  be  estimated.  They  are  dependent  on  the  type  of  material, 
manufacturing  process,  rate  of  loading  (strain  rate),  and  temperature 
(Reference  34,  page  56). 
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Little  is  presently  known  about  the  shock-absorbing  behavior  of  highly 
plastic  interlayers  and  their  effectiveness  in  the  laminated  arrangement. 

More  research  has  to  be  undertaken  to  understand  the  behavior  of  the  PVB 
interlayers.  The  type  of  test  data  needed  to  simulate  are  those  of  Reference 
35. 

To  simulate  the  viscous  damping,  the  elastic  bending  moment  - slope 
relationship  (the  first  term  of  Equation  24  in  finite  difference  notation 
along  a distance  Ax  at  a location  n)  is  modified  to  include  rates  of  change 
of  deflection  slope  (second  term  of  Equation  24)  accounting  for  the  viscous 
damping  component  of  bending  moment  distribution  in  the  pane. 

Mn  = + EIn  + + km]/^xn 

Material  data  on  viscous  damping  due  to  plate  bending  dependent  on  actual 
bending  curvature  rate,  6,  were  not  found  in  the  literature.  The  damping 
coefficient,  kvg,  in  Equation  24  represents  a viscous  moment  of  inertia.  It 
is  composed  of  strain  rate  force  coefficients,  k^n,  in  the  individual  layers 
acting  at  distances,  d^ , from  a mid-layer  axis  such  that  the  moments  of  upper 
and  lower  layers  balance;  that  is, 

n=upper  layers  n=lower  layers 

*vs  = b(Ji  ktn  fcn  n + £ # (25) 

The  effective  influence  of  this  definition  of  damping  coefficient 
upon  plate  dynamics  has  been  observed  in  the  follow-on  analysis,  consider- 
ing i^g  as  a bulk  parameter  constant  throughout  the  finite  difference  model 
of  plate. 
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SECTION  III 

FINITE  DIFFERENCE  ANALYSIS  METHODS 
1.  INTRODUCTORY  COMMENTS 

At  times,  generation  of  closed  form  solution  to  the  governing  differen- 
tial equations  is  not  feasible.  In  case  of  structural  problems,  the  reason 
for  this  may  either  be  the  complexity  of  the  boundary  condition  of  realistic 
structures,  or  the  necessary  description  of  nonlinear  material  behavior. 

Under  such  circumstances,  however,  an  approximate  solution  may  be  obtained 
using  numerical  techniques.  One  of  these  techniques  is  the  finite  difference 
method.  This  is  a procedure  that  transforms  a continous  system,  described 
by  differential  equations,  into  a discrete  one  consisting  of  algebraic 
equations.  The  Independent  variables  are  divided  into  a number  of  segments 
of  suitable  length  and  the  differential  operators  are  replaced  by  suitable 
finite  difference  operators,  expressed  in  terms  of  the  displacements  at 
pivotal  points.  The  choice  of  the  number  of  pivotal  points  and  the  order  of 
the  operators  affects  the  accuracy  of  the  solution. 

The  following  three  ways  are  available  to  improve  the  accuracy: 

(1)  Increasing  the  number  of  divisions  improves  the  quality  of  results; 
however,  with  larger  number  of  divisions  inaccuracies  due  to  round-off  errors 
tend  to  be  significant.  In  addition,  computational  time  increases  rapidly 
with  the  number  of  divisions. 

(2)  Accuracy  of  the  results  can  also  be  improved  by  using  higher  order 
difference  operators.  This  may  result  in  unwieldy  equations  making  the 
numerical  operations  tedious. 
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(3)  Another  way  to  obtain  better  accuracy  is  to  use  Richardson's 
extrapolation  technique  (Reference  36)  which  consists  of  determination  of 
the  displacements,  etc.  for  three  different  values  of  the  number  of  divisions 
and  extrapolation  of  the  results.  This  technique  presupposes  that  the  suc- 
cessive values  of  calculated  displacements,  etc.  approach  monotonically  to 
the  correct  solution  from  one  side. 

In  dynamic  problems,  the  solution  procedure  is  started  by  assuming  all 
time  dependent  variables  to  be  zero.  Then,  loads  are  applied  and  accelera- 
tions are  calculated.  These  are  numerically  integrated  to  give  local  dis- 
p’  ■■cements.  The  stresses  in  the  structure  due  to  the  displacements  are  used 
to  compute  new  loads.  The  procedure  is  then  repeated  continuously.  Transient 
and  stationary  dynamic  as  well  as  static  problems  may  be  solved  this  way. 

In  the  static  case,  the  load  has  to  be  slowly  applied  to  keep  accelerations 
causing  vibrations  at  a minimum. 

The  very  efficient  marching  solution  procedure  used  herein  does  not 
require  matrix  operations.  Calculations  are  straight  forward  and  results  are 
achieved  in  split  seconds  computer  time.  Accuracy  of  results  Increases  with 
smaller  increments  of  time  and  finite  difference  divisions. 

Another  advantage  for  using  the  finite  difference  methods,  as  described 
in  the  sub-sections  below,  has  been  found.  One  ready-to-use  finite  element 
analysis  routine,  NASTRAN,  for  determining  the  transient  dynamic  response  was 
found  to  use  10  to  100  times  more  computer  time  per  run  than  the  finite  dif- 
ference analysis.  It  may  be  pointed  out  that  results  obtained  with  finite 
elements  are  more  precise  than  with  finite  difference  technique.  However, 
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high  accuracy  is  not  required  at  the  preliminary  design  state.  Dozens  of 
exploratory  computer  runs  are  required  for  initial  understanding  of  the 
unknown  dynamic  behavior  of  the  new  structure.  In  addition,  little  is  known 
at  this  time  about  the  proper  choice  of  sophisticated  finite  elements  from 
the  many  available  to  achieve  the  accuracy  that  suits  the  realistic  structure. 

Based  on  the  objectives,  preliminary  dynamic  impact  load  design  methods 
for  bird  proof,  shock  absorbing,  high  visibility  windshield/mounting  frame 
systems  of  advanced  AF  vehicles  were  developed.  These  methods  are  both 
general  and  systems  oriented,  employing  realistic  material  and  structural 
data  of  the  window/mounting  frame  systems  as  the  variable  input.  Bending  of 
the  two-dimensional  structure  was  analyzed  using  elementary  cross  beam  arrange- 
ments with  a variety  of  window  mounting  conditions.  Different  cases  and 
end  attachments  were  considered  in  the  beam  element  equations. 

The  static  straight  beam  model  has  been  expanded  into  a dynamic  (Reference 
37,  Chap. 8,  Dynamic  Analysis),  two-dimensional  model  with  structural  curvature 
added  in  one  dimension  (References  37  and  38) . Computer  routines  were  developed 
and  analyses  conducted  for: 

(a)  the  flat  plate  (square  or  rectangular) , 

(b)  the  ideal  conical  shell  segment, 

(c)  the  windshield  portion  of  a F-lll  type  crew  station,  employing 
straight  and  curved  finite  difference  equations  as  derived  from  the  exact 
representative  differential  equations  for  beams,  and  selectively  dimensioned. 

Relevant  information  for  the  use  of  these  routines  is  provided  in  Appen- 
dix III.  The  routines  are  listed  in  Appendix  IV.  CDC  6600  Fortran  Extended 
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Version  3.0  was  used  as  the  computer  language.  Memory  size  required  is  50K  to 
100K  octal.  The  programs,  as  available  on  cards,  use  a card  reader  and  a list 
device.  Evaluated  results  of  these  routines  are  reported  in  Section  IV. 


2,  FINITE  DIFFERENCE  FORMULATION  OF  DYNAMIC  BEAM  EQUATIONS 

To  illustrate  the  use  of  the  finite  difference  method  for  solving 
dynamic  problems,  consider  the  flexural  response  of  a straight  beam  under 
the  action  of  time  varying  loads.  The  differential  equation  describing  the 
flexural  motion  of  the  beam  is 

EI  _4laL_-  <,(,)  + m (26) 

dX*  dtz 

where  E,  I,  and  m are  the  Young's  modulus,  section  moment  of  inertia,  and 
mass  respectively.  In  general,  E,  I,  and  m can  vary  with  position  along 
the  beam.  Moreover,  q(t-)  is  the  externally  applied  load. 

The  formulation  of  the  problem  is  completed  by  specifying 
initial  conditions  and  boundary  conditions  at  the  support  of  the 
beam.  Letting  subscript  b refer  do  the  support  at  x = b,  the  boundary 
conditions  can  be  written  as 


(V)b  = - Kwb  • (W)b  (27) 

(M)b  = - K0b  • (0)b  (28) 
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where 

0 .JL, 

(29) 

M = El  ~~  . and 

(30) 

v --  1M-. 

dX 

(31) 

are  the  slope,  bending  moment  and  shear  force,  respectively. 

Equation  27  is  the  boundary  condition  of  an  elastic  transverse 
support,  while  Equation  28  is  the  boundary  condition  for  a rotational 
support.  In  particular,  Kwb=  °»  Kwb=  oo,  K ^ 0,  and  Ke^=  correspond 
to  free,  rigid,  pinned,  end  clamped  supports,  respectively.  Finally, 
assuming  that  the  beam  is  undeflected  at  time  t = 0 and,  that  q(0)  = 0, 
the  initial  conditions  are 

w = M=e=V=~=0  at  t = 0 (32) 

dt 

Even  though  Equation  26,  with  boundary  conditions  given  by 
Equations  27  and  28,  and  initial  conditions  specified  by  Equations  32, 
can  be  solved  by  classical  methods,  the  finite  difference  approach 
will  be  used,  which  will  prove  necessary  when  attempting  to  solve 
the  curved  windshield  problem.  In  the  finite  difference  approach, 
the  partial  differential  equation  (Equation  26)  is  replaced  by  an 
approximate  difference  equation  which  is  then  solved  numerically. 

'Hie  accuracy  of  the  solution  depends  on-  the  accuracy  of  the  finite 
difference  approximation  of  the  differential  equation. 


x 
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In  applying  Che  finite  difference  method,  it  is  convenient  to 
replace  Equation  26  by  a sequence  of  first  order  differential  equations 
which  are  then  replaced  by  finite  difference  equations.  Upon 
introducing  the  auxiliary  variables  V,  M,  ©,  a and  v.  Equation  26 
can  be  replaced  by  the  following  six  first  order  equations: 


QL.  - q ( t ) * m a 

ax 

(33.1) 

(33.2) 

M=  E«  -Sr- 

(33.3) 

q=  & 

(33.4) 

> «- 
■a  xj 

ii 

o 

(33.5) 

v - -1»— 
dr 

(33.6) 

Upon  dividing  the  beam  into  n Increments  and  the  time  variable 
into  equal  increments  of  duration,  At,  Equations  33  replaced  by 


finite  differences  are 

Vn  4-jit 

( v + • ”«.t)(AX)„ 

(34.1) 

Mn+ 1 , t 

" Mn , t = 

V"+i’l  (AX)n  + £ 

(34.2) 

®n  -f  ^ » t 

®n  - •j  » t " ;l 

(34.3) 
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l.t  -*n,t 

= ©n+5't  (Ax)n  +± 

(34.4) 

< 

1 

C 

> 

II 

+ an , ♦ 

(34.5) 

<1 

c 

* 

II 

♦ vn,t  ^ 

(34.6) 

where  subscripts  n n,  n+J$,  n+1  refer  to  distinct  locations 
of  these  variables  in  increasing  order  of  steps,  Ax/2,  along  the 
beam  (see  Figure  9).  Equations  34.5  and  34.6  represent  the  Euler 
backward  integration  procedure  for  solving  Equations  33.5  and  33.6. 

It  should  be  noted  that  the  finite  difference  operators  used 
in  Equations  33.1  thru  33.4  are  of  the  second  order  type  with  symmetrical 
placement  of  the  pivotal  locations  at  n-1/2,  n,  n+1/2,  and  n,  n+1/2,  n+1, 
respectively.  Equations  of  this  type,  yielding  an  error  proportional 
to  (Ax) 2,  could  be  replaced  for  higher  accuracy  by  the  fourth  or  sixth 
order  first  - difference  operators  lD(y0)]  of  Table  9.27  of  Reference  28. 

Boundary  conditions,  Equations  27  and  28,  are  not  sufficient  in 
the  finite  difference  solution  scheme.  They  must  be  supplemented  by 
additional  conditions.  Mac  Neal  (Reference  37)  indicates  that.assumed 
edge  distributions  one-half  finite  difference  interval, Ax/2,  beyond 
the  ends-  of  the  beam  should  be  used.  The  two  required  edge  distributions, 
based  on  the  assumption  of  linear  change  of  variables  beyond  the  beam 
boundaries,  are 

Vb+^  .1  = 2 Vb,»  " (35.1) 


n * 
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0U  , = 20w  -Ou  , 

b+ ^ » t b,  t b - j • t 


(35.2) 


Combinati.  of  Equations  28,  34.3  and  35.2  yields  the  slope  at  edge  b: 


0b+i 


Kat  (AX); 


The  assumed  linear  variation  in  Equations  35  could  be  replaced  by 
higher  order  polynomials  using  additional  values  of  the  variables 
at  b + 1.5,  b + 2.5,  etc.  for  the  left-hand  and  b - 1.5,  b - 2.5, 
etc.  for  the  right-hand  boundary  of  the  beam  (Figure  9). 

An  example  solution  for  a beam  of  cell  size,  n,  demonstrates 
the  computation  procedure.  For  simplicity,  Ax,  EX,  and  m are  assumed 
constant  at  any  pivotal  location.  The  following  numbers  and  types  of 
finite  difference  equations  exist  for  the  solution:  n+1  each  of  Equations 
34.1  and  34.3,  n each  of  Equations  34.2  and  34.4,  and  two  (2)  each  of 
Equations  27,  28,  35.1  and  35.2,  for  the  left  and  right  boundary  of  the 
beam.  These  algebraic  equations  may  be  assembled  in  matrix  notation 
in  the  format: 


i r i<»> 


EI 

m (Ax)4 


with  [A]  being  a symmetric  n + 1 by  n + 1 coefficient  matrix. 
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Ac  the  end  of  time  step  N = 1:  |q(l)|n+2  * 0>  | w<l) | n+l  a 0 anci  t*'e 
calculated  accelerations  |w(l)|n+^  + 0.  Repeated  application  of 
n+l  equations  each  (Equations  34.5  and  34.6),  with  a = w,  yield 
new  values  j w(N  + 3-^  | n-f-i  t0  be  inserted  again  into  the  computational 
scheme  of  Equation  37. 

The  following  computational  approach  via  the  analog  computer 
represents  a continuous  solution  in  time  of  the  dynamic  structural 
problem.  The  finite  difference  equations  for  the  straight  beam  shown 
in  Figure  9 may  be  analoged  by  an  electrical  wiring  diagram  (Figure  10) 
representing  an  actual  computer  network  which  will  give  the  desired  solu- 
tions. Limited  computer  hardware  available  in  the  Analog  and  Hybrid 
Computer  Facility  at  WPAFB  is  a restriction  on  this  method.  This  is 
indicated  by  the  number  of  analog  components  necessary  for  the  fairly 
simple  beam  finite  difference  model  for  impulse  loading  under  symme- 
trical conditions  (Figure  10). 

3.  CURVED  BEAM  EQUATIONS 

The  curved  beam  differential/difference  equations  are  needed  for 
modeling  the  curved  geometry  of  the  windshield  in  one  or  two  directions. 

In  case  of  the  F-lll  windshield  of  single  curvature,  the  curved  beam  equa- 
tions are  applied  only  in  one  direction.  The  finite  difference  analysis 
procedure  is  similar  to  the  one  for  the  straight  beair.  However,  the 
static  equations  were  not  generalized  to  include  the  dynamic  term, 
(d^w/dt^) •m(x) , as  in  Equation  26.  The  reason  for  this  is  that  the 


\ 
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finite  difference  mass  inertia  of  the  flat  plate  and  single  curved  plate 
(Sections  III-3  and  III-4,  respectively)  is  applied  only  once  at  the 
intersections  of  the  cross  beam  arrangement. 

The  static  equation  for  the  curved  beam  comparable  to  the  one  for 
the  straight  beam  may  be  written  either  as 


# r - El  (<\  , d2“n  ) 

,n  F’  ■ ,»  v ^ ) 

d<*n  El  ( d**l  d*“t  \ 

" ‘ r3  \ d*6  d*4  / 

(38) 

dF.  «- 

(39.  D 

dF, 

d<£  = "Fn"q» 

(39.2) 

m,-rr 

n 

(39 .3) 

de  _ Mr 
dtf>  ” El 

(39.4) 

dun 

if" 

(39.5) 

du. 

T$-~"n 

(39.6) 

(F»)»+r(F»)„  =[(F,)„+i-H+J(2*ln  ¥%*i 

(40.1) 

(F.)„+i-(F.)n-r[(-F4-(l.)J(211"  ^)„ 

(40.2) 

(40.3) 
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(®)n  + | " (®)n  =(  El  )n+-j 

(40.4) 

("n)„+i  -M„.£  =[(“t)„  +(®r)„]  (2 

(40.5) 

(“t)n  + | -(ut)„  =(-un  )n+i(2sin  2 

(40.6) 

The  variables  in  the  curved  beam  equations  (Figure  11)  are: 

(1,)  dependent  variables , 

normal  force,  Fn  (pound) 
tangential  force,  Ft  (pound) 
bending  moment,  M (pound-inch) 

(2)  external  loads, 

normal  load,  qn  (pound/rad) 
tangential  load,  qt  (pound/rad) 

(applied  along  centerline  of  ring  of  radius,  r (in), 
per  angular  increment, A 0 (rad)  ) 

(3)  local  dependent  variables, 
normal  displacement , u^  (in) 
tangential  displacement,  ut  (in): 
slope  of  deflection,  9 (rad) 

These,  the  basic  equations,  and  the  corresponding  network  analogy 
for  rings  of  constant  or  variable  radii  of  curvature,  loaded  in-plane, 
are  given  in  Reference  38.  Although  expressions  for  forces  and  deforma- 
tions in  circular  rings  are  widely  publicized  (Reference  39) , the  validity 
of  the  curved1  beam  analysis  to  be  set  up  was  initially  verified.  The 
finite  difference  Equations  39  were  checked  by  comparing  them  with  the 
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Figure  11.  Variables  of  the  Curved  Beam  Equations 
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explicit  Equations  41  for  the  special  case  of  a circular  ring  subjected 
to  a pair  of  concentrated  radial  loads  in  equilibrium  (Reference  37 , 
Chap.  4 "Analogies  for  One-Dimensional  Elements"),  for  qn=qt=0: 


Fn  = P sin 

(41.1) 

F|  = Pcos  <f> 

(41.2) 

M = Pr  (cos  $■--£-) 

(41.3) 

h = sin  * + * cos*  ~ ^ 

(41.4) 

ut  * + Tr)*~  sin*  ) 

(41.5) 

e = -|*-(sm*  - w ) 

(41,6) 

Table  4.2  of  Reference  37  compares  the  numerical  results  of  the  exact 
solution  based  on  Equations  41  and  those  of  the  approximate  finite  dif- 
ference equations  for  four  (4)  cells  per  quadrant  of  ring.  The  agreement 
between  the  two  solutions  is  shown  in  Figure  12. 

The  normal  deflection,  Uj,,  at  any  point  i,  of  the  curved  beam  bent 
in  the  plane  of  initial  curvature  (Equation  41.4)  was  analytically  derived 
for  proof  of  accuracy.  The  changes  of  the  horizontal  and  vertical 
diameters,  at  ^=0  and  90°  respectively,  were  checked  against  the  solution 
derived  by  the  method  of  unit  loads  (Reference  40,  p.  94,  Eq.  (6)  )=, 

>.  * /-TT  '**  <42> 
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VARIABLES: 

fn/p 

Ft/P 

H/PR 

un/pr3/ei 

Ut/PR3/EI 

0/PR2/EI 


CIRCULAR  ORDINATE,  0 


Figure  12.  Variables  for  One  Quadrant  of  Circular  Ring  Problem 
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with  the  bending  moment,  M,  in  the  ring  section  (Equation  41.3)  due  to  the 
actual  loading,  2Q,  and  the  equation  of  a virtual  bending  moment,  m,  due  to  a 
virtual  load,  2P,  acting  normal  at  the  section  where  un=yD/2  is  to  be 
found.  Instead  of  calculating  un  as  a function  of  the  angular  ordinate,  <f>, 
use  of  the  transformation,  p'  = i +^72,  is  convenient  with  M having 
no  abrupt  changes  in  the  assumed  interval  of  Integration,  OSjS'sir.  The 
variable  of  integration,  Y' . is  introduced  in  both  the  intervals  0 
and 

M“Qr(sin  \)  (43) 


Tliis  is  not  the  case  with  the  virtual  bending  moment,  m (Figure  13) . 

At  m is  2ero  and  changes  its  sign. 

The  results  of  products  of  moments,  Mm,  and  their  integrations,,  for 
Osy's^'are 


Mm  * 2 P Q r*(-f  cos  <p‘  sin2 y + z sin  <p'sin y'-cpsY+'k cos p'sin y-^  sin  p' cosy') 
J Mm  ky  +-^co5  p'+y  sinp'-^p'  cas<p‘) 


For  p'£y'£.TTi 

Mrnm  ZPQfCk  sin  f’ cosy- J 


cos  $ sin  y'-  zsinp'siny'cosy'-t-  jr  c osp's(nly) 


If- 

\ Mm  Jy  °2PQrl(-w"f'  cosp'+^sinip'-  ^ <p'coS<t>'+%  cos  <p') 

a*  ‘ 


(44) 


(45) 


The  normal  deflection,  un  (Equation  41.4),  is  given  by 


uft- dy') 

-I?  (?£ir!  r^'f~  i) Cos Y~v-' 

“If  $'sin  rco5!f~  ii) 


(46) 
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which  is  in  agreement  with 


pr» 

U-nlf- 0) 


(47) 


taken  from  Reference  40,  chapter  8,  table  VIII  -formulas  for  changes  in 
horizontal  and  vertical  diameters  of  circular  rings  and  arches. 

In  addition  to  the  six  finite  difference  Equations  39,  six 
boundary  conditions  and  assumed  edge  distributions  are  required  for 
solution  of  the  curved  beam  problem.  Four  of  the  boundary  conditions  and 
edge  distributions  are  identical  with  those  required  for  the  straight 
beam.  Equations  27,  28,  35  do  apply,  where  V has  to  be  replaced  by  the 
normal  internal  shear,  Fn,  and  M,  o are  variables  along  the  curved  beam. 

Two  more  equations  arc  needed.  The  tangential  load,  Ft,  at  the  upper 
and  lower  boundary  was  selected  such  that  the  internal  load,  Ft,  at  the 
center  of  external  load  application,  fn,  is  zero,  that  is 


(^t)centerload  = "0  (48) 

Furthermore,  fixed  tangential  supports  were  assumed  at  the  two  boundaries, 
that  is 

(ut)b  = 0 (49) 

4.  FINITE  DIFFERENCE  FORMULATION  FOR  FLAT  PLATE 

In  finite  difference  formulation,  the  flat  plate  is  replaced  by  a 

grid  of  straight  beams  in  rectangular  x-  and  y-  directions.  The  cross 

beams  are  interconnected  by  forcing  the  local  deflections,  w„  ,,  ..  (instan- 

A >y  > L 

taneous  results  of  numerical  double-integration  of  Equations  34.5  and  34.6) 
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to  be  the  common  deflections  of  the  two  intersecting  beams.  In  this 
idealization,  torsional  rigidities  of  the  beams  are  neglected. 

Replacing  the  pivotal  indices  of  Equations  34,  such  as  n-1/2,  n,  n+1/2, 
n+1  by  combinations  of  x-1/2,  x x+1/2,  x+1  and  y-1/2,  y,  y+1/2,  y+1,  the 
left  side  of  Equation  34.1  may  be  replaced  by  the  finite  differences  of 
the  shear  forces,  V,  in  the  two  directions.  Assuming  that  Ax  = Ay  = As  * 
constant,  the  finite  difference  equations  for  the  shear  forces  are 

W^x+i,y,t  ~ hicVi,y,l  + ~ Wy)»,y-i,t  = 

where  Q and  AM  are  loads  and  masses,  respectively. 

Each  of  Equations  34.2,  34.3,  and  34.4  may  be  assumed  to  represent 


and  y-  direction.  Thus, 

•AS 

(50.2) 

(M?Wt  = " ^Vx,y+|r,t 

• AS 

(50.3) 

^x+£,y,t  ' 

AS 

(50.4) 

AS 

(50.5) 

wx+ilYlt  (e*K+i,y,t  ’AS 

(50.6) 

Wvy+i;t  - WX/y/t  = (©y)x;y+i,£  *AS 

(50.7) 
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when  as  is  a grid  parameter.  As  in  the  straight  beam  problem,  pairs  of  boun- 
dary conditions  and  distributions,  for  left  and  right  edges,  in  x-  and  y- 
directions  are  needed. 

Equations  50  could  be  written  for  variable  finite  difference  grid 
spacings,  (Ax)x>y  and  (Ay)x,y.  In  this  case,  any  rectangular  plate  could  be 
treated  with  stepped  changes  from  grid  line  to  grid  line  of  plate  proper- 
ties or  geometrical  irregularities  (as  reinforcements  or  cut-outs) . 

The  bending  stiffnesses  may  be  determined,  as  for  the  homogeneous 
plate  of  thickness,  t: 


t*  (Ay)x,v 
12 


(51.1) 


(ei 


(51.2) 


Using  Equations  50.2  and  50.3,  stresses  at  the  surfaces  of  the  pane 
are  given  by 


6'" 


6 -M 

A5‘t2 


(52) 


where  M is  the  calculated  bending  moment.  This  equation  is  valid  for 
6’-Fyieid.  For  higher  stresses,  a fictitious  stress,  Ffe,  referred  to  as 
the  bending  modulus  of  rupture  is  defined  by 

Fb=  (53) 

where  M is  the  ultimate  bending  moment  (determined  from  tests  of  similar 
beams),  I is  the  moment  of  inertia  of  the  cross  section,  and  c is  the  dis- 
tance from  the  beam  neutral  axis  to  the  extreme  fiber. 
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Throughout  the  analysis , external  shock  loads  are  introduced  at  the 
nodal  points  of  the  finite  difference  grid.  The  dynamic  leading  condition 
applied  analytically  is  a schematic  distribution  of  the  impact  load,  Qtot, 
representing  a pressure  mound  that  varies  with  time.  Methods  for  deter- 
mining Q are  given  in  References  41  and  42.  Motion  pictures  of 
experimental  tests  make  it  possible  to  estimate  the  windshield  area 
covered  by  the  traveling,  decelerating  mass  of  the  bird  (Appendix  II-l) . 

The  mass  (Equation  54)  is  distributed  similarly  to  the  stiffnesses 
(Equations  51)  between  the  beams: 

(AM>x,y= r(Ax)X,y‘ (Ay)x,y-t 

The  marching  procedure,  starting  out  at  one  edge  of  the  plate  and 
considering  one  grid  line  after  the  other,  is  as  follows: 

Calculate  all 

(1)  ev  and  ev,  Inside  of  plate  and  at  the  edges  (Equations  50.6  and  50.7), 

a y 

(2)  ox  and  ey,  at  boundaries  (Equations  28,  35.2,  50.4  and  50.5), 

(3)  Mx  and  My,  at  supported  edges  (Equation  28), 

(4)  Vx  and  Vy,  at  supported  edges  (Equation  27), 

(5)  Mj^  and  My,  inside  of  plate  and  at  the  boundaries  (Equations  50.4  and  50.5) 

(6)  Vx  and  Vy,  inside  and  at  boundaries  (Equations  50.2  and  50.3), 

(7)  w,  inside  and  at  boundaries  (Equation  50.1) , 

(8)  w,  w (Equations  34.5  and  34.6,  respectively), 

(9)  Repeat  (1)  thru  (8)  for  next  time  step,  at  time  interval,  At,  etc. 
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This  marching  procedure  is  sketched  in  Appendix  III-l  (Figure  42, 
Sequence  in  Cyclic  Loop  for  Calculating  Finite  Differences).  It  may  be 
noted,  that  the  program  variables  are  discussed  in  Appendix  III-l. 

The  following  parametrical  example  illustrates  the  use  of  the  flat 
plate  analysis  method.  Requirements  for  mounting  a flat  rectangular 
windscreen  were  determined  on  the  basis  of  a parametric  screening  of 
mounting  conditions.  A 2'  x 2'  window  pane  (Figure  14)  of  variable 
thickness  (1"  average)  having  the  stiffness  of  polycarbonate  at  200°F  and 
up  to  5 times  that  stiffness  was  assumed  to  be  impacted  by  a pyramidally 
distributed  load  impulse  (1/2  millisecond),  Qt(Jt^  100,000  pounds.  The 
flat  plate  routine  of  Appendix  IV-2  was  used. 

With  the  arbitrary  grid  size  of  6 x 6 * 36  (13  x 13  including  mid 
mode  points)  for  a square  panel  of  sides  of  2 ft,  the  grid  parameter  is 
as=*4  inch.  At  a pane  thickness  of  t*l  inch,  the  bending  moment  of  inertia, 
I*=t^*  As/12,  is  0.5  in^. 

Two  principal  design  parameters,  Qtot,  and  plate  thickness,  t,  can 
be  varied  so  that  stresses  due  to  ^ and  My  (in  this  case  +3467  pound- 
in)  , remain  below  the  assumed  shatter  strength  of  the  polycarbonate 
(Ftu=5.2  ksi  at  200°F  (Table  III)).  In  addition,  the  deflection  is  res- 
tricted to  plus  or  minus  1 to  2 inches  at  any  location.  The  effects  of 
flexible  supports  along  the  frame  boundaries  were  investigated  by  assuming 
parametric  ranges  of  technically  achievable  values  for  the  frame  stiffness 
and  the  other  boundary  parameters.  Extreme  values  occurred  at  the 
center,  mid  edge  and  corner  locations  for  the  square,  symmetrical  case, 

Figure  15  shows  the  effects  of  both  a hard  and  a soft  shock- 
absorbing window  mount  on  the  bending  moment  at  the  center  of  the  plate 
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Y 


Figure  14.  Symmetrical  Impact  Loading  of  Square  Plate 


55 


AFFDL-TR-74-1 55 


MATERIAL  PROPERTIES  ( 


polycarbonate 
RT  2 oo  °P 
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(+  5.2  ksl  (Reference  44)) 


TABLE  III 
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50,000  POUNDS 
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SPRING  CONSTANT 

10,000 

25  POUND/ IN 

ABSORBED  ENERGY 

28 
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— 
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• o5  .1 

SECONDS 
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Figure  15.  Central  Bending  Moment  for  Hard  and  Shock- 
Absorbing  Window  Mount 
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(Figure  14) . The  dashed  lines  are  the  upper  and  lower  bounds  on  shatter 
strength  indicating  that  Qtot  = 50,000  pounds  is  close  to  the  failure 
load.  Figure  16  shows  the  mid  pane  deflection,  wQ,  and  mid  frame  deflec- 
tion, W3,  for  the  shock-absorbing  case  (soft  mounting)  of  Figure  15. 

Using  a Young's  modulus,  E»0.2*10^  psi,!  representative  of  the  pane 
stiffness  for  polycarbonate  at  elevated  temperature,  the  pane  bending 

stiffness,  EI(l-V^)  'vEI  =EI ,,  in  the  longitudinal  and  transverse  direction, 

x y 

is  of  the  order  of  10^  pound-in2.  The  bending  stiffness  of  the  frame, 

EIj,  is  106  pound-in2.  Finite  difference  distribution  of  bending 
stiffness  is  shown  in  Appendix  IIL-l  (Figure  43) . 

Total  mass  and  mass  distribution  for  a pane  volume  of  24"  x 24"  x 
1"=  576  in-*  and  specific  mass  density  of  polycarbonate  (6.2  lb/ft2)  is 
(4) (6 . 2) / ( (32 .2) (12) ) = .064  slugs.  Twenty- five  mass  allocations, AM,  are 
required  at  25  nodal  points  of  the  interior,  and  24  allocations  of  A M/2  at 
the  boundary  in  case  of  the  6x6  mesh  resulting  in  AM=. 064/37  slugs. 

An  effective  frame  mass  totaling  AMf>AM/2  has  to  be  added  toAM/2  at  the 
edge  of  the  frame. 

Simulation  of  dynamic  mounting  must  be  as  realistic  as  possible. 

In  this  plate  model,  mounting  design  is  represented  by  the  four  properties, 

Ky,  Ke  and  Ke  varying  with  impact  behavior  of  the  window/frame  system 
(Figure  43  of  Appendix  III-l) . Spring  constants,  Ky,  distributed  at  the 
nodal  points  along  the  frame,  may  act  continuously,  or  as  1^(1),  while  the 
frame  is  deflected  downward,  and  as  1^(2)  during  the  upward  mode  of  deflection. 
The  latter  scheme  simulates  the  design  mentioned  above.  K^(T)  models 
effective  viscous  damping  during  the  first  quarter  cycle  (Figures  1 and  7) . 
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SECONDS 


Figure  16,  Deflections  Wg  and  W3  for  Shock-Absorbing  Design 
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After  release  of  the  damping  fluid,  less  effective  damping,  K^(2), 
takes  over.  The  clamp-in  conditions,  Kg  and  K^,  are  assumed  different 
during  downward  and  upward  deflection.  Typical  values  for  the  mounting 
system  are  ^(1)=  25  pound/in  (soft  spring),  ^(2)=1, 000,000  pound/in 
(hard  spring  back),  K^(l)=  10  pound-sec/in  (highly  efficient  viscous 
damping),  K^(2)=.5  pound-sec/in  (less  efficient  damping),  and  Kex=K0y- 
100,000  pound-in/rad  (stiff  clamp-in).  An  actual  mounting  system  may 
have  a variable  distribution  of  these  four  mounting  properties  to  achieve 
an  improved  shock  absorption. 

Initially,  the  pyramidal  impulse  load  distribution  (Figure  43.3, 
Appendix  IIT-1)  yields  unit  loads  of  Qj=Qtot/35,  and  Q2=2-Q^, 

Qtot  can  be  in  the  105  pound  range  equivalent  to  103  psi  average 
pressures  if  equally  distributed.  This  load  occurs  during  5 time  steps, 

N,  of  duration  At=  .0001  sec  each.  After  impact,  a nominal  unit  load  of 

O. 1  pound  or  any  other  constant  pressure  load  is  assumed  to  act  at  the 
previously  impacted  nodal  points.  Time  steps  at  a rate  of  less  than  10 
per  mllli  second  were  utilized  after  the  duration  of  impact  to  decrease 
the  overall  computational  effort. 

Using  the  above  input  data,  the  dependent  variables  have  values 
in  the  following  ranges: 

slopes  of  deflection,  ©x  and  ©y  +I0"2(rd)  or  +10®  (deg) 
bending  moments,  and  My  +103 (pound-in) 


shear  loads,  V„  and  V„ 

x y 

normal  accelerations,  w 


+10^ (pound) 
+.10"5(in/sec2) 


velocities  of  displacement,  +10^  - 103 (in/sec) 


deflections,  w 


+l(in) 


N 
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Energies  in  the  finite  difference  scheme  for  the  6x6  grid  size 
plate  are  calculated  as  follows.  For  convenience,  the  node  and  mid  node 
points  are  numbered  consecutively.  Thus  the  node  points  have  odd  subscripts. 
At  the  beginning  of  time  history  at  which  impact  starts,  the  three  energies 
considered  (Equation  6)  are : 


Ekin(t=0)=Epot(t=0)=Eabs(t=0)=0  (55) 

Ekin  is  newly  determined  at  every  time  step,  N: 

Ek£n=i  ' Z ‘ X (56) 

(=1,3,5-13  j»  1,3, 5 -13 

Ep0t  and  EakS  are  supplemented  at  every  time  step,  N: 

Epot(t+At)  = Epot(t)  + ^ Epot(t+  At)  (57.1) 

Eabs(t+At)  = ^abs(t)  + A^abs(t+  At)  (57.2) 

Epot  consists  of  the  following  components: 

^Epot  = ^EpotEI  i"  ^^potke  * ^Epotkw  (58) 


AE  is  the  change  of  elastic  potential  energies  in  the  x-  and  y- 

potE  L 

directions  due  to  pane  as  well  as  window  frame  bending.  AEpotke 
and  AEp0tkw  are  elastic  energies  stored  in  the  rotational  and  trans- 
verse spring  attachment  systems  around  the  perimeter  of  the  window. 


They  are  given  by: 


AEpetEKi)*  Z ' Z (Oxi+l,j«)  -0xH,ja))[(©xi+l,j(t)-©xH,jU)) 

(=3,5,7-11  1,3,5-13  - E1 

-(©x  £tl,j(-t-At)  _©x£-|,j(t-Afc))J 

+ Z • Z (0yijH(t)  -©yy-i(t))[(0ytjj*Kt)~©jri|j-i(t)) 

Is  1)3,5"*  13  Js3 

"(0y  £/j"H  (t -At)  - ©y  i,  j-Kt-At))]-^- 


(59) 
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AEpot  k®(t)  Z * Z kexj  0*Lj(l)  (©xij(i)-©xij(i-Al>) 

i=  I ami  13  J « Ij3;5«-|3 

+ Z • Z k©yt  ©ytj(4)  (Oyij(t)-Oyij(-t-At)) 

I*l;3;5—13  j*laml  |3 

^Ep.t  kw(t)  * Z ' Z kwj  Wij(t)  ( WiljCtJ  - WlJ  (t  +Al)) 

tm|a«4  13  J- 1,3,5-13 

+ Z * Z kwi  Wij(t)(  Wcj(i)  - Wij(t+At)) 

j»  | ami  |3 


(60) 


(61) 


Bending  stiffnesses  of  the  pane  or  the  frame,  El  and  El  , were 

x y1 

assumed  constant  with  time,  t.  Spring  constants,  1^,  ke  and  ke  , may 

x y 

vary  with  time.  AEgbg  consists  of  the  following  components: 


AEabs  ~ AEabskv  + AEabsk6  + AEabsk* 


(62) 


AEabskvis  the  visco-elastic  energy  in  x-and  y-  direction  absorbed  in 
the  damping  interlayers  of  the  pane  as  well  as  frame  due  to  stretch- 
damping occurring.  AEabgk^  and  AEabgk^  are  visco-elastic  energies 
absorbed  in  the  built-in  damping  system  of  the  frame.  The  energies  are 
given  by 

A Eabs kv(-t)  - Z ' Z ( ®*c+|,j(t)  " ©xM,j(l)).[(®xi+l;j(t)  “Oxt-ljjU)) 
t*3;5^7**,ll  ,1*1,3,5—13 


Ml 

AX 


+ Z z ( e>-y+t(t)  - <=>yi,j-l(t))  [(9y l;jf|(t)  -©yiji-Ki)) 

t-l,3,SH3  J«3,5,7...|' 

-(©yi,j  + l(t-At)  -©yi^-Kt-At))]-^ 

where  viscous  damping  constants  of  the  pane  or  the  frame,  kVgx  and  k^y 
were  assumed  constant  with  time. 


(63) 
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A Eabsko(i)  ' Z ' Z k 6xj  0x  ij  (t)  ( ©xij  (t)  -Qxdj  (i-At)) 

L-  I QnJ  13  j=  1,3,5-13 

+ Z * Z k&yi  0ytj(t)  ( ©y  ij(t) -Qyij(i-At)) 

i*  Ij3;5-I3  j * I Qiiel  13 

AEabsMt)  = Z ' Z (<Wj  U/tj  (1)  ( Wijtt) 

1*1  and  13  j*  1,3,5—13 

+ Z 1 Z kwj  U/£j{£)(  U/tj(i> -l^tj(t+Atj) 

I*l,3,5—I3  j=lan<i|3 

Constants,  kgx,  kay  and  k^,  can  vary  with  time. 

Stepwise  summation  yields  potential  energy  (Ep0t),  kinetic  energy 
(Ekin) , and  absorbed  energy  (Eabs)  in  the  system.  Epot,  E^q  and  EajjS 
are  2ero  before  impact.  At  the  end  of  impact,  Ep0t,  the  total  energy 
received  is  at  its  maximum  and  0,  since  the  motion  of  the  window/ 

frame  system  has  barely  started.  Several  different  computer  runs  showed 
Epot  + Ekin  + Eabs  t0  1)6  a^most  constant,  in  agreement  with  a true 
structural  system,  and  that  Ea{,s  really  picks  up  during  the  first 
quarter  cycle  using  effective  viscous  damping.  After  the  first  quarter 
cycle,  Epot  + Ej^  is  the  total  energy  remaining  in  the  window  frame 
system. 

During  incorporation  of  the  energy  equations  into  the  computer 
routine,  Ea^g  was  assumed  to  be  absorbed  by  the  frame  at  the  24  nodal  points 
at  a viscous  damping  of  K^,(l)=10  pound-sec/in.  Typical  frame  deflection  was 
.054"  and  the  time  was  .03  sec  for  the  first  quarter  cycle.  The  estimated 
result  was  EajjS  *v24  pound-in  (one  case  of  low  energy  absorption) . By 
comparison,  the  "flat  plate  routine"  resulted  in  16.8  pound-in.  EpQt  and 
E^in  results  are  similar. 
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Mean  intervals  of  oscillations  (usually  hundredths  of  seconds) 
or  frequencies  may  be  obtained  by  plotting  time  histories  of  variables. 

The  state  of  forced  vibration  is  experienced  during  the  short  time  of 
impact  and  periodic  free  vibration  thereafter.  Zero  crossings 
of  time  histories  of  bending  moments  at  the  center  of  the  plate,  at 
the  frame,  and  elsewhere  show  similar  occurrences  (although  at  different 
times)  as  observed  in  the  case  of  one  computer  run  of  a 1/4  sec  real  time 
event.  First,  the  bending  moment  changes  from  negative  to  positive 
at  .067  sec  after  Impact,  the  next  crossing  from  positive  to  negative 
is  at  .115  sec;  the  third  crossing  is  at  .1885  sec  and  the  fourth  one 
at  .2468  sec.  A full  oscillatory  cycle  occurs  within  .12  sec.  This 
is  equivalent  to  ~8  cys/sec. 

Dynamic  behavior  of  the  typical  window  pane  with  shock-absorbing 
support  structure  is  presented  in  Figures  17  thru  24  and  Tables  IV  thru 
VI.  The  following  constants  were  used  in  the  parametric  study:  size, 

s=2  ft  (square);  load  Qtot= 100, 000  pounds  (pyramidal);  time, 
sec;  mass,  m=.77  nlugs/ft2/in;  frame  stiffness,  El£=5’EI  (pane)  pound-in2. 

For  the  original  grid  size,  As*4",  and  a time  step,  At=.0001  sec, 
no  solution  instability  occurred.  After  completion  of  the  original  study, 
the  erroneous  mass  assumption  became  obvious  (with  m=. 77/12=. 064  being 
realistic,  rather  than  .77  slugs/ft2/ln) . In  addition,  a lower  Qtot= 
15,000  pounds  (pyramidal)  was  taken  to  be  more  typical  of  an  inclined 
panel  impact.  Using  the  new  m and  Qtot,  a high  percentage  of  the  para- 
metric series  of  runs  became  unstable  under  certain  conditions  of  El, 

1^,  k0,  k^,  etc.  Quantitative  stability  limits  of  these  conditions 
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Figure  21.  Maximum  Mid  Corner  Stresses 
of  Plate 
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TABLE  IV 

EDGE  DEFLECTIONS  OF  HARD  SPRING  SUPPORTED  PLATE 


E 

th 

kw 

k0 

w3inax 

w3min 

.2  lo6 

.9 

0 

1 00000 

.o33 

— 

BBI 

.9 

0 

1 00 

. oo34 

— 

1o6 

1.5 

0 

1 00000 

.0122 

-.o132 

1.5 

0 

1 00 

. ol  ol 

- . oo99 

mm 

.9 

lo 

1 00000 

,oo54 

-.oo32 

mm 

.9 

lo 

loo 

.oo33 

— 

1.5 

lo 

1 00000 

. ol  21 

-.0123 

^o5 

1.5 

lo 

1 00 

.oo94 

-.oo93 

K 

i 

i 


6 


AFFDL-TR-74-155 


TABLE  V 


EXTREME  STRESSES  AND  DEFLECTIONS 


1 00000 


<x  (time) 

i8r  !ksc 

cr~l  timeT 

for  k ~ 0 

vsc 

<r3 

-2.00  (.oil) 
_K37  _(_.  o25_). 

-1.73  (.008) 
1.88  ( .c?5) 

stresses 
neglegibly 
small  at 
k©  = loo 

-2.00 

2.o7 

(.oo7) 
( .o2o) 

1.63  ( .008) 
-1.47  (.0I8) 

-2.17  (.oo7) 
2.84  ( .o2o) 

1.74  (.008) 
-1.71  (.0I8) 

-1.67 

.89 

( .oo7) 
(.o23) 

v 

-1.77  (.oo7) 
1.19  (.c-25) 

-1.98  (.oo7) 
1.99  ( • o2o) 

1.63  (.oo7) 
-1.36  (.o17) 

-2.15  (.oo7) 
2.71  ( .o2o) 

1.72  (.oo7) 
-1.62  (.018) 

1°  -2.81  (.022) 

loo  -5.17  (.o13) 

-3-16  (.023) 

-2.57  (.olo) 

1 ooooo  -2.59  (-oil) 

1 .60  ( .o25) 

1 o -2.66  (.013) 

1nA  -2.51  (.021) 

100  ~2"758~(  ToVJT 

1o  -2.0o  (.023) 

1o  -2.53  (.olo) 

, 1.3o  (.025)  - 

. -2.54  (.oil) 

10  1.55  (.o25)  - 


comments: 
for  t=s.o25, 

o 3 every  loth 

time  interval 

1.68  (.019)  2.12  (.025)  looks  favor- 

_aole,  however 

T.To  ( .o19)  I 2. 12  (.025) 

deflections  \ 

smaller  than  \ , 

those  at  k0  = loo  \ / 

1.69  (.o25)|2.o4  (.o25)  \ / 

7.66  T-  o25T  zTo~(  ^25l  \ / 


■ 

I 


3.1o  (.008)  1.19  (.olo)  .38  (.0I8) 

iL'li  L22H  IL 

3.15  ( . 00O)  1 . 2o  (.olo)  .39  (.0I8) 
-1.53  (.025)  -.23  (.o25) 


3.o8  (.008)  1.18  (.olo)  .30  (.0I8)  trn  dampens  out 

:ll5£_L°_25i  too  high 

3.13  (.008)  1.19  (.Olo)  .39  (.018)  / tolerable 
-1.52  (.025)  -.22  (.025)  Wo  toierable 

w-  is  small 


PARAMETRIC  INVESTIGATION  OF  MOUNTING  CONDITIONS  FOR  PLATE 
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could  have  been  easily  investigated.  Instead,  At=. 00001  sec  was  arbi- 
trarily used,  and  all  parametric  runs  investigated  became  stable  again. 

The  maximum  center  deflections,  womax(in),  which  occur  first  after 
impact  were  calculated  and  are  selectively  reported  in  Figure  17  for 
both  a soft  pane  material  (E=.2*10^psi)  and  a very  stiff  one  (E=10^psi) 
for  pane  thicknesses  between  .6"  and  1.5".  The  following  mounting  para- 
meters, 1^,  k^,  and  ke,  were  assumed  equally  distributed  along  the 
boundary: 

1^=10  thru  100,000  pound/in, 
k^=  0 thru  10  pound-sec/in, 
ke=100  thru  100,000  pound-in/rad . 

The  viscous  constants,  k^,  ky^  and  kvgc£,  initially  were  assumed  zero 
to  demonstrate  the  influencing  effect  of  k&  alone.  Curves  plotted  in 
Figure  17  are  the  peak  deflections  taken  from  Figure  18. 

Maximum  deflections,  w3max,  for  the  pane  at  the  middle  of  the  rein- 
forced edge,  equaling  the  amount  of  the  local  spring  deflection,  occur 
at  time  >.025  sec  in  case  of  the  soft  spring  of  k^,  = 10  pound/in  (Figure 
19),  except  for  cases  (E) , (G)  and  (H).  1^=100,000  pound/in  causes 

negligible  maximum  deflections,  w^,  at  the  middle  of  the  edge. 

Surface  stresses  occurring  at  the  center  and  the  mid  edge  of  the 
plate  are  shown  in  Figures  20,  21,  22.  For  the  typical  ranges  considered, 
time  histories  of  kinetic  energies,  E^n,  are  shown  in  Figu.  23. 

The  viscosity  effect  of  the  interlayer  represented  by  kVsc>  has 
been  investigated  for  32  different  combinations  of  parameters  represented  by 
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1^=10  and  100,000  (pound/ln) , 
k^=0  and  10  (pound-sec/in) , 

ke=100  and  100, 000 (pound-in/ rd ) , 
k0=O  and  10  (pound-in-sec/rd) , 

kvsc=0  and  10  (pound-in  -sec/rd) . 

It  was  observed  that  use  of  ky^-lO  tends  to  flatten  peak  variations 
of  stresses  in  the  original  curves  of  k^^O.  In  addition,  curves  for 
kyg^lO  lag  behind  those  for  k^^O  (Figure  24).  The  set  of  parameter 
combinations  was  screened  for  extreme  values  of  stresses,  SQ  and  63, 
and  deflections,  wc  and  W3. 

The  effectiveness  of  shock  absorption  in  dampening  the  oscillating 
stresses  is  shown  in  Figures  20  and  21.  The  peak  stress  of  the  first 
dycle  for  the  damped  system  usually  is  slightly  lower  than  the  maximum 
stress  for  the  undamped  one. 

The  undamped  stress  reversals,  rather  than  the  one-time  stress 
peak  of  the  damped  case,  may  be  the  reason  for  fatigue  failure.  For 
the  undamped  system,  the  possibility  of  local  build-up  of  peak 
stresses  exists,  as  the  sum  of  kinetic  and  potential  energy  in  the 
system  theoretically  remains  constant.  As  time  elapses  after  impact, 
energy  dissipation  by  natural  causes,  not  accounted  for  in  this 
simplified  analysis,  could  have  occurred. 

The  largest  tolerable  deflections  during  the  time  history  after 
impact  (Figures  17,  18,  19)  are  another  design  restriction.  A trade- 
off decision  has  to  be  made  between  soft  springs  (yielding  higher  de- 
flections and  lower  pane  stresses),  and  harder  ones.  A matrix  of  criteria 


80 


AFFDL-TR-74-155 


for  design  of  the  flat  impact  resistant  windshield  has  been  prepared 
using  the  limited  number  of  results  and  is  presented  in  Table  VI.  Low 
maximum  dynamic  stresses,  acceptable  frame  and  window  maximum  deflections 
and  high  impact  energy  absorption  result,  from  using  technically 
feasible  soft  translational  springs,  hard  rotational  clamps  and  effective 
viscous  damping  units. 

The  influence  of  the  shape  of  impact  pressure  mound  (such  as  constant 
vs.  pyramidal  distribution),  the  time-load  profile  (such  as  one-minus- 
cosine  vs.  square  wave),  the  effect  of  geometrical  eccentricity  and 
severity/rapidity  of  impulse  (Qtotl' timpl  = Qtota'^lmpZ  = const) » may 
be  investigated  for  "hard,  medium  and  soft"  boundary  parameters. 

Equally  distributed  loads  cause  higher  stresses  at  clamped  edges  than 
pyramidal  loads.  Eccentricity  of  impact  may  create  higher  stresses 
(shifted  from  the  center  of  the  pane  towards  one  edge) , along  with  overall 
lower  deflections.  Rapidity  of  impact,  t^p}  < timp2  causes  higher 
stress  concentrations  for  very  short  periods  of  time.  The  advantage  of  a 
less  severe  Qtot2  < Qtotl > however,  should  be  determined. 


5.  FINITE  DIFFERENCE  IDEALIZATION  OF  WINDSHIELD  (F-lll) 
a.  The  F-lll  Windshield 

The  structural  model  for  the  F-lll  Crew  Module  considers  the  wind- 
screens to  the  left  and  right  held  together  by  the  head  beam,  with 
sills  acting  as  lateral  supports,  a bow  in  the  front  attached  to  the  bulk- 
head, and  an  unsupported  arch  in  the  rear,  which  simultaneously  acts  as 
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part  of  the  frame  of  the  hatch.  The  stiffness  and  mass  properties  of  these 
parts  and  the  edge  reinforcements  of  the  window  panes  were  estimated  from 
the  details  of  the  drawings.  The  head  beam  is  2 times  stiff er  than  the 
arch.  The  bow  is  8 times  stiffer  than  the  arch  and  the  sills  are  over  20 
times  stiffer  than  the  arch.  The  elastic  support  of  the  remaining  module 
structure  relative  to  the  windscreens  can  only  be  guessed,  and  production 
runs  show  how  important  this  input  is. 

b.  The  General  Equations 

As  in  the  case  of  finite  difference  formulation  for  the  flat  plate, 
the  F-lll  windshield  is  replaced  by  a cross  beam  arrangement  consisting 
of  curved  beams  in  x-direction  and  straight  beams  in  y-direction. 

Further,  it  is  an  applied  case  of  the  basic  monolithic  conical  shell 
segment  with  right  and  left  window  panes,  and  support  members  consisting 
of  front  bow,  right  and  left  sill,  aft  arch  and  center  head  beam. 

The  finite  difference  equations  valid  in  both  the  general  shell  and 
windshield  case,  contain  the  following  variables  in  the  radial  x-direction: 

(1)  bending  moment  (pound-in),  Mx, 

(2)  normal  deflection  (in),  un, 

(3)  tangential  internal  force  (pound),  Ft, 

(4)  slope  (rd),  9X, 

(5)  normal  force  (pound),  Fn, 

(6)  tangential  deflection  (in),  ut. 

They  are  related  by  six  (6)  first-order  curved  beam  difference  equations. 

Variables  in  the  longitudinal  y-direction,  related  by  four  (4) 
first-order  straight  beam  difference  equations,  are: 
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(1)  moment  (pound-in),  My, 

(2)  normal  deflection  (in),  un, 

(3)  slope  (rd),  0y, 

(4)  normal  force  (pound) , V. 

Tangential  force,  Ft,  is  equivalent  to  circumferential  tension,  and 
normal  force,  Fn,  to  radial  shear.  Fn  could  be  used  as  an  indicator  of  the 
shear  stress  distribution  in  the  window  pane.  V is  the  shear  force  in  the 
longitudinal  direction. 

The  following  are  the  general  finite  difference  equations: 

(F*)x+±,y,t -&)**,#*  (2sillf)x/(VW'(V\ y-ht  = 

~ A M \y  ^ ) %,y)t  (66  • 1) 


~ = ^ SOl^x+i,y,t  W>y>t 

(66.2) 

$ 

^^x+l,y,t  " = “ y,t  & r stnT^^+i,y ,t 

(66.3) 

^Vxjy+l/t  ” ” ~^K,y +£;t  ^^x,y+i 

(66.4) 

l 

(e*^,y,t  ~ (0^x-i,y,t=  (l fei,y,t  (r'A^)x,y 

(66.5) 

f 

£ 

i 

i 

(66.6) 

f 

t 

i 

j 

j 

(^x+l^t  “ ^■»*^y,t  " SirL~2^\+^>y 

(66.7) 

i 

i 

i 

^x^t/"  (UA,y,t  = ^A^x,y 

(66.8) 

! 

i 

\ 

t 

(uAc4fcy,t  r(tt*^-i,y,t*  ”^BW  5Ca^y 
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The  following  equivalencies  exist  between  the  flat  plate  Equations  50 
and  Equations  66: 

66.1  <-*•  50.1 

66.3  50.2 

66.4  <-*  50.3 

66.5  +-■»  50.4 

66.6  50.5 

66.7  <— »•  50.6 

66.8  > 50.7 


Equations  66.2  and  66.9  account  for  the  effects  of  curvature,  tangential 
load,  Ft,  and  tangential  displacement,  ut,  respectively. 

Equation  66.1  is  the  balance-of-force  equation  for  calculating  the 

I 

normal  accelerations,  un.  un  is  integrated  (Equation  34.5)  to  find 

normal  velocities,  un,  and  integrated  again  for  displacements,  un  (Equa-  j 

tlon  34.6) . | 

As  in  case  of  the  flat  plate,  boundary  conditions  and  edge  distri- 

' I 

butions  are  satisfied  by  Equations  27,  28,  35.1  and  35.2.  I 

The  details  of  the  algorithm  (marching  procedure)  for  calculating 

the  time  dependent  variables  of  Equations  66,  are  similar  to  the  ones 

for  the  flat  plate  (page  53). 

c.  The  Geometry 

The  principal  geometrical  dimensions  of  the  F-lll  windshield  were  j 

i 

either  found  in  original  drawings  or  estimated  as  closely  as  possible.  The  j 

actual  radii  for  bow  and  arch,  Rmin  and  \aK>  are  given  in  Figure  25 
(Notations  in  capital  letters,  such  as  RMIN  and  RMAX,  are  used  for 

these  and  all  the  other  variables  in  the  FORTRAN  program  of  Appendix  IV).  I 

The  central  angles  of  the  cone  segment  (Figure  25) , minimum  at  bow  and 
maximum  at  arch,  ^min  and  s^iax*  are: 
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0min  = 180°-2  arc  sin  ((26.1-10.9)/26.1)  * 108.9°  (67.1) 

0max  * 180°-2  arc  sin  ((30.5-21.8)/30.5)  = 146.8°  (67.2) 

The  geometry  of  the  window  is  further  described  by  the  thickness, 
t (in),  of  the  window  pane.  Lbea  is  the  length  of  the  head  beam  (in), 
and  Ls:j^  the  length  of  the  sill  (in)  (Figure  26) . 

Estimated  effect ive  widths  for  generation  of  the  finite  difference 
grid  for  the  edge  members  are: 

^hea=^’  ^sil=3<5,  ^bow=''’  ^arc=2  inch. 

d.  Stiffnesses 

Pane  stiffness  is  reflected  in  Young’s  modulus,  Epane (pound/ 
in2),  and  bending  stiffness  per  inch  of  width  in  the  longitudinal 
direction  for  head  beam,  sill,  bow  and  arch,  EIhea,  Elsii,  EIbow  and 
Elarc  (pound-in)  and  transverse  per  inch  of  length,  Elthdf  Eltsl> 

EItbw  and  EItac  (pound-in)  (Figure  27). 

Estimated  effective  longitudinal  and  transverse  bending  stiffnesses 
(Reference  45)  for  the  frame,  including  the  titanium  sheet  reinforcement 
strips  (.03"  x 2.2"  at  upper  and  .045"  x 2.2"  at  lower  side)  with  a pane 
thickness  of  .4"  are  EIbea=6,040,000;  EIsi^=66,000,000;  EIbow=20,000,000 
EIarc=2,760,000  (pound-in2)  along  the  frame.  The  bending  stiffness 
effect  of  edge  reinforcements  with  the  titanium  stripping  accounts  for 
160,000  pound-in2  along  the  frame  and  is  included  in  the  above  values. 
Table  III  presents  the  material  properties  for  polycarbonate  (Reference 
46).  Values  of  EIs_q  and  EIbow  are  estimates  rather  than  precise  data, 
because  of  the  unknown  supporting  stiffnesses  provided  by  the  structural 
opening  of  the  module. 
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ALUMINUM  2024 


Figure  27.  F-1U  Windshield/Mounting  System 
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The  EIt  values  In  the  transverse  direction  are  much  smaller  than  the 

El  values  in  longitudinal  direction  and  represent  partly  the  transverse 

bending  stiffnesses  of  the  flanges  to  which  the  reinforced  glazing  edges 

are  mechanically  fastened.  Stiffness  values  of  critical  cross  sections 

( I/in=(flange  thickness)  / 12  ) are  thought  too  low  to  be  representative 

of  mounting  stiffnesses.  Considering  one  half  the  thickness  of  the 

2 2 

standard  glazing  of  0.4  in  , I/in“flange  thickness* (0.2")  was  assumed. 
The  resulting  values  are:  EIthd  = 73,500;  Eltsl  = 84,000;  Eltbw  =84,000; 

Eltac  = 64,000  (pound-in). 

Bending  stiffness  in  both  the  circular  and  longitudinal  direction 
of  the  conical  windshield  are: 

EIX  = Epane  (t3/12)  (Ay)  (68) 

Ely  - Epane  (t3/12)  (As)  + AEIy  (69) 

AEIy  = (Curv)*Epane  • t (1-cos  (A^/2))2-r2  • (As)  (70) 

where  As  (=r*A0)  is  the  finite  difference  along  the  arch,  and  AEIy 
is  an  effect  of  curvature  contributing  to  the  gross  longitudinal 
stiffness.  Considering  the  finite  difference  element  as  an  arch,  the 
neutral  axis  is  located  parallel  to  the  base  of  the  arch  at  a distance 
less  than  the  height  of  the  arch.  Ignoring  the  first  term  in  Equation 
69,  the  fraction,  Curv,  was  graphically  determined  to  be  or  less. 

Irt  the  analysis,  a variation  of  Curv  between  0 and  0.2  has  been  used. 
Additionally,  Curv  is  assumed  to  vary  as  a function  of  arch  deflection 
to  represent  the  realistic  situation. 
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e.  Boundary  Conditions 

The  attachment  conditions  of  the  window  pane  to  the  frame  and  of  the 
frame  to  the  support  opening  of  the  module  structure  is  simulated  by 
elastic  spring  type  boundary  conditions  for  the  sill,  the  bow,  and  the 
arch  (excluding  the  head  beam).  The  sill  and  the  bow  are  considered  to 
be  rotationally  as  well  as  transversely  very  stiff  compared  to  the  arch. 

A 100:1  stiffness  ratio  was  assumed  with  reasonably  flexible  values  of 
1000(pound/in)  and  1000 (pound/rd)  for  the  arch  resulting  in  the  following 


values : 


Kunsl  = 100 »000 
^unbw  = 100j00y 
^unac  = 1.000 


K0Sil  - 1°°.000 
Kobow  = 100,000 

K0arc  * I*000 


In  this  preliminary  study,  these  parameters  were  varied,  and  the  effects 
on  amplitude  and  frequency  responses  of  stresses  and  deflections  were 


analyzed. 

While  Curv  accounts  for  the  effect  of  the  curvature  on  bending  stiffness 
in  the  window  pane,  a correction  factor  (Cor)  is  used  for  the  actual  in- 
creased bending  stiffnesses  at  the  windshield  corners. 

The  viscous  type  translational  damping  at  all  four  edges  is  negligible 
in  the  conventional  windshield  frame  of  the  F-lll  fleet;  therefore  0.1 
(pound-sec/in)  is  considered  a reasonable  value.  Variation  of  the  effect 
of  this  property,  as  well  as  in  rotation,  has  been  explored. 


^unsl  ” 0*1 


^osil  = 0*1 


Kunac  = 0*1 
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f.  Mass  Properties 

All  mass  properties  are  based  on  the  following  specific  densities, 

<> (pound-sec/in4),  for  titanium  (0.000414),  aluminum  (0.000258)  and  poly- 
carbonate (0.000112). 

Estimated  mass  properties  of  frame  per  running  inch  (pound-sec2/ 
in2),  based  on  the  dimensions  on  original  drawings  (Reference  45)  and 
specific  densities,  for  the  head  beam,  sill,  bow  and  arch,  (with  the 
mass  of  the  reinforced  window  edges,  ~0. 0013,  added  in  each  case)  are: 

Mhea=  °-00054«  Msil"  °*00071>  0.00043  and  Marc  « 0.00038.  The 

types  of  mass  elements, AM,  at  the  corners  are: 

AM*  Ax  • Mbow+Ay  • Msil  (71) 

Along  the  frame, 

AM  ■=  Ax  . Mfaow  (72) 

Mass  properties  of  windshield  (window  pane)  calculated  from 
windshield  dimensions  are: 

(4M)w»s-(iS)*,/(‘y)x,xt  <»> 

with  y = 0.000112  and  t = 0,9. 

g.  Computer  Program 

The  computer  program  developed  for  special  features  of  the  windshield 
in  the  F-lll  module  (Appendix  III-4)  incorporates  the  following  possibi- 
lities: (a)  a choice  of  loading  history,  (b)  any  location  of  impact, 

and  (c)  any  boundary  conditions,  except  for  the  fixed  tangential  support. 
Typical  cases  for  the  F-lll  windshield  may  be  found  in  Section  IV-2. 
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SECTION  IV 
APPLICATIONS 

This  section  contains  the  numerical  results  obtain.  for  square 
panel  test  WT-18  and  the  F-lll  windshield  test  FM-1  using  the  finite 
difference  analysis  method.  These  results  are  compared  with  available 
experimental  data. 

1.  ANALYSIS  OF  SPECIFIC  SQUARE  PANEL  TEST  (WT-18) 

One  phase  of  a contractual  effort  by  PPG  Industries  (Reference 
47)  was  devoted  to  preliminary  impact  testing  of  18  flat  26"  x 26"  panels 
of  14  types  of  pane  construction.  The  panels  were  bolted  into  a test  frame 
at  an  angle  of  23°.  Except  for  two  failures  caused  by  the  frame  during 
sorting-out  runs,  results  were  predictable. 

The  following  is  a summary  of  the  18  test  cases  (Reference  47) . 

Two  of  the  specimens,  WT-1  and  WT-2,  were  monolithic  polycarbonate  and 
both  failed  by  shearing  of  the  panes  along  the  top  edges  of  the  frame. 

The  others  were  laminated:  three  of  three,  one  of  four,  two  of  five, 
six  of  seven,  one  of  eight,  and  three  specimens  of  nine  layers  each.  The 
materials  used,  in  order  of  application  and  with  range  of  thickness 
indicated  were:  polycarbonate  (0.093"  - 0.766"),  PPG  proprietary  interlayers 

(0.09"  and  0.12"),  stretched  acrylic  (0.06"  - 0.125"),  Herculite  II  (0.11"), 
and  as-cast-acrylic  (0.06"  - 0.125").  Total  pane  thickness  was  0.75"  thru 
1.041".  Pane  edges  were  simply  stepped  down  to  thicknesses  from  0.359" 
thru  0.451",  In  seven  panels,  the  pane  edges  were  composed  of  monolithic 
polycarbonate  and  eleven  panels  of  two  layers  of  polycarbonate  with  one 
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PPG  interlayer.  Total  panel  weights  varied  from  20.25  thru  30.5  pounds 
(specific  weights  of  4.68  thru  7.48  pound/ft^).  Weights  of  birds  were 
3.9i  to  4.18  pounds.  The  actual  speed  varied  between  96.5%  and  114.7% 
of  the  intended  speed.  Range  of  final  velocities  was  755.1  thru  908.2 
ft/sec  (447  thru  538  kt).  The  kinetic  energies  at  impact  were  from 
35,398  pound-fc  in  shot  WT-2  to  53,496  pound-ft  in  shot  WT-5. 

Eleven  failures  occurred  and  were  typically  described  as:  panel 
completely  or  center  of  panel  blown  out;  bird  penetrated  from  top  edge 
shear;  bird  bounced  but  had  45°  yaw  and  top  edge  of  frame  suffered 
permanent  set;  several  cracks  in  polycarbonate;  and  panel  bulged.  On 
nine  of  these  panels  the  permanent  set  was  from  1/16"  to  1/2".  Of  the 
seven  tests  with  no  penetration  occurring,  shots  WT-9  at  44,077  pound-ft 
total  kinetic  energy  of  impact  and  WT-18  at  45,353  pound-ft  were  accept- 
able with  good  residual  visibility.  The  critiques  of  the  remaining  five 
tests  are:  structural  deficiency;  improper  damping  frame;  9/16"  permanent 

set  and  0 degree  residual  visibility;  crack  in  polycarbonate  middle  ply 
and  7/16"  permanent  set;  spall  of  stretched  acrylic. 

Figure  28  gives  one  of  the  original  data  forms  prepared  for  each 
individual  test.  The  corresponding  set  of  deflection  histories,  cal- 
culated frommotion  picture  frames,  at  the  center  of  the  panel  and  4 
inches  back  of  the  center  are  shown  in  Figure  7.  Test  shot  WT-18, 
presented  in- Figure  28,  was  one  of  the  moot  successful  tests  and  has 
been  simulated  for  this  reason. 
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Sample  Code  : -7a 

Construction  • .ps'wut-.BoV.  .t8fiV-.o*VL-  .isa’pc- 

*ROvl/L-  .KS*f> C.  1 

Thickness  •.  ,95s"  ArealOensity  :5sd  ib./pt1  Panel  Wt.*.^.75ls. 
Edge  Section  : _j__  1 ~ 


A 


.444-, 461 

T 


.953 


(FKANT) 


( REAR.) 


Results  on  3 - 30-73  @ 4f*P  and  ^installation  an6le 
Shot  no. ••  wr* is  Spec d l Isj, s* Bird wr. ; &00 1*> ■ 
Remarks  * j - o < goc  kf)  ft-lb) 

S>*6  6t>o/vic.e6  . No  Pouy  cAhS-Bow/vr**  cpacic/axs  back 
Pnf . ^ 2>rteAWK/<;  «-P ackWg  ow  rop  cdce-, 

Pe^Pe^MOcE-  »P  An.  "PANELS  s«or  ???  AtPftc+.ZoVo  OP 
OorfcoAP-D  Plof  8w>u»n  ©F£.  Soot  AOHkz/oa. t t£UE<*»HcAE . 
■^ar^lvoAt-  OiSiBiuir^  E^csluckst, 


(FROM  REFERENCE  47) 


Figure  28.  Original  Test  Data  of  Shot  WT-18 
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a.  FINITE  DIFFERENCE  SOLUTION 

Using  Che  material  properties  given  in  Table  VII,  the  stacking 
sequence  of  the  laminate  is  0,125"  acrylic,  0.12"  interlayer,  0,188"  poly- 
carbonate, 0.09"  interlayer,  0.188"  polycarbonate,  0.12"  interlayer,  0.125" 
polycarbonate,  totaling  ■~0.956".  The  effective  bending  stiffness  of  the 

panel  of  shot  WT-18  may  be  estimated  using  Equations  22  and  23,  El2am  per 

2 

inch  width  ^24,000  pound-in  at  a corresponding  modulus  of  the  laminated 
panel,  EM0D  * 0.33  x 10^  psi,  was  used  in  the  analysis. 

Applying  the  equations  for  plate  vibration  (Equations  3,  4,  5),  the 
theoretical  vibration  frequency , (*)th,  is  compared  with  the  actual  fre- 
quency, (Jac,  as  determined  from  Figure  7.  The  plate  naturally  vibrates 

after  the  impact  disturbance  has  diminished.  An  estimate  is  2 7T/ 

_1  2 

(0.015  sec)  = 420  sec  . With  D>vEl2am  * 24,000  pound-in  , £ * 0.00011  pound- 
sec^/in^,  and  s * 26  in,  " 790  sec“l  (88%  higher  than  the  actual 
vibration).  The  vibrational  behavior  of  the  heavy  frame,  not  taken  into 
account,  may  explain  the  actual  slow  down  from  0.008  to  0.015  sec  per  cycle. 

Proper  choice  of  grid  size,  As,  and  time  step,  At,  according  to  the 
condition,  (At/As)*  (EA>)^^  >1,  to  achieve  accuracy  of  prediction  by 
finite  differences,  was  based  upon  Reference  48.  Using  a time  step, 

At*  0.0001  sec  and  an  applicable~|E/j  * 50,000  in/sec.  As  has  to  be  5" 
or  smaller.  The  square  or  rectangular  flat  plate  program  (Appendix  IV-1) 
originally  using  a grid  size,  NG  = 13,  and  in  consequence  the  mass 
quotient,  NM  ■ 37,  to  calculate  equal  mass  distribution,  DMM  = MPANE/NM, 
as  well  as  the  load  quotient,  NQ  = 70,  to  determine  unit  load,  QTtfT/NQj 
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I 


TABLE  VII 

MATERIAL  PROPERTIES  CONSIDERED  IN  THE 
ANALYSIS  OF  PANEL  SHOT  WT-18 


64  °P 


transparent 
polycarbonate 
plastic  sheet 
MIL-P-83510 


stretched 

acrylic 

sheet 

MIL-P-5425 


polyvinyl 
butyral  base 
flexible 
interlayer 


tensile 

strength  (ksi) 


10.3 


10. 


3. 


tensile 
modulus 
x 1ob  (psi) 


0.34 


0.43 


0.24 


Poissons 

ratio 


specific 
mass  density 
(lb  sec^/in^) 


(0.3) 


0.000112 


0.35 


0.000111 


0.46 


0.000102 
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for  pyramidal  impact  loading,  was  modified  for  selective  grid  sizes  up 
to  NG  = 33.  Table  VIII  contains  the  values  of  NG,  NQ  and  NM  for  the 
computer  input.  Grid  size  chosen  for  analysis  of  shot  WT-18  is  NG  ■ 21. 
Table  IX  shows  results  for  an  unrelated  case  of  the  occurrences  of  maxi- 
mum/minimum pane  center  stresses  and  deflections  for  different  grid  sizes. 
Values  can  be  seen  to  converge  at  finer  grid  sizes.  For  grids,  NG2:25, 
the  program  was  found  to  be  unstable.  The  cause  of  this  instability  was 
not  investigated. 

The  test  frame  for  the  26"  square  panel  shots  was  made  of  1/2" 
cold  rolled  steel  with  32"  x 32"  outer  and  24"  x 24"  inner  dimensions 
(Figure  29  a).  Attachment  bolt  holes  were  located  1/2"  aft  of  the  inside 
edge.  1/4"  diameter  bolts  of  60  ksi  strength  at  1 1/4"  center  spacing 
were  used.  Four  (4)  clamps  were  used  to  hold  this  frame  with  the  mounted 
pane  on  a very  rigid  support.  Positioning  of  the  clamps  was  varied  ini- 
tially from  test  to  test  until  the  best  possible  clamping  positions  were 
established,  causing  no  plastic  deformation  of  the  frame  under  impact  to 
avoid  subsequent  refurbishment  of  ^the  frame  after  every  shot.  Clamp 
locations  were  fixed  from  shot  WT-6  on  (Figure  29  b) . Locations  of  clamps 
were  taken  from  the  impressions  of  the  clamps  on  the  used  hardware. 

At  a grid  size  of  NG  ■ 21,  the  clamping  idealization  is  shown  in 
Figure  29  b where  the  bolt-center  lines  (25"  x 25")  representing  the  lines 
of  support  are  equivalent  to  one-directional  support  of  the  frame  only. 
This  arbitrary  decision  for  simplicity  of  support  condition  neglects  the 
rotational  behavior  of  the  frame. 
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TABLE  VIII 

EQUAL  MASS  AND  PYRAMIDAL  LOAD  DISTRIBUTIONS  FOR  DIFFERENT 
PRACTICAL  GRID  SIZES 


Grid  Size 
NG 

Quotient  NM 

MPANE 
DMM  = m 

!_  - - 

Quotient  NQ 
Unit  Load  = 

NQ 

5 

5 

2 

9 

17 

2o 

13 

37 

7o 

17 

65 

168 

21 

1.0 1 

33o 

25 

145 

572 

29 

197 

91  o 

33 

257 

11 6o 
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TABLE  IX 

COMPARATIVE  RESULTS  AT  DIFFERENT. GRID  SIZES 


grid 

size 

max/mi n 
center 
stress+ 
(ksi) 

time  of 

occurence 

(sec) 

max/min 

center 

deflection 

(in) 

time  of 

occurence 

(sec) 

5 

-6o.  7 

o.  oo5 

33.3 

o.oo5 

61  ,o 

o.  ol  5 

-33.2 

o.  ol  5 

9 

-38.6 

0.005 

16.0 

o ,co4 

4o.  2 

o.ol  2 

-16.3 

o . ol  2 

13 

-31.4 

o.  oo5 

13.4 

o.oo4 

34.6 

o.  ol  1 

-13.5 

o.  ol  1 

r? 

-31.7 

o . oo5 

12.3 

o.oo4 

31.4 

0.010 

-12.1 

o.olo 

21 

-31.6 

o . oo3 

11.6 

o.oo4 

31.5 

o • ol  1 

-11.7 

o.olo 

(discontinued  because  of  instability  of  parameters  of 
computer  program  (blow-up)) 

(+print-out  recorded  every  10th  step) 

panel  site  25"  squared 
time  step  0,0001  sec 
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The  calculated  bending  stiffness  of  the  frame,  using  the  formula  EIF  * 

Esteel  x wi£jth  x thickness3/12 , is  approximately  1.3  x 10®  pound-in2. 

Total  mass  of  the  pane  without  frame  based  on  the  weight  given  in  Figure 

2 

28  is  0.0615  pound-sec  /in. 

-3 

Total  mass  density  of  the  frame  is  £ steel  x volume  * 0.72  x 10  x 
2 

224  ~0.16  pound-sec  /in.  The  overall  spring  stiffness  of  the  frame  in 
the  downward  mode  of  deflection  is  assumed  to  be  ky  * 100,000  pound/in 
and  a very  small  number  in  the  upward  mode.  Spring  stiffness  in  the 
upward  mode  due  to  clamping  at  the  eight  (8)  idealized  locations  (Figure 

29  b),  is  100,000  pound/in.  For  a bird  velocity  of  854.7  ft/sec  for  WT-18, 
vertical  impact  and  parallel  sliding  velocities  at  an  inclination  of 

23°  are  334  ft/sec  and  786  ft/sec,  respectively. 

Using  '.he  normal  velocity  and  the  mass  of  a 4 pound  bird  in  the  force 
equation  of  Reference  42,  the  calculated  force  is  16,120  pounds.  The  im- 
pact time  is  0.00133  sec  cr  about  13  time  steps  of  At«  0.0001  sec  each.  The 
effective  bird  dimension  is  ^5".  Due  to  central  Impact,  the  sliding  time 
is  0.00125  sec  or  longer.  This  yields  the  acting  time  of  the  impact  pres- 
sure mound  almost  as  long  as  it  takes  to  slide  across  the  pane  towards  the 
upper  edge  of  the  frame.  A two  step  pyramidal  load  distribution  was  used 
in  the  analysis  to  account  for  the  traveling  pressure  mound.  It  has  a 
pyramidal  central  distribution  lasting  for  time  steps,  NNNN«1  thru  7,  and 
this  distribution  shifts  midway  between  the  center  of  plate  and  the  upper 
frame  (shift  J«J+6)  for  NNNN*8  thru  12. 

The  following  are  the  input  data  used  to  calculate  the  deflections 
at  12  different  locations  simulating  shot  WT-18  (Figure  30) : 
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Grid  parameters  - NG=21,  NM=10].,  NQ=330 

Panel  dimensions  - X=Y=25 

Panel  thickness  ~ THICK=.956 

E-module  of  laminated  panel  - EM0D=.329  E+6 

Bending  stiffness  of  frame  - EIXF=EIYF=,129  E+6 

Mass  of  pane  element  - DMM*»0.0615/NM 

Mass  of  frame  element  - DMF*«0.1613/(2*NG-2) 

Total  impact  load  - QTOT=25,000 
Time  step  - DT=0.0001 

Spring  stiffnesses  chosen  representing  downward  fixed  support  and 
upward  free  support  at  the  eight  clamping  locations  were:  KW=100,000 

(fixed)  or0.1~10  (free).  Deflection  data  were  calculated  for  400  time 
steps  totaling  0,04  sec.  The  data  are  plotted  in  Figure  30.  The  "Y"  and 
marked  curves  in  Figure  30  correspond  closely  with  the  test  curves  in 
Figure  7. 

b . PRELIMINARY  OPTIMIZATION  OF  MOUNTING  FRAME 

A scalar  analysis  of  windshield/mounting  frames  had  been  conducted 
by  means  of  the  analog  computer  (or  MIMIC)  simulation  language  (References 
49  and  50)  as  discussed  in  Section  II-l.  Series  (A)  and  (B)  runs  for 


r.n 

Pl  =! 

f.i' 

i 

•,  Kp/Ky=<  1 } and  Cw/CmH 

14’ 

k 

UoJ  1 

(totaling  twenty-seven  (27)  runs  each)  were  conducted  to:  (1)  find  the 

range  of  parameters,  Mtf/Mp,  I'p/lfy  and  Cy/C^,  representing  quasi-rigidiy 
mounted  conditions;  (2)  find  a difference  between  calibration  assumptions 
(A)  and  (B) , and  (3)  optimize  mounting  for  lowest  pane  deflections  relative 
to  mount  deflection,  which  is  equivalent  to  lower  dynamic  stresses  in  the 
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pane.  Follow-on  runs  were  made  for  further  optimization,  an  assessment 
of  FQtj/My  = const,  and  vibration  of  the  coupled  system. 

Using  the  values  M;,j/Mp=  j^j  and  Kp/K^lO  at  Cp/C^*'  I resulted  in 
extreme  frame  deflections,  Xp,  of  only  10%  of  the  pane  deflections  relative 
to  ground  reference,  Xw.  Therefore,  deflections  Xy-Xp,  are  approximately 
the  same  as  in  the  unmounted  case  of  infinite  stiffness  and  mass  of  the 
frame  (i.e.  ratios  Kp/Kw  =00  and  M^/Mp  = 0). 

Figure  31.1  shows  the  eight-track  histories  of  three  cases  of 
calibration  (A),  while  Figure  31.2  shows  the  histories  to  calibration  (B) . 
Only  slight  individual  differences  in  the  two  series  of  traces  are  ob- 
served. Table  X identifies  these  analog  records  at  the  common  tracing  time 
factor  of  2mm  * 1 ms.  As  an  example,  deflection  variables,  Xp  (track  2), 
%-Xp  (track  51,  and  Xw  (track  7)  are  at  an  identical  scale  and  double  that 
of  Xw  (track  1). 

Figure  32  represents  deflections,  20Xp  at  2 V/L  and  10  (Xw-Xp)  at 
1 V/L,  on  tracks  (2)  and  (5),  respectively,  for  the  mounted  case  with 
calibration  (B) . Run  "1"  in  Figure  32.2  at  Kp/K^.l,  Mw/Mp=°10,  and 
Cp/CvplO  gives  the  lowest  window  deflections  (stresses)  at  these  parameters. 
Runs  "2",  "3",  "5"  and  "6"  reveal  that  the  frame  limit  deflection  criterion 
arbitrarily  set  at  approximately  2 inches  has  been  exceeded  in  early  time 
periods  with  the  traces  invalidated^  via  built-in  analog  circuitry  . 

These  analog  runs  were  followed  by  further  adjustments  of  the  three 
parameters  (Table  XI) . Use  of  the  digital  "MIMIC"  .code  (References  49, 

50)  was  supplemental  (Table  XII),  and  applied  to  cross-check  results. 


For  the  entire  series  of  runs  of  Table  XI,  the  values  10  (Xw  - XF)  at 
1 V/L  for  the  flexible  frame  remained  practically  constant.  10  (Xy  - Xp) 
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(1)  CALIBRATION  (A) 

Figure  31.  Eight-Track  Histories  of  Three  Analog  Runs 
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TABLE  X 

IDENTIFICATION  OF  EIGHT-TRACK  ANALOG  COMPUTER  RECORDS 


track  number 

scale  factored 
variable 

recording  parameter,  V/L* 
analog  volt  output  per  line 

1 

10  XW 

2 

2 

20  XF 

2 

X) 

0) 

4J  QJ 

S 11 

0 u 
B ^ 

3 

4 

.05  Xw 
.05  XF 

5 

2 

5 

io(xw-xF) 

1 

v .. 

6 

5 ‘ 10“5  Fq/Mj, 

2 

Xi 

•H 
60 
•H  1 
U ** 

u 

j « 
-1 
-4 

7 

8 

M , 

10  xw 

.05  Xw 

1 

5 

(See:  Figure  31) 
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TABLE  XI 

EIGHT-TRACK  PARAMETERS  FOR  SERIES  (B)  ANALOG  COMPUTER  RUNS 


Run 

KF/Kw 

mw/«f 

cF/cw 

II 

10  (xw-xF) 

at  1 V/L 

1 

0.1 

1 

100 

7 

{lines] 

13.4 

2 

0.1 

100 

8 

/ 

k 

3 

0.1 

100 

9 

4 

0.03 

100 

9 

5 

0.01 

100 

9 

6 

0.1 

50 

9 

7 

0.03 

50 

9 

8 

0.01 

50 

9 

9 

0.1 

10 

9 

0.05 

10 

9 

> 

i 

11 

0.01 

10 

9 

14.8 

Series  (B) : (Ky/My)^^=  467  sec“^ 

Cy/My  “ 76  sec"1 

Fq.%  = 1.2x10^  in/sec^  at  = 10“^  sec 


AFFDL-TR-74-155 


TABLE  XII 

MAXIMUM  DEFLECTIONS,  XWF,  FROM  MIMIC  X-Y  PLOTS  FOR  SERIES  (B)  PARAMETERS 


* 


» 


*»*  MlMjC  SOURCE  PROGRAM 


CONSTANTS  are  M«MW/mF/  K2*KF/KWj  C2*CF/CW 
C0N(M,K2*C2) 

CONSTANTS  AR£.K1«(Kw/MW)**2#  C1«CW/MW,  ZR»ZER0«Q 
- " C0N(Kl-iCliZR) 

par  < Fw/ t l ) 

*DT -♦0005' 

F ( T ) *FW*MW  T< T I A «ZR  T»  T I A «ZR  T>TI 

F0T«.Fm- 

F0T  FSW(T-TliFW*ZS/ZR> 

K3  K1#m 

K4  K3*(l#*i<H) 

*C3 Ci*M - - 

C4  C3* ( 1 « +C2 ) 

the  mounted  frame 

D2XW  F0T-K1#(XW-XF)-C1*(D1XW»DXXF) 

D2XF  K3*XW«iW*XF4-C3*0lXW«C4*0lXF 

01XF  lNTU)2XF#0t) 

-01-XVf lNT-(-^Xw,0»'>  ■ -•  - 

XW  INT ( OlXto#Q# ) 

XF  INT(D1XF<0.) 

XWF  XW-XF 

FIN( T*  #04 } 

OUT 

eUT-(-fr-i  XWF-i  -rXf-> - 

PL0(T,XWF,XF/F0T#FW#TI) 

* - END 


Run 

VKw 

Vi 

VCW 

^WF=^W~Xp 

fin.] 

1 

.1 

10 

10 

1.5 

2 

9 

1.5 

3 

8.5 

1.45 

4 

8 

1.4 

5 

7 

1.38 

6 

5 

1.2 

7 

1 

.6 

8 

.1 

.7 

Series  (B) : FW  = FQ/MM  * 1,200,000  in/sec^ 
K1  * Ky/My  = 2.18089x10’’  sec“: 
Cl  = CW/MW  = 76  se<H 


TI  = = .001  sec 

ZR  = 0. 

M = My/Mp 

K2  = Kf/Kw 
C2  = CF/C„ 
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for  the  rigid  frame  is  about  1.53  times  higher  than  the  flexible  values 
10  (Xw  - XF).  X - Y plots  of  runs  "1",  "5",  and  "7"  (Table  XII,  Figure  33) 
reveal  the  penalties  of  high  frame  deflections,  XF,  for  low  window  pane 
deflections,  XWF=Xy-XF;  for  the  latter  two  runs,  MIMIC  run  "1"  of  Figure 
33  is  identical  with  analog  run  "1"  of  Figure  32.2. 

The  results  of  this  preliminary  optimization  study  of  the  frame  pro- 
perties, Kf/Km,  M^/Mj,  and  Cp/Cw,  are  useful  in  conducting  a detailed  finite 
difference  optimization  of  the  window  mounting  system  that  best  suits 
impact  case  WT-18. 

A study  of  varying  F0  and  tj  (keeping  Fotj/My  = const)  was  conducted 
by  analog  and  duplicated  by  MIMIC  runs  (Table  XIII).  Resul.s  reveal  a 
decrease  in  maximum  excited  amplitude  response,  (Xw-Xp)maxI  bY  lowering 
Fo/My  and  elongating  tj.  This  simulates  heavier  birds  at  lower  impact 
velocities.  Table  XIII  shows  that  the  effort  of  mounting  the  window  for 
shock  absorption  In  case  of  the  heavier  birds  becomes  inefficient. 

9 

It  is  important  to  note  that  1/ 2 • • (X^-Xp)  mav  is  the  maximum 
potential  strain  energy  per  unit  window  mass  accumulated  during  the  first 
downward  mode  of  deflection  after  impact. 

Vibration  of  the  coupled  system,  with  the  traces  Xf  and  X^-Xp  nearly 
in  phase,  may  be  seen  in  Figure  32.1.  At  2 lines  per  ms,  runs  "2"  and; 

"3"  (KF/Ky=l,  M^/Mf=10)  indicate  38  1/2  lines  per  cycle  or  about  52  Hz. 

Run  "6"  (KF/Kw“My/MF  =*1)  yields  46  Hz'v289  rad/sec.  Using  the  system 
constant,  (Kw/M^)^/^  ■ (218089)^^  ■ 467  sec-^,  the  upper  A2_lioe  in 
Figure  4 (eigenvalues),  and  the  (h)  - line  (of 64)3)  in  Figure  5,  frequen- 
cies of  the  quadratic  eigenvalue  problem,  satisfy  the  conditions  of  run 
"6"  of  Figure  32.1. 
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Figure  33.  MIMIC  X-Y  Plots  of  Corresponding  Window 
Deflections,  X,—  , and  Frame  Deflections 
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TABLE  XIII 

SERIES  (B)  ANALOG  AND  MIMIC  RUNS  FOR  SQUARE  IMPULSE, 
(VMW>  •’ri=CONSTANT 


Constants: 


(%/My) 1/2  = 467  sec--*' 

F0ti  ■ 1.2*10^  in/sec 

CW/MW  = 76  sec  ^ 

Kp/Ky  = 0.1 

s 

» 

O 

Cp/ Cy  = 9 

Analog  Results  (duplicated  by  MIMIC) : 


Run 

mounted  case 

rigid  case 

10(Xw~Xp)max 

10xWmax 

at  1 V/L 

at  1 V/L 

[lines] 

[lines] 

1 

0.1 

15 

23 

2 

1 

14.8 

22.8 

3 

5 

12 

18 

4 

10 

7 

9.8 

(10  lines  ■ 1 inch  deflection) 
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2.  THEORY/TEST  CORRELATION  OF  WINDSHIELD  (F-lll) 

The  finite  difference  analysis  method  (see  Appendix  IV-4  for  computer 

program)  was  used  to  simulate  the  experimental  strain  gage  readings  of  the 

windscreen/canopy  of  a bird  strike  exposed  F-lll  module,  with  ftages  mounted 

at  designated  locations  in  the  vicinity  of  the  T-intersection  of  arch  and 

head  beam  (Reference  30) . The  glazing  material  of  both  F-lll  windshields  was 

monolithic  polycarbonate  (0.9"  thick),  impacted  by  a 4 pound  bird.  The 

records  of  the  following  shots:  (1)  on  the  center  of  one  of  the  windshields, 

FM-1,  (2)  one  at  the  corner  of  the  windshield,  FM-2,  and  (3)  one  at  the  center 

of  the  head  beam,  FM-7,  were  available  for  analysis  (Table  XIV).  Shot  FM-1 

is  compared  herein  with  the  theoretical  results. 

Theoretical  impact  loads  for  shots  FM-1,  FM-2,  and  IM-7  were  estimated 

using  the  empirical  power  laws  of  References  33  and  34.  They  are  presented 

in  Table  XIV.  The  calculations  for  shot  FM-1  are:  Impact  velocity, 

VEM-  j”932  ft/sec,  and  angle  of  inclination, ^22° , give  normal  velocity, 

Vn*349.5  ft/sec,  and  parallel  velocity,  Vt“865  ft/sec.  The  bracket  of 

probable  impact  loading  [ (using  the  less  severe  MacAuley  equation) 

2/3  2 

F»  0.05835  x W (pound)  x V (ft/sec),  and  the  sinusoidal  average  of  the 
original  Mitchell  equation,  2 Fj/ir"  0.1156  x W2^3  (pound)  x V2  (ft/sec)]  is 
17,970  thru  47,900  pounds.  Duration  of  impact  is  given  by  bird  dimension/ 

Vn  ■ 0.00122  sec  (MacAuley)  and  3/8  bird  dimension/Vn  ■ 0.00045  sec  (Mitchell). 

The  common  impulse,  force  x time,  is  21.85  pound-sec.  Time  of 
travel  is  distance/Vt«. 00218  sec.  The  bird  disintegrated  after  having 
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Li 


covered  .0012/. 00218-55%  (MacAuley)  and  .00045/. 00218-20. 6%  (Mitchell)  of 
the  distance  traveled  over  the  windshield. 

High  speed  photography  of  test  FM-1  indicated  that  the  impact  load  1 

corresponds  closer  to  the  MacAuley  load  (17,970).  The  load  initially  was 
set  at  an  arbitrary  127,200  pounds  times  a factor  FRC  (see  notation  FRC  in  » 

Appendix  IV-4).  FRC  « .15  was  chosen,  yielding  19,200  pound  of  load  at  ; 

i 

initial  Impact,  to  match  the  experimental  with  the  theoretical  curves  j 

i 

(Figures  35  and  36).  j 

The  idealized  centers  of  impact  in  the  nodal  grid  and  the  subsequent 
spread  of  the  traveling  debris  over  the  windshield  area  were  determined 
approximately  using  pertinent  data  from  Table  XIV.  They  are  given  in 
Table  XV.  Table  XVI  gives  the  input  time  histories  of  the  normal  and  j 

tangential  load  components,  Qn,  and  Qt,  respectively  for  the  "F-lll  wind-  1 

* 

shield"  subroutine,  i.e.,  J 


«Ux,y,i  * a*,y,t 

= 2 Qx,y<i  0K+t,y 


Figure  34  taken  from  Reference  35  shows  the  strain  gage  location 
for  shot  FM-1  and  the  corresponding  idealized  nodal  point  coordinates  at 
which  time  histories  of  bending  moments  were  determined.  These  bending 
moments  were  transformed  into  strains  for  the  purpose  of  comparison. 

The  FM-1  case  of  the  F-lll  windshield  is  analyzed  at  fixed  grid  size 
(17  x 17)  and  with  the  window  pane  and  frame  parameters  fixed  at  all  trans- 
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Figure  34.  Strain  Gage  Locations  for  Case  EM-1 
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lent  time  steps  of  comoutation.  The  input  data  for  the  computer  program 
had  to  be  assumed  at  various  levels  of  confidence,  by  trial  and  error. 

The  mass  parameters  of  the  bow  and  the  sill  (MB0W  and  MSIL) , the  constant 
representing  viscosity  of  the  pane  (VISC) , and  most  of  the  stiffness  and 
damping  properties  of  the  boundaries  at  bow,  sills  and  arch  locations  had 
to  be  experimented  with.  Other  local  adjustments  necessary  to  match  the 
theoretical  results  with  the  corresponding  oscillograph  traces  from 
Reference  35,  transcribed  in  Figure  35,  were  directly  programmed  in  the 
computer  code  (Appendix  IV-4:  HY(9,7J‘‘EIY(9,7)*2.  thru  EIX(5,17)= 

EIX(5,17)/2.) . 


Figure  35  a shows  three  strains  in  the  y-direction  (direction  of 
flight):  £y  9,9,  £y  9>17  and  £y  ? >1?.  The  last  trace,  £y  ? 1?,  should 
be  considered  partially  invalid  since  the  strain  gage  failed  during  the 
test.  These  three  traces  correspond  to  the  series  of  curves  shown  in  the 
upper  portions  of  Figures  36.  In  Figure  35  b,  four  strain  gage  readings  perpen- 
dicular to  the  longitudinal  axis  of  the  airplane  are  recorded.  These  strains, 
£x  9,13’  ^x  13,17’^x  9,17  ant*  £x  5 17 » correspond  to  the  four  curves  shown 
in  the  lower  portion  of  Figure  36. 

As  can  be  seen  from  these  figures,  discrepancies  between  the  experimental 
(Figure  35)  and  theoretical  results  (Figure  36)  were  experienced.  The  inac- 
curacy of  the  model  may  be  one  reason  for  this. 

In  case  of  corner  shot  FM-2  at  very  low  yield  of  the  stiffer  corner, 
acting  load  and  time  come  close  to  the  critical  Mitchell  condition.  In 
case  of  the  strain  gage  data  of  the  center  beam  test,  FM-7,  damage  from  pre- 
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vious  tests  and  necessary  strapping  may  have  obscured  performance  of  the 
module. 

The  computer  input  data  for  case  FM-1  are  shown  in  Table  XVII. 

These  data  are  in  computer  notation  and  the  reader  is  referred  to  Appendix 
1V-4  for  their  identification. 

A total  of  twenty-two  stiffness  and  dampening  parameters  are  considered 
in  the  computational  schema  of  FM-1.  A parametrical  study  with  just  three 
value  changes  (low,  medium,  high)  of  each  of  these  parameters  would 
require  about  3.135  x 1010  individual  runs.  It  was  determined,  how- 
ever, that  local  vibration  characteristics  can  be  influenced  easily  by  local 
as  well  as  overall  parametrical  changes.  The  parameters  of  the  individual 
frame  components,  for  example,  the  transverse  and  rotational  dampening  pro- 
perties of  the -arch,  KUNDAC  and  K0DARC,  would  have  local  influence  upon  the 
arch  area.  The  parameters,  EPANE  and  VISC,  on  the  contrary,  represent  pane 
material  characteristics  and  have  overall  influencing  effects. 

The  strain  gage  data  of  case  FM-1  was  the  only  one  readily  available 
for  use  in  the  comparison  of  test  data  with  theory.  For  this  reason,  only 
a small  series  of  specific  parameters  was  chosen  to  investigate  their  influ- 
ence upon  the  strains  at  locations  at  which  the  strain  gage  data  were 
recorded.  Representative  dampening  parameters  were  selected  to  in- 
vestigate the  shock  absorption  provided.  Parameters  for  the  arch,  repre- 
senting a part  of  the  frame,  were  chosen  because  it  has  much  smaller 
stiffness  than  the  head  beam  and  therefore  displaces  more  easily  than  the 
other  frame  components,  and  dampening  is  effective  only  due  to  displacement. 


i 
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TABLE  XVII 

INPUT  DATA  CHART  FOR  CASE  FM-1  OF  BIRD 
IMPACTED  F-lll  WINDSHIELD 


AFFDL-TR-74-155 


The  parameters,  in  the  nomenclature  of  the  computer  program  (Appen- 
dix IV-4) , are  EPANE,  VISC,  KUNDAC  and  K0DARC.  They  are  assumed  equally 
distributed.  EPANE  = 315,000  pound/inch^  is  the  calculated  modulus  of  the 
pane,  and  the  second  chosen  value,  200,000  pound/inch^,  represents  a 
hypothetical  low-modulus  pane  material.  All  assumed  values  of  the  remaining 
three  parameters  are  simply  "low-medium-high";  VISC  - .1  x 10“-^,  5,  10; 
KUNDAC  = .1,  5,  10;  K0DARC  = .1  x 10“ 10  , 5,  10. 

A four  do-loop  set-up,  a computer  program  modification,  was  used 
to  generate  all  combinations  of  EPANE,  VISC,  KUNDAC  and  K0DARC.  The  diffe- 
rent graphical  plots  of  the  strains  in  Figure  36  contain  only  five  of  these 
combinations.  The  pane  stiffness  for  all  six  plots  is  315,000  pound/inch^. 
The  remaining  three  parameters  in  order  are:  (,1*10“10,  .1,  ,ri0“10), 
(.1-10“10,  .1,  10),  (.1-10“10,  10,  .1>10“10),  (10,  .1,  ,1*10"10).,  (10,  5,  5), 
and  (100,  .1,  .l^lO"^0) , equivalent  to  parametrical  runs  "1",  "3",  "7", 

"19",  "23"  and  a "hypothetical"  one.  The  time  span  of  T=.05  sec  was 
covered  by  time  intervals  of  DT  = .0001  sec  with  the  final  time  step  NEND= 
500.  The  strains  at  every  fifth  step  were  plotted  (NP=5) . Similarly, 
tables  of  strains  as  well  as  deflections  of  all  54  runs  were  printed.  The 
total  of  nine  strains  contain  the  seven  at  the  recorded  edge  locations,  and 
£ co  and  r , „ at  the  initial  impact  location.  The  latter  two  strains 
were  selected  assuming  that  they  represent  the  .Jiaximum  stretch  occurring 
in  the  pane.  The  print-outs  of  deflection  histories  consist  of 
three  groups  of  three  columns  each:  UN(5,9),  UN(5,13)  and  UN (5, 17)  tracing 

the  route  of  impact  on  the  impacted  side,  UN(9,9),  UN(9,13)  and  UN (9 ,17) 
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along  the  head  beam,  and  UN(13,9),  1315(13,13),  Uh'(13,17)  at  symmetrical 
locations  of  the  opposite  windshield.  A matrix  print-out  of  all  variables 
(QN,  QT,  0X,  0Y,  MX,  MY,  V,  FN,  FT,  UN,  UND,  UNDO,  UT)  was  generated  for 
run  "1"  at  every  50th  time  step  (NP=50)  to  check  the  numerical  ranges  of 
all  time-dependent  variables.  The  original  17x17  matrices  were  printed  in 
a reduced  (clear)  format  containing  no  mid -way  zero  values. 

These  are  general  comments  about  the  original  strain  data  as  compared 
with  the  six  series  of  simulated  strains.  The  transcripts  of  oscillograph 
traces  seem  to  have  much  natural  dampening  when  matched  by  the  parametrical 
means  provided  in  the  model  of  FM-1.  The  strains  in  x-  and  y-  direction 
oscillate  quite  smoothly  with  lew  peaks,  and  smooth  vibration  at  the  edges 
and  corners  is  required  not  to  crack  and  destroy  the  attachments.  The 
simulation  with  VISC=100,  KUNDAC=.l  and  K0DARC".l  x 10--^  (~0)  looks  favor- 
able. A high  value  of  VISC  (Figure  36  f)  has  an  effect  on  both  the  traces, 
longitudinal  and  lateral.  However,  it  is  felt  that  high  overall  dampening 
represented  even  by  VISC-10  does  not  exist  in  the  pane  material. 

Another  initial  observation  is  that  the  y-group  and  x-group  of  strains 
go  independently  thru  similar  changes  (or  one  group  is  not  affected  at  all) 
due  to  the  individual  change  of  one  parameter.  The  reason  that  the  four 
lateral  strain  gages  oscillate  at  higher  modes  of  vibration  than  the  three 
longitudinal  ones,  an  observation  from  all  the  charts  of  Figure  36,  is  a 
greater  stiffness  in  flight  direction  than  in  lateral  direction.  Two  of 
the  longitudinal  strain  gages  bend  up,  while  one  bends  down,  and  three  lateral 
gages  bend  up,  while  one  bends  down.  This  is  due  to  the  gage  location. 
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A comparison  of  run  "1"  with  "19"  (Figure  36  a and  d)  shows  Chat  VISC=10 
has  little  effect  on  changing  strains  without  dampening  provided  by  KUNDAC  and 
K0DARC.  The  three  y-strains  in  "19"  approach  the  original  recordings  of  "1" 
(Figure  35) . High  translational  rather  than  high  torsional  dampening  (see  case 
"7"  vs.  case  "3")  is  more  favorable  for  matching  the  strain  gage  recordings 
in  the  x-direction.  The  torsional  dampening  in  "3"  endangers  the  x-direc7 
cion  with  high  poi.k  cyclic  oscillation.  The  peaks  are  removed  in  run  "7" 
with  KUNDAC=10,  but  the  frequency  is  still  too  high  to  match  the  experiment. 
Comparison  of  run  "7"  with  run  "23"  shows  that  the  x-strains  in  both  the 
cases  are  almost  identical,  Indicating  that  similar  effects  can  be  achieved 
via  different  means.  Using  run  "23"  as  a new  start  for  curve  fitting, 

VISC=20  Instead  of  10,  KUNDAC=10  instead  of  5,  and  K0DARC=2  instead  of  5 
can  be  used  for  comparison. 

The  following  comments  result  from  an  evaluation  of  the  complete  matrix 
print-out  for  run  "1".  All  the  matrices  reflect  non-symmetry  due  to  the 
unsymmetrical  loading  condition.  The  matrices,  printed  every  50th  time  step 
for  500  steps,  showed  stability  of  solution  during  this  time  span.  The 
results  are  of  the  same  order  of  magnitude  for  all  54  different  runs, 

The  slopes,  0X  and  0Y,  show  maximum  values  of  about  -.2,  initially,  at  the 
impact  location.  Later,  at  the  450th  time  step,  the  £.1  limit  for  0X  is 
exceeded  at  the  opposite  windshield  area,  This  is  due  to  the  fact  that 
run  "1"  contains  practically  no  dampening  of  oscillations.  On  the  second 
sheet  of  the  print-out  format,  the  matrices  of  moments,  MX  and  MY,  were 
calculated  with  values  up  to  -.2  x 10^  pound-in.  These  were  also  computed 
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at  the  opposite  side  of  impact  from  .03  sec  on.  The  highest  values  typically 
occur  at  the  sills  for  MX  and  at  the  front  bow  and  arch  for  the  MY  moments. 
Vertical  shear  forces,  FN  and  V,  generate  up  to  -.2  x 10^  pound,  and  tangen- 
tial shear,  FT,  up  to  ^.5  x 10^  pound.  The  third  type  of  print-out  sheet 
(UN,  UND,  UNDD,  UT)  initially  shows  the  high  values  on  the  side  of  impact. 
Later,  these  are  dispersed.  High  deflections,  UN,  in  the  pane,  head  beam 
and  arch  area  are  -(1  to  2)  inches;  velocities,  UND,  range  le9s  than  -10,000 
inch/sec;  accelerations,  UNDD,  less  than  ^10^  inch/sec^.  Tangential  dis- 
placements were  calculated  in  the  il  inch  range. 

The  strain  and  deflection  histories  of  the  27  runs  each  of  EPANE=315,000 

pound/inch^  (high)  and  200,000  (low)  have  been  compared.  Both  series  were  at 

o 

250  time  steps  with  every  5th  one  (50)  recorded.  EPANE=200,000  pound/inch 
caused  overall  lower  frequencies  of  vibration.  The  first  seven  strains  in  the 
nine  column  strain  tables  pertaining  to  the  metal  edge  frame  members  of  the 
windshields  indicate  values  with  all  strains  up  to  -.002  in/in.  The  two- 
directional  strains  in  the  pane  at  the  point  of  impact,  the  last  two  columns, 
are  ten  times  those  in  the  frame,  in  the  range  of  ^.02  in/in. 

It  Is  interesting  to  note  the  number  of  sign  changes  of  strains,  SX(5,9) 
and  SY(5,9),  in  the  tables  as  measures  of  frequency.  SX(5,9)  in  the  high 
stiffness  case  is  both  negative  and  positive  at  three  different  times,  SY(5,9) 
is  negative  seven  times  and  positive  six  times  at  higher  frequency.  During  the 
first  .025  seconds,  strains  SX(5,9)  and  SY(5,9)  are  a bit  slower  in  the  low 
stiffness  case;  negative  three  times  and  positive  two  times,  negative  six 
times  and  positive  five  times,  respectively.  The  four  x-strain  gages  oscillate 
at  about  the  same  frequency  as  SX(5,9);  the  three  y-strain  gages  at  about  1/3 
of  the  frequency  of  SX(5,9). 
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The  material  response  at  impact  grid  location  could  be  characterized 
as  cyclic  strain  both  of  high  frequency  and  magnitude.  Peak  values  for 
case  "1"  are: 

SX(5,9)  = -.0269,  +.0118,  -.0189,  +.00804,  -.00808,  +.00686  and  | 

SY (5,9)  = -.0275,  +.00448,  -.00563,  +.00553,  -.00757,  +.00081,  -.014 
+.0128,  -.0045,  +.00593,  -.00133,  +.00692,  -.00281  in/in. 

The  first  peak  values  of  strains,  SX(5,9)  and  SY(5,9),  are  practically 
unchanged  throughout  all  high  stiffness  cases.  The  peak  strain  reversals  j 

that  follow  the  initial  ones  are  very  irregular  from  case  to  case.  Higher  ! 

cyclic  stretch  may  be  found  in  the  x-direction  rather  than  in  the  y-direction  1 

of  windshield  location  (5,9).  j 

I 

Typical  first  strain  reversals  from  minimum  to  maximum  values  for  case  J 

"24"  (high  stiffness  and  combination  10,  5,  10)  are:  ASX(5,9)  = .03579  in/in 
and  ASY(5,9)  = .02895  in/in.  Second  strain  reversals  (max/min)  areASX(5,9) 

= .02479  in/ln  andASY(5,9)  = .00611  in/in.  The  same  strain  reversals  for 

the  low  stiffness  case  equivalent  to  "24"  are;  .0448  and  ,029  in./in  for 

SX(5,9);  .03538  and  .00744  in/in  for  SY(5,9),  This  shows  lower  cyclic  stretch 

for  the  stiffer  pane.  Here,  stretch  ratios  of  the  two  ca.ses  do  not  compare  j 

linearly  with  EPANE  stiffness  ratios, 

Of  further  interest  are  the  cases  of  cyclic  stretch  at  impact  location 
which  do  not  repeat  into  the  •>'.01  range.  This  is  the  case  for  SY(5,9)  in  runs  j 

"12",  "19",  "20"  and  "21“  (high  stiffness).  All  SX(5,9)  exceed  £.01  in/in  at  j 

least  twice  after  the  first  exceedance.  A similar  investigation  could  be  , j 

conducted  for  the  low  stiffness  cases.  > 

t 

i 
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The  high  and  low  modulus  def J ection  histories  (27  cases  each)  were 
investigated.  Every  deflection  table  contains  three  groups  of  three  columns 
each,  as  mentioned  previously.  The  first  group  of  deflections  on  the  impact 
side  are  characterized  by  higher  values  than  the  remaining  deflections}  up 
to  two  inches  in  the  downward  mode  (positive  values).  The  second  group  (head 
beam)  possesses  maxima  up  to  1.3  inch  downward  with  1 inch  at  the  opposite 
side  (third  group) . j 

The  histories  of  the  local  deflections  throughout,  the  series  are  gene-  ' 

rally  repetitive.  These  are  the  first  deflection  peaks  and  times  of  occurrence 

(for  runs  "1",  "19",  and  "23")  that  have  been  found  to  differ  to  the  largest  > 

| 

extend  in  the  first  high  deflection  group:  j 

UN (5, 9)  = 1.9K.002)  / 1.85 (.002)  / 1.85(.002);  . \ 

f 

UN (5, 13)  » 1.9 (.0045)  / 1.8(.0045)  / 1.8(.0045); 

UN(5 ,17)  = -.399(.0055)  / -.368(.0055)  / -.145(.005). 

The  second  peak  deflections  are  somewhat  similar.  The  first  peak  deflections 
for  the  corresponding  cases  of  lower  pane  stiffness  are: 

UN(5,9)  = 2. 45 (.0025)  / 2,36(.0025)  / 2.36(.0025) ; 

UN(5 ,13)  = 2 . 36 ( . 006)  / 2.21(.0055)  / 2.2(.0055); 

UN (5 ,17)  = -.396 (.0065)  / -. 36(.0065)  / -.124(.006). 

Here,  EPANE  = 200,000  pound/in^  can  be  seen  to  cause  higher  deflections 
than  315,000  pound/in^.  Typical  deflections  at  impact  location  are  1.91" 
and  2.45"  for  the  high  and  low  stiffness  cases,  respectively,  yielding  a 
ratio  of  1:1.28  while  the  pane  stiffness  ratio  is  1.58:1. 

I 

i 

It  may  be  concluded  that  the  edge  member  strains  in  the  impact  environ-  ! 

ment  of  PM-1  undergo  greater  changes  within  the  series  of  parametrical  runs 
(see  graphs  of  Figure  36)  than  the  corresponding  deflections. 
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SECTION  V 

RESULTS,  CONCLUSIONS  AND  RECOMMENDATIONS 

The  feasibility  of  high-visibility  bird-proof  windshield  design 
with  existing  glazing  materials  has  been  demonstrated.  Novel  window 
frame  concepts  with  sufficient  shock  absorption  capacity  to  accommodate 
part  of  the  impact  energy,  thereby  isolating  the  impacted  window  from 
damage,  were  examined. 

The  report  presents  the  finite  difference  method  of  designing 
bird-proof  windshields  efficiently.  It  shows  the  parametrical  design 
characteristics  when  using  the  developed  computer  routines  and  verifies 
the  analysis  method  by  simulating  actual  case  data.  The  finite  difference 
program  modeling  monolithic  linear  material  under  pure  bending  could  be 
replaced  by  models  of  other  material  behavior  including  laminated  ar- 
rangements. Modification  of  the  computer  routines,  however,  is  not 
justified  unless  actual  test  data  (combined  deflection  histories, 
strain  gage  readings,  and  details  describing  the  material  failure  modes) 
are  available.  A more  detailed  bird  impact  model  could  be  developed 
based  on  the  results  of  realistic  impact  pressure  distributions. 

A computer  program  was  developed  using  finite  difference  equations 
fo'r  the  square  window  pane  with  attachment  frame  and  shock  absorbing  mount- 
ing system  under  symmetric  loading  condition.  Requirements  on  mounting 
conditions  were  determined  by  screening  the  appropriate  parameters. 

Flat  plate  analysis  under  realistic  flight  conditions  of  1.2  Mach, 
an  impact  of  a 4 pound  bird,  and  technically  achievable  shock-absorption 
indicated  that  the  dynamic  peak  stresses  in  the  pane  could  be  reduced 
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by  30%  or  more.  This  can  be  done  by  using  a highly  viscous  damping  sys- 
tem in  connection  with  soft  springs. 

Damping  represented  by  coefficients  of  10  pound-sec/in  at  the  boun- 
dary nodal  points  and  spring  constants  of  25  pound/in  acting  at  the  same 
positions  were  found  to  be  technically  acceptable  for  the  impact  regime  of 
up  to  100,000  pound  of  0.5  milli-sec  duration.  Spring-back  deflections 
of  the  window  frame  went  up  to  one  inch  (a  technically  acceptable  compro- 
mise) . Up  to  90%  of  the  energy  was  absorbed  immediately  after  impact. 

Using  the  principle  of  conservation-of-energy  and  recording  the  time  history 
of  energies,  it  was  found  that  the  lowest  dynamic  stresses  in  the  pane 
occurred  for  the  mounting  case  for  which  highest  amount  of  energy  was 
absorbed. 

The  computer  program,  written  to  perform  stress  analysis  of  flat  wind- 
screen due  to  bird  impact,  was  modified  for  conical  shell  analysis  of 
fighter,  bomber  or  attack  aircraft  windscreen  canopies.  A change  from  four 
to  six  first-order  dynamic  difference  equations  was  necessary  to  account  for 
the  effect  of  curvature. 

The  input/output  formats  of  the  computer  routines  were  modified  to 
yield  more  fruitful  Information  for  design,  and  to  avoid  unstructured  runs 
in  the  hit  and  miss  experimentation  for  optimum  systems  performance.  Energy 
absorbed  and  retained  and  extrema  of  stresses  as  functions  of  frequency 
response  were  accounted  for  by  employing  available  mechanical  properties 
(such  as  those  for  transparent  polycarbonate  plastic  sheet) . 

Photomechanically  recorded  deflection  data  of  the  PPG  laminated  poly- 
carbonate 26"  x 26”  test  pcnel  (shot  WT-18)  were  successfully  simulated. 
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The  mass,  stiffness  and  dampening  parameters,  Mf/My,  K£/Kw  and  C^/Cy,  for 
optimum  mounting  conditions  were  determined  for  this  particular  case  using 
analog  and  . rMIC  analysis  techniques.  This  scalar  approach,  as  described 
in  the  report,  is  an  example  of  preliminary  optimization  of  the  best  suit- 
able window  mount  for  an  existing  windshield. 

Test  data  of  the  experimental  strain  gage  readings  of  the  bird  strike 
on  windscreen/canopy  of  an  F-lll  crew  module  (tests  FM-1,  FM-2  and  FM-7) 
generated  in  a bird  shot  program  at  Arnold  Engineering  Development  Center 
(AEDC) , were  provided  by  the  Improved  Windshield  Protection  ADP  office. 

The  conical  Linite  difference  computer  program  was  used  to  analyze  some  of 
the  data.  Results  have  been  matched  successfully  with  strain  gage  readings 
for  test  FM-1.  It  was  decided  not  to  analyze  data  of  test  FM-2.  In  case 
of  center  beam  hit  FM-7,  accumulated  damage  from  previous  tests  and 
necessary  strapping  obscured  the  performance  of  the  module. 

The  scope  or  recommended  experimental  work  is  to  develop  and  physically 
demonstrate  a shock  absorbing  window  screen  to  ensure  availability  of  the 
technology  for  future  systems.-  The  work  should  include  three  steps, 
namely  to  devise  innovative  concepts,  to  identify  absorbed  impact  energies 
and  to  reduce  maximum  stresses  of  window  pane  as  functions  of  the  design 
parameters.  Substantial  effort  should  be  devoted  towards  identification 
of  those  problem  areas  requiring  additional  research. 


142 


AFFDL-TR-74-155 


APPENDIX  I 

ASSESSMENT  OF  WINDSHIELD  STRUCTURAL  MATERIALS 

The  purpose  of  material  assessment  is  to  stimulate  bird  proof  design  by 
appropriate  stacking  of  highly  shock  proof  glazing  materials  and  viscous, 
plastic  interlayers  at  maximum  allowable  pane  thicknesses  for  optical 
qualification.  This  consideration  plus  feasible  end  attachment  to  the 
mounting  frame  contributes  to  the  accommodation  of  shock  energy,  thereby 
isolating  the  impacted  window  from  damage. 

An  additional  goal  is  to  obtain  transparency  cross  sections  possessing 
the  lowest  life  cycle  cost  that  will  meet  the  transparency/support  structure 
system  requirements.  Sierracin's  economical  and  innovative  approach  (Ref- 
erence 51)  of  retro-fitting  the  upgraded  Lockheed  P-3  and  Cessna  T-37  wind- 
shields should  be  mentioned.  The  knowledge  gained  from  these  military  pro- 
grams is  currently  being  applied  to  commercial  aircraft  windshields,  such 
as  the  Boeing  747  and  Lockheed  L-1011  at  substantial  savings  to  the  carriers. 

The  ability  of  different  glazing  materials  and  combinations  of  materials 
to  withstand  bird  strike  penetration  and  the  failure  modes  that  can  be  iden- 
tified in  case  of  damage  varies  widely  (Reference  4).  Actual  windshields, 
of  the  F-4  for  example,  are  known  to  have  experienced  initiation  of  breakage 
in  the  corners,  which  may  suggest  designing  very  stiff  corners  for  specially 
vectored  flexibility  for  stress  relief.  A suitable  approach  could  be  pro- 
perly applied  plastic  paddings  (Reference  27)  and/or  gussets. 

The  design  variables  for  bird  impact  resistance  are  thickness  of 
material,  type  of  construction,  edge  attachment  design,  and  supporting 
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structural  frame  (Reference  52)  within  the  required  geometry.  A factor 
of  safety  to  account  for  uncertainties  „.iOuld  be  applied  to  demonstrate 
enough  strength  of  the  windshield  for  thermal  stresses  since  they  occur  due 
to  thermal  differential  expansions  or  contractions.  Lack  of  rigidity  of 
the  windscreen  support  structure  generally  results  in  lower  bird  impact 
resistance.  Other  poor  design  details  which  should  be  avoided  are  sharp 
corners,  pre-existing  material  deficiencies,  and  poor  quality  of  holes. 

Too  many  laminas  are  disadvantageous  for  windshield  strength  (Ref- 
erence 33).  Windscreen  panels  laminated  with  PVB,  with  four  .24"  main 
plies  of  glass,  instead  of  two  .5"  glass  plies,  tested  at  +50eC,  indicate 
a reduction  of  failing  speeds  of  approximately  20%.  The  shatter  strength 
of  laminated  plate  is  highest  at  maximum  specific  strength  of  the  laminas 
at  optimum  stiffness  assuming  the  adhesives  in  the  laminating  process 
do  not  degrade  the  basic  materials  and  weaken  the  laminated  panel. 

The  bending  stiffness  of  laminated  plate  with  highly  plastic  inter- 
layer^) is  reduced.  Shear  load  distribution  to  account  for  shear  defor- 
mation in  addition  to  bending  deflection  should  be  considered  in  this 
type  of  laminate  arrangement. 

Proper  design  and  fabrication  of  laminates  must  consider  the  compa- 
tibility of  adherends  having  different  thermal  expansions.  Optimally  de- 
signed adhesive  edge  and  corner  bonds  with  the  adhesive/glazing  material 
interface  behavior  would  distribute  stresses  evenly.  This  would  eliminate 
fastener  holes  and  corresponding  stress  concentration  factors. 
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i 

The  major  requirement  for  cockpit  front  windows  to  resist  impact 

| 

by  bird  and  hailstones,  makes  some  of  the  organic  materials  well  suited 
to  economical  designs.  Conforming  to  this  requirement,  the  windows  ! 

not  only  stop  the  missiles  but  retain  a high  degree  of  visibility  after 
impact  (Reference  53) . A glass  facing  is  often  mandatory  for  the 
outside  surface  of  these  materials  from  the  standpoint:  of  durability 

I 

(especially  abrasion  protection)  and  positioning  of  anti-icing  films  i 

close  to  the  outside  for  effectiveness. 

The  location  of  impact  is  another  very  important  design  factor  < 

(Reference  54) . The  failure  strength  of  windscreens  with  four  main  ply  * 

laminations  appears  to  be  less  critical  for  corner  shots  than  those  with 

1 

only  two  main  ply  laminations.  This  supports  arguments  for  tailored 

i 

flexibility  (and  strength)  where  required  (Reference  41,  42).  No  attempt  ■ 

j 

should  be  made  to  design  no-deformation  characteristics  into  areas  of  ! 

probable  bird  strike.  This  would  require  immense  strength  (and  weight) 
and  it  would  not  be  practical  except  for  small  birds  and  low  impact 

speeds.  It  is  more  practical  to  provide  for  yielding  characteristics  j 

such  that  the  impact  is  deliberately  cushioned.  The  authors  of  References 
41  and  42  mention  other  good  design  practices.  Taper  angles  should  be 
selected  as  sharp  as  practical.  Surfaces  should  be  kept  perfectly  smooth 
offering  no  projection  against  which  the  bird  might  pile  up. 

The  bird  Impact  problem  is  accentuated  by  the  temperature  effects  i 

| 

due  to  kinetic  heating  of  the  interlayers  (Reference  32) , which  on  many 
present  windscreen  designs  of  sandwich-type  assembly  is  taken  into  j 
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I 

I : 

consideration.  Temperature  of  interlayers,  particularly  at  their  edges, 
where  the  panel  is  joined  to  windscreen  frame  have  to  be  kept  generally 
below  130°C,  which  is  the  upper  thermal  capability  of  organic  materials. 

t 

Proposed  and  on-going  programs  in  the  AF  Materials  Laboratory  (AFML) 
include  (Reference  1) : efforts  to  develop  a new  adhesive  and  interlayer 

enabling  the  fabrication  of  polycarbonate  composite  windows;  evaluation 
of  abrasion  resistant  coatings  for  transparencies;  and  work  on  high 

temperature  transparencies  of  polyarylsulfone  for  400-500°F  use.  | 

An  analysis  must  be  conducted  to  determine  the  thermal  windshield  I 

I 

environment  of  subsequent  design  for  measures  to  control  temperature.  J 

1 
I 

The  surface  conditions  of  windshields  at  supersonic  speed  may  be  estimated 
using  shock  tables  of  air  flowing  around  a cone.  Assuming  a cone  opening 
angle  of  22°  ( the  approximate  angle  of  inclination  of  a F-lll  windshield) , 

the  steady  flow  condition  for  a low  altitude  case  at  Mi=1.2,  ambient  air  j 

pressure,  p-j=14.7  psi,  temperature,  Ti=520#Rankine,  and  density  of  air  \ 

3 1 

§1=0. 07651  pound/ft  (density  altitude>*0) , is  as  follows.  A separated  nor-  j 

mal  shock  exists  for  these  conditions  with  P2/pi=1.513,  T2/Tj=1.128,  and  M2=  j 

0.8422.  Indices  1 and  2 indicate  conditions  upstream  of  and  after  the  shock,  j 

respectively.  At  a Prandtl  number,  Pr=.7,  for  air  at  600°R,  the  recovery  ! 

1/3 

factor  is  r**0.7  =0.888  (turbulent  case).  The  conditions  affecting  the  j 

windshield  at  this  low  altitude  case  are:  wall  temperature,  T2=198°-F, 

pressure  normal  to  surface,  p2=22.2  psi  and  (negligible)  viscous  drag  \ 

friction  of  0.026  psi.  Similar  calculations  for  a high  altitude  case  of 

I 
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Mj=2 . 5 at  50,000  ft,  pL=1.69  psi,  ^=-67°?  and  &=. 01161  pound/ft3 * *, 
indicate  a conical  shock  wave  of  34.3°  opening  angle,  and  the  windshield 
conditions:  T2=33R°F  and  P2=4.3  psi. 

Reference  32  contains  a chart  of  skin  temperature  as  a function  of 
Mach  number  and  altitude.  At  0 degree  incidence,  .9  emisslvity,  and 
2 ft  distance  from  edge,  the  skin  temperature  in  turbulent  flow  is 
fairly  constant  (100°C)  for  40,000  to  120,000  ft  altitude.  In  heat  flow 
and  thermo-stress  analyses,  discrepancies  in  the  thermal  properties  must 
be  considered.  Typical  values  for  glazing  and  frame  materials  are  : 
coefficient  of  expansion  for  soda  lime  glass  ~9.;  PVB  (+25°C)  400.; 
and  stainless  steel  ~20  in/in6C  x 10“^;  and  thermal  conductivities  of 
~15.;  ~2.5;  ~90  Btu/ft1 2/h/°C/ln,  respectively  (Reference  32). 

A full  spectrum  of  fracture  toughness  data  of  glazing  materials 
does  not  exist  (Reference  16).  The  k-value  of  MIL-P-5425  and  MIL-P- 
25690  acrylic  sheet  varies  between  1 to  4 ksi^  inciT  as  a function  of 
temperature  between  -40°F  and  150°F. 

1.  WINDOW  PANE  MATERIALS 

a.  LAMINATED  GLASS 

Basic  glazing  is  characterized  (Reference  32)  as  the  one  having  the 

best  optical  properties,  consistent  strength,  and  stiffness  over  a wide 
temperature  band,  being  unaffected  by  solvent  attacks,  and  capable  of 

maintaining  strength  when  laminated.  Glass  is  highly  sensitive  to 

fabrication  and  service  induced  scratches.  High  thermal  conductivity 
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and  low  coefficient  of  expansion  should  be  matched  to  that  of  the  metal 
retaining  frame. 

The  inherent  brittleness  of  glass  usually  relegates  it  to  a sand- 
wich construction  consisting  of  two  or  more  layers  of  glass  with  a 

flexible  plascic  interlayer  in  between.  Therefore,  the  thickness  of 
the  glass  faces  of  the  laminate  have  little  effect  on  impact  strength 
within  reasonable  limits,  while  the  plastic  interlayer  at  an  optimum 
temperature  has  a dominant  effect  (Reference  16) . Thin  glass  face  sheets 
laminated  to  acrylic  structural  plies  do  not  increase  the  penetration  speed. 
They  break  at  much  lower  speeds  than  the  calculated  failure  speed  of  the 
composite  (Reference  19).  Laminated  aircraft  glazing  is  specified  in 
MIL-G-25871A(1). 

The  average  modulus  of  rupture  of  soda  lime  toughened  glass  is  27,000 

psi,  and  of  chemically  toughened  glass  is  35,000  psi.  200°C  is  the  suggested 

maximum  temperature  for  operational  use  of  glass.  Cockpit  enclosures 
for  Mach  2.5  to  3 flight  vehicles  (370^sT  s 560° F)  of  the  non-laminated 
flat  glass  design  showed  the  disadvantages  of  poor  visibility,  heavy 
weight  (16  psf)  and  environmental  control  system  (ECS)  penalties. 

Thermally  toughened  glass  contains  a built-in  parabolic  stress 
distribution  through  its  thickness.  The  ratio  of  compression  to  central 
tension  is  2 to  1 with  the  neutral  layer  at  21 7,  of  thickness  below  the 
surface  (Reference  32).  Moduli  of  rupture  at  rates  of  loading  from  10,000 
to  60,000  psi/min  do  not  reflect  strength  behavior  in  rapid  bird  strike 
environment.  The  strength  of  flat  glass,  as  represented  by  the  modulus 
of  rupture,  is  affected  not  only  by  the  various  strengthening  treatments, 
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but  also  by  abrasion,  which  reduces  the  modulus  drastically  (Reference 
16).  In  addition,  the  measured  strength  of  a batch  of  specimens  shows 
a wide  scatter  under  identical  test  conditions. 

Failure  of  combination  glass  windows  is  unpredictable.  The  707 
windshields  are  a combination  of  .2”  semi-tempered  glass  outer  ply,  . 38,! 
hard  PVB  membrane  interlayer  with  an  electrical  conductive  coating,  and 
.5"  fully-tempered  glass  inner  ply,  sandwiched  together  and  mounted  in  a 
rigid  frame.  Loading  is  a repetitive  cycling  of  heat,  cold,  pressurization 
and  pressure  relief.  High  window  heating  is  always  required  prior  to 
descending  from  higher  cruise  altitudes.  Thermal  gradients  thru  the 
laminas  are  induced  by  temperatures  of  an  average  -40°F  in  outer  atmosphere, 
110°F  in  the  PVB  and  70°F  inside.  Holding  pressure  differential  is  8.6  psi. 

Bullet  resistant  glass  is  one  of  the  recently  developed  cockpit 
survivability  features  of  a close  air  support  aircraft  (Reference  55).  The 
bullet-proof,  bird-proof  front  windshield  successfully  defeated  a 7. 62mm 
armor-piercing  projectile  without  generating  any  backface  spallation. 

b.  FLEXIBLE  INTERLAYERS 

Interlayers  of  windshields  must  have  good  optical  properties  and 
adhesion  to  the  basic  glazing,  chemical  stability,  solvent  resistance, 
and  high  ductility  over  a wide  temperature  range  (Reference  32).  Bird 
impact  tests  of  laminated  windshields  made  with  PVB  plastic  have  shown 
impact  resistance  to  be  primarily  a function  of  the  interlayer  thickness, 
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temperature,  and  the  percent  plasticizer  content  in  the  interlayer. 

PVB  interlayers  are  most  effective  for  penetration  protection  at  21% 
plasticizer,  and  heated  to  a temperature  of  about  90°  to  100°F. 

At  lower  or  higher  temperatures,  the  impact  strength  decreases  rapidly 
(References  16  and  34).  Differences  in  the  behavior  of  glass  or  acrylic 
windows  with  thin  (1  to  2 mm)  PVB  or  cast-in-place  (CIP)  interlayers  were 
found  to  be  insignificant  in  the  temperature  range  from  0°  to  30°C  (Refer- 
ence 19).  Windows  with  PVB  interlayers  of  5 mm  thickness  had  to  be  at 
18°C  for  comparative  results. 

Thick  PVB  interlayers  play  the  part  of  an  elastic  member  at  30°  to 
40°C.  Temperature  seems  to  have  no  influence  on  glass  strength,  only 
on  PVB  strength.  The  results  of  a theoretical  study  (Reference  56) 
indicate  that  for  certain  types  of  aircraft  the  time  required  for  the 
entire  vinyl  layer  of  the  windscreen  to  reach  its  optimum  temperature  may 
be  greater  than  that  required  for  the  aircraft  to  clear  the  bird  inhabited 
zone. 

Under  favorable  conditions,  ultimate  strength  of  PVB  interlayer  material 
is  about  2.7  ksi  at  a stiffness  of  10^  psi,  and  modulus  of  rupture  is  3000 
psi  at  70°C.  Changes  of  impact  velocity  up  to  100  mph  were  experienced  by 
modification  of  the  flexibility  of  the  supporting  structure  (Reference  34). 

In  cases  of  bird-proof  w i dshield  designs  (Reference  16) , the  thickness  of 
the  PVB  interlayer  was  reduced  20%  at  the  aft  edges  and  corners  to  counter 
reduced  flexibility  there. 
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PVB  is  visually  objectionable  in  thermal  flight.  It  develops  inter- 
nal bubbles  if  exposed  to  temperatures  above  160°F  for  very  long  periods, 
or  230°F  for  a few  hours.  The  refractive  index  of  plastic  (as  in  PVB)  is 
nearly  the  same  as  of  glass  (1.5  to  1.52).  Haze  in  PVB  at  2"  thickness  is 
less  than  2%. 

c.  ACRYLIC  COMPOSITE  WINDSHIELD 

In  most  glass  windshields,  the  interlayer  is  structural,  serving  as 
the  primary  means  for  resisting  bird  impact,  and  also  as  the  fail-safe 
member  to  prevent  explosive  decompression  if  the  glass  breaks  (Reference 
46).  Where  aircraft  profiles  require  curved  panels,  plastic  materials 
become  an  obvious  choice  (Reference  16).  The  first  flat  composite 
windshield  was  delivered  to  Douglas  Aircraft  Company  as  a part  of  a 
study  of  comparative  bird-impact  strength  performed  in  1961.  In  the  case 
of  the  stretched  acrylic  plastic,  its  superior  structural  performance 
relegates  the  interlayer  to  merely  a flexible  bonding  layer. 

The  effect  of  temperature  in  the  range  -40°C  to  +50°C  on  the  failure 
speeds  of  acrylic  windscreens  is  generally  small.  The  stretched  acrylic 
materials  are  more  critical  than  glass  at  the  lower  temperatures,  whereas 
the  reverse  is  true  for  higher  temperatures.  Monolithic  acrylic  panel 
results  show  failure  speeds  approximately  15%  lower  than  multi-laminated 
panels  of  the  same  total  stretched  acrylic  thickness  (Reference  54) . When 
as-cast  and  stretched  acrylic  are  bonded  together  to  a flat  window,  they 
behave  as  if  the  window  were  completely  made  of  as-cast  acrylic  (Reference  19). 

Boeing's  first  choice  for  the  747  airplane  was  a pair  of  curved 
electrically-heated  composite  panels  with  less  aerodynamic  drag  and  cockpit 
noise  (Reference  46).  The  main  structure  of  the  windshield  is  composed  of 
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two  layers  of  stretched  acrylic  .9"  thick  providing  a fail-safe  dual  load 


path.  Flat  .05"  thick  front  and  back  surface  plies  of  chemically  streng- 


thened glass  are  flexed  to  the  windshield  for  abrasion  resistance,  adding 


bending  stiffness  to  the  assembly.  Surprisingly  small  peeling  forces 


were  found  in  these  surface  plies.  The  laminae  are  held  together  with 


3 PVB  interlayers.  The  windshield  developed  for  the  L-1011  of  the  same 


materials  may  be  considered  a novel  design  to  carry  heavy  hoop-tension  and 


shear  loads  as  a part  of  the  fuselage  surface.  At  a curvature  of  59", 


it  is  substantially  larger  (3x4  ft)  and  thinner  than  the  747  windshield. 


One  of  the  most  advanced  designs  of  stretched  acrylic  (Reference  51) 


is  the  windshield  of  the  F-5E.  It  is  a wrap-around  configuration  with  a 


shallow  viewing  angle  that  requires  an  exceptional  amount  of  care  in  pro- 


cessing. The  processing  machine  stretches  biaxially  squares  of  acrylic  from 


six  feet  to  ten  feet  while  it  slowly  cools  in  an  environmentally  controlled 


room.  This  is  followed  by  a large  scale  optical  grinding  and  polishing  pro- 


d. HIGH  STRENGTH  TOUGHENED  GLASS 


Windshield  design  without  birdproofing  interlayers  is  possible  if  the 


plies  carry’’  g the  inflight  loads  can  be  provided  with  adequate  impact 


strength.  This  is  the  case  for  stretched  acrylic  at  lower  speeds.  Chemically 


toughened  glass  could  be  used  (Reference  32)  if  it  can  be  developed  to  a 


modulus  of  rupture  of  50,000  pal.  The  latest  types  ot  toughened  glass 


have  gone  a long  way  towards  reducing  the  inherent  weight  advantage  of 


plastics  (Reference  57) . Triplex  10-20  high  strength  glass  is  approximately 


a 


n 
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30%  more  resistant  to  bird  impact  than  thermally  toughened  glass  (Reference 
54) . The  relationship  between  failing  speed  and  total  glass  thickness 
was  found  to  be  linear  and  independent  of  the  main  ply  composition. 

The  technique  of  defeating  a bird,  using  inclined  heavy  gage  glass 
panels,  could  be  called  bird-bouncing.  The  bird  is  deflected  off  the 
windscreen  and  the  debris  retains  a high  proportion  of  the  initial  energy 
(Reference  18) . This  is  the  design  philosophy  for  glass  windscreens  as 
defined  in  Reference  58.  The  prime  function  of  the  glass  is  to  resist 
bird-impact  and  pressure  loads.  When  two  or  more  sheets  of  heat-tempered 
plate  glass  are  bonded  together  by  flexible  interlayers  of  PVB,  the  PVB 
provides  a reserve  load  for  pressurization.  The  layers  of  glass  resist 
loads  through  bending,  with  the  frame  around  the  panel  acting  as  a simple 
support.  It  is  designed  so  that  airframe  loads  are  not  transmitted  through 
the  transparency.  The  high  modulus  of  rupture  required  for  the  impact 
resistance  is  achieved  by  artificially  introducing  a stress  distribution 
in  the  glass  (Figure  37) . In  the  unloaded  state  the  center  is  in  tension 
while  the  surfaces  are  in  compression.  A large  bending  load  has  to  be 
applied  before  the  curface  changes  from  compression  to  tension.  The 
annealed  strength  of  the  glass  has  to  be  added  to  the  compressive  stress 
to  calculate  the  moduli  of  rupture  in  bending. 

Triplex  has  gained  acceptance  for  manufacturing  curved  cockpit  glazing 
for  Boeing  747  windscreens.  It  Is  rated  for  impact  resistant  to  a 4 pound 
bird  at  427  kt  (Reference  58) . The  windshields  are  composed  of  two  panes  of 
thermally  toughened  "10-20"  glass,  each  12  mm  thick.  In  the  A300B  Airbus  pro- 
duction (Reference  57)  flat  "10-20"  windscreens  will  be  used.  They  will  be 
composed  of  two  main  plies,  each  10  mm  (.394")  thick,  and  a facing  ply  of 
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flat  .05"  Chemcor  glass.  Corner  bird  strike  under  low  temperature  soak 
conditions  is  to  be  considered  critical  for  the  A300B  windshield. 

e.  POLYCARBONATE  PLASTIC  SHEET 

The  most  outstanding  property  of  polycarbonate  is  its  ductility.  The 
stress-strain  curve  shows  elastic-perfect  plastic  behavior  with  small 
instability  at  initial  yielding.  Its  low  abrasion  resistance,  however, 
requires  surface  plating  with  other  materials.  As  a consequence,  the 
impact  strength  of  fusion-clad  polycarbonate  is  much  lower  than  that  of 
unclad  polycarbonate  (Reference  46).  This  is  caused  by  the  brittle  nature 
of  the  as-cast  acrylic  cladding  and  propagation  of  cracks  through  the 
fusion  bonded  surface  into  the  polycarbonate.  Sierracin  developed  a 
manufacturing  process  to  overcome  this  problem.  Instead  of  cladding,  a 
thin  glue  line  (.01")  of  flexible,  high  temperature  resistant  silicon 
resin  bonds  the  polycarbonate  with  thin  and  flexible  .01"  - .025"  inter- 
layer acrylic  sheet. 

The  "230  series"  of  ductile  protective  coatings  provide  increased 
weathering  resistance  of  polycarbonate  as  compared  to  hard  Sierracote 
220  coating  (Reference  59).  Specimens  were  impact  tested  with  a specially 
shaped  40  pound  weight  of  1"  diameter  from  a height  of  7,5'.  Impact 
strength  was  not  adversely  affected  by  the  thin  coating  in  the  range  of 
-65°F  thru  room  temperature. 

Another  peculiarity  of  polycarbonate  has  been  found.  Notched  Izod 
impact  strength  vs.  material  thickness  data  (Reference  60)  indicates  exis- 
tence of  a sharp  ductile  to  brittle  transition  in  polycarbonate  at  a 
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sheet  thickness  of  about  .18".  The  measured  specific  strength  was  16.8 
ft-pound/in  for  a specimen  thickness  of  .125",  and  only  3.9  ft-pound/in 
for  .1875".  For  this  reason,  it  is  desirable  that  polycarbonate  in  one 
eighth  inch  sheets  be  used  as  the  primary  structure. 

Polycarbonate  sheet  is  the  most  heat  resistant  transparent  plastic 
currently  available  for  production  programs.  Internal  bubbling  defines 
the  upper  limit  for  operating  temperature  and  occurs  after  moderate 
(short)  exposure  at  315°F  (350°)  (Reference  46).  Polycarbonate  unfortu- 
nately, exhibits  low  light  transmission  and  high  haze  content.  Quality 
surface  polishing  of  extruded  sheet  has  to  be  extensively  pursued  during 
the  windshield  manufacturing  process  and  is  difficult  to  achieve.  As  a 
result,  a change  in  the  prime  specification  to  a minimum  of  acceptance  of 
65%  light  transmission  has  been  recommended  by  Rockwell  International 
(Reference  44) . 

The  most  suitable  grade  polymer  found  for  aircraft  glazing  is  flat  sheet 
"3000  S.L.".,  ultra-violet  stabilized,  water  white  grade  polycarbonate  from 
Farbenfabriken  Bayer  A.G.  (Reference  61).  Mechanical  properties  are  identi- 
cal to  those  of  the  other  grades,  Fty=64.8  MN/m^  (1  psi  = 6895  N/m^)  at 
20°C,  and  60%  reduction  at  maximum  operating  temperature  of  140°C.  Stressed 
and  unstressed  polycarbonates  were  found  resistant  to  many  aircraft  fluids, 
but  were  immediately  attacked  by  Skydrol  500  and  bromochloromethane.  MIL-P- 
83310  is  the  military  specification  for  transparent  polycarbonate  plastic  sheet. 

In  one  test  series,  the  response  to  bird  impact  of  27  different  flat 
panels  of  polycarbonate  (and  acrylic)  windshield  materials  was  evaluated  for 
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impact  energy  absorbing  potential  at  non-penetration  vs.  pane  thickness  and 
angle  of  impact  (Reference  62) . In  another  impact  experiment  with  a 4 pound 
bird  at  a Canadian  test  site  (Reference  63) , a 36"  x 36"  curved  panel  of 
polycarbonate  in  a representative  windshield  frame  survived  corner,  center 
and  edge  shots  at  a top  speed  of  650  mph  (<v500  kt).  Composition  of  the  panel 
was:  0.125'1  as-cast  acrylic  outer  layer  and  ".1"  plastic  interlayer  dis- 

continued at  the  edges;  0.875"  polycarbonate  main  structural  ply,~.l" 
plastic  interlayer;  0,15"  polycarbonate  splinter  shield.  Four  (4)  other 
concepts  and  one  (1)  of  polycarbonate,  were  tested.  In  all  penetration 
occurred. 

In  a more  recent  test,  a two  ply/interlayer  polycarbonate  pane 
(3'  x 4'  and  1"  total  thickness)  was  tested  at  140°F  with  a 4 pound  bird 
at  450  kt.  The  specimen  showed  an  impression  due  to  permanent  plastic 
deformation  approximately  2 1/2"  maximum  in  depth  and  somewhat  larger  than 
the  projected  area  of  bird  (Reference  64). 

Rockwell  International,  Los  Angeles  Division,  conducted  an  evaluation 
of  glazing  materials  for  the  B-l  (Reference  44) , and  selected  monolithic 
polycarbonate  as  the  prime  candidate,  obviously  due  to  its  exceptionally 
high  impact  resistance  over  a wide  range  of  temperature,  and  its  strength 
at  high  temperature  (Ftu  = 5.2  ksi  at  220°F).  The  main  right  and  left  wind- 
shields of  the  B-l  have  the  geometry  of  a single  curved  truncated  cone  with 
approximately  a 50"  average  radius,  covering  an  area  of  2500  ln^  each,  with 
a ratio  of  width/height  of  2,  and  angle  of  inclination  approximately  23°  to 
25°.  Swedlow  Inc.  has  been  selected  as  subcontractor  (Reference  65)  to 
develop  these  bird-proof  windshields  for  speeds  over  650  mph.  The  operational 
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temperature  requirement  is  in  the  -40  to  125°F  temperature  range  and  requires 
a thermal  gradient  study  to  be  conducted.  The  structural  build-up  of  the 
laminate  in  sequence  of  arrangement  is:  (1)  .09"  abrasion  resistive  acrylic 

ply  at  the  outside,  (2)  vacuum  deposited  gold  film  underneath  for  de-icing 
and  radar  sensitivity,  (3)  .16"  cast-in-place  silicon,  (4)  .87"  polycar- 
bonate (the  load  carrying  ply),  (5)  .05"  cast-in-place  silicon  layer,  and 
(6)  .15"  polycarbonate  inside  layer.  This  lamination  totals  1.398"  thick- 
ness, providing  70%  transparency  at  3%  haze.  The  side  and  upper  windows 
are  claimed  (Reference  66)  to  be  bird-proof  and  temperature  resistant.  Side 
windows  of  the  B-l  will  be  a laminated  structure  of  polycarbonate  and 
stretched  acrylic  plastic,  while  the  upper  windows  - above  the  crew  - will 
be  monolithic  polycarbonate.  Special  abrasion  resistant  coatings  will  be 
applied  to  the  polycarbonate  exteriors. 

A more  recent  report  on  the  B-l  windshield  (Reference  15)  states  that 
it  is  by  far  the  largest  aircraft  project  ever  undertaken  and  the  most  tech- 
nical demanding.  The  composite  transparency  consists  of  two  polycarbonate 
plies,  two  cast-ln-place  silicone  interlayers,  and  a glass  outer  ply.  It 
contains  a Sierracote  anti-ice  heating  element  and  ''oatings  to  provide  radar 
reflection  and  static  drain.  The  use  of  normal  interlayer  materials  is 
precluded  by  the  operating  requirements  of  the  B-l,  particularly  the  high 
temperatures  generated  by  near-supersonic  speeds  at  low  altitudes,  which 
let  to  new  processing  techniques.  Exceptional  optical  requirements  are  im- 
posed on  the  windshield  by  its  shallow  viewing  angle.  It  is  further  stated, 
that  the  light  weight  fighter  YF-16  canopy  is  the  largest  and  deepest  drawn 
polycarbonate  component  ever  formed. 
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2.  FRAME  DESIGN  DATA 

Edge  attachment  designs  (Figures  38  and  39)  have  become  more  complex 
as  the  demands  of  enclosure  requirements  have  increased  (Reference  16) . 

The  development  of  the  method  of  fixation  and  the  rigidity  of  supports 
for  a particular  aircraft  make  it  impractical  to  specify  a certain  design 
for  a given  set  of  parameters. 

A manufacturer’s  handbook  on  glass  aircraft  windshields  (Reference  67) 
states  that  some  glass  parts  with  no  edge  attachment  are  mounted  in  the  air- 
plane with  a rubber  channel  similar  to  that  used  in  automobiles  some  years 
ago.  Others  consist  of  metal  inserts  with  molded  vinyl  and  glass  fabric- 
phenolic  reinforcement,  metal  "Z"  frames  and  metal  retaining  rings.  The 
complexity  of  the  edge  attachment  is  a function  of  the  structural  and  oper- 
ating design  parameters  of  the  individual  aircraft.  One  common  mounting 
design  concept  for  glass  is  "slip  plane  parting  medium"  to  reduce  delamina- 
tion, edge  chipping  caused  by  temperature  and  pressure  induced  stresses.  It 
distributes  the  stress  over  a larger  area  and  minimizes  the  possibility  of 
failure  due  to  point  stress  loading. 

Conflicting  recommendations  were  found  concerning  the  effect 
of  frame  stiffness  upon  the  bird  impact  resistance  of  windshields.  A 
Civil  Aviation  Agency  (CAA)  report  (Reference  16)  states  that  windshields 
mounted  in  more  flexible  supporting  structures  resisted  penetration  at 
velocities  about  100  mph  higher  than  did  identical  panes  mounted  in  the 
more  rigid  cockpit  structure. 

In  Reference  18,  the  problem  is  seen  as  follows:  the  mounting  attach- 
ment and  frame  must  transmit  the  impact  load  in  such  a way  that  undue  local 
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Figure  39.  Calculated  Static  Flight  Load  Deflections 


for  Different  Mounting  Concepts 


161 


HOOP  TENSION  WITH 
BUILT-IN  SILL 
RESTRAINT 


SUPPORTED  WITH 
I"- RILL  RESTRAINT 


V 


AFFDL-TR-74-155 


stress  concentrations  do  not  result  in  failure  of  the  screen.  Similarly, 
the  frame  should  be  able  to  accept  these  loads  without  local  deformation 
causing  stress  concentration  in  the  screen.  A recent  experimental  investi- 
gation (Reference  53)  indicated  that  windscreens  should  be  mounted  in  a 
frame  as  rigid  as  possible  to  obtain  the  highest  impact  resistance.  This 
is  very  important  for  thermally  toughened  glass  windscreens  which  were 
shown  to  be  particularly  sensitive  to  reduction  in  frame  stiffness.  In 
Reference  58,  however,  use  of  the  frame  to  absorb  energy  thru  distortion 
during  the  bird  strike  is  cited. 

The  points  of  contact  between  the  glazing  material  and  the  enclosure 
frame  have  to  be  designed  as  load  bearing  to  transmit  all  principal  in- 
flight loads,  and  at  the  same  time  non-load  bearing  to  accommodate  the 
thermal  expansion  of  the  glazing  material.  Reference  16  concludes  that  the 
glazing  material  should  be  securely  fixed  within  the  supports,  but  should 
be  free  to  expand  or  contract  due  to  changes  in  temperature  and  with  aging. 
The  coefficient  of  thermal  expansion  of  the  frame  should  be  matched  as 
nearly  as  possible  to  that  of  the  glazing  (Reference  32) . If  this  is  not 
possible,  the  frame  design  should  be  such  that  compression  stresses  are 
created  in  the  transparency. 

The  edge  attachment  for  glass  is  different  from  that  for  a plastic 
glazing  material.  Reference  46  states  the  careful  isolation  of  flat  glass 
windshields  from  the  fuselage  loads  (Figure  38.1)  is  required.  Plastic 
windshields,  especially  the  curved  ones,  can  be  bolted  directly  into  the 
metal  structure  (Figure  38.2)  so  chat  cabin  pressurization  loads  are  carried 
by  hoop  tension,  while  bird  impacts  are  absorbed  by  the  windshield  acting 
as  an  arch  in  compression.  This  permits  the  windshield  and  the  surrounding 


i 


AFFDL-TR-74-155 


cabin  structure  to  be  loaded  structurally  in  the  most  efficient  manner. 

The  747  cockpit  structure  does  not  take  full  advantage  of  this  "bolt-in" 
capability,  but  the  Lockheed  L-1011  does.  Monolithic  polycarbonate 
one-half  inch  thick  with  a bolt-through  attachment  has  had  success  for 
particular  aircraft  applications  (Reference  16).  A clamped  configuration 
gives  generally  better  results  than  a bolted  one  (Reference  19) . A well 
designed  bolted  attachment,  however,  can  give  practically  the  same  results 
as  a clamped  design  solution. 

The  method  of  fitting  and  attachment  can  present  maintainability  pro- 
blems (Reference  68)  because  of  residual  strain  due  to  initial  fitting  or 
subsequent  change.  For  this  reason,  windshield  panels  are  designed  and 
developed  already  mounted  in  their  frame  for  use  in  the  Royal  Air  Force 
The  windshield  assembly  of  the  747  airplane  is  similarly  preinstalled  in 
an  aluminum  frame  (Reference  46).  Installation  into  the  aircraft  is  a 
simple  metal-to-metal  bolt-up  with  dry  pressure  seals,  performed  entirely 
from  outside  the  airplane.  Use  of  edge  attachments  (a  means  of  fastening 
the  side  edges  of  the  glazing  to  the  enclosure  and  consisting  of  side  beams, 
transverse  frames  and  other  reinforcements)  may  easily  lead  to  greater  stress 
concentrations  than  those  the  edge  attachments  were  intended  to  eliminate 
(Reference  16).  To  minimize  these  stress  concentrations,  stresses  in 
the  materials  should  be  uniformly  distributed.  In-flight  experimental 
testing  using  strain  measurements  of  instrumented  panes  and  frames  of  cock- 
pit windows  is  one  way  of  achieving  this  (Reference  14) . 

Chapter  5 of  Reference  69  provides  information  on  fastening  and  joining 
techniques  of  reinforced  plastics.  It  also  furnishes  data  on  the  effects 
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of  the  use  of  different  materials  and  joint  parameters  upon  failure  modes 
and  fatigue  strength.  The  consequence  of  bolting  alone  can  be  seen  in 
higher  stress  concentration,  vs.  a greater  risk  of  catastrophic  operational 
failure  in  case  of  adhesive  bonding.  The  following  lap  joint  data  for 
0.063  in.  thick  7075-T6  clad  are  selected  to  indicate  two  joining  methods; 
however,  they  do  not  represent  typical  glazing/frame  material  joining 
characteristics.  Lap  joint  strength  for  3/16  in.  diameter  tapered  inter- 
ference fasteners  at  a spacing  of  4 diameters , is  1700  pound/in.  The  stress 
concentration  factor,  kc,  is  3.  Lap  joint  strength  of  an  adhesive  bond  at 
a length/thickness  ratio  of  40  is  1900  pound'/ in;  kt=1.5  (Reference  70). 

The  following  are  a few  results  of  efforts  undertaken  to  provide 
improved  windshield  protection  for  the  F-lll.  After  a series  of  tests  at 
AEDC/Tullahoma  TN  using  a F-lll  crew  escape  module  with  the  original  wind- 
shields replaced  by  several  monolithic  polycarbonate  windshields  and  those 
of  modified  laminated  polycarbonate  Pittsburgh  Plate  Glass  (PPG)  proprietary 
plastic  interlayers,  the  damage  to  the  windshield  frames  was  assessed. 

The  part  of  the  ledge  on  the  right  hand  arch  (which  had  completely 
failed  during  one  of  the  tests)  supporting  the  windshield  behind  the  impact 
location,  was  fractured  in  a combination  of  shear,  bending  and  tension.  The 
rest  of  the  ledge  failed  primarily  in  shear  except  for  4-5  inches  along  the 
bottom  of  the  arch.  The  arch  had  a slight  rearward  set.  The  center  post 
had  received  moderate  deformation  due  to  a direct  bird  impact.  As  a 
result,  a decision  was  made  to  reinforce  the  windshield  support  structure 
(Reference  47).  A static  stress  analysis  was  conducted  employing  a NASTRAN 
computer  model  of  the  F-lll  windshield/frame  arrangement.  The  applied  pressure 
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load  simulated  the  bird  strike  test  condition  FM-2  of  Reference  35.  In 
two  computer  runs,  a structural  configuration  with  a six  inch  corner  radius 
added  to  the  upper  aft  corner  of  the  windshield  arch  and  beam  intersection 
was  compared  with  the  original  design.  The  addition  of  the  radius  resulted 
in  a reduction  of  20%  in  the  peak  tension  stresses  in  the  mold  lint  surface 
of  the  edge  of  the  transparency. 

Structural  failures,  as  they  occurred  on  the  bird  impacted  F-lll 
windshield/canopy  with  modified  laminated  windshield  and  the  reinforced  aft 
arch  manufactured  by  McDonnel  Aircraft  Company,  was  as  follows.  Two  arch 
reinforcing  members  (1  right  and  1 left)  were  tested  first  as  designed  and 
manufactured  (Reference  71).  Attachments  indicated  permanent  set,  such  as 
flattening  and  torsional  set,  in  the  aft  arch  windshield  mounting  flange. 

Cast  aluminum  corner  fittings  on  the  ait  arch  were  broken.  Bolts  sheared 
and  kept  Che  mounting  flange  from  failing  at  that  location.  The  badly  damaged 
windshield  and  aft  arch  both  held  together  and  no  material  penetrated  the 
cockpit  in  these  tests.  It  was  decided  that  the  reinforcements  were  too 
strong.  The  amount  of  support  provided  by  the  arch  reinforcing  members  was 
reduced  by  machining  of  the  members  and  by  use  of  a different  fastener  pattern. 
There  was  no  bird  penetration,  but  the  arch  and  reinforcing  member  were 
too  badly  damaged  to  allow  installation  of  further  windshields  for  testing. 
Work  was  continued  with  further  arch  modifications.  The  results  showed 
better  bird  resistance  for  the  F-lll  in  terms  of  penetration  at  the  pre- 
established  goal  of  bird  strike  performance. 

3.  RELEVANT  LEADING  EDGE  DATA 

Non-penetration  is  recommended  as  the  design  requirement  for  leading- 
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edge  structures  for  bird-strike  resistance  (Reference  72).  Bird  impact 
resistance  is  proportional  to  leading-edge  skin  thickness.  An  empirical 
design  formula  was  developed  based  upon  the  penetration  velocity  of  typical 
designs  using  L73  aluminum  alloy  and  riveted  construction.  This  velocity 
was  measured  as  a fu;  ction  of  bird  weight  (up  to  4 pound),  skin  thickness, 
nose  radius,  rib  thickness  and  pitch,  and  sweep  angle.  Conventional  butt 
joints  reduce  the  penetration  velocity  by  15%.  Reduction  of  skin  thickness 
between  rib  or  joint  attachment  lines  (up  to  15%)  by  means  of  chemical 
etching  or  machining  is  ineffective  in  reducing  penetration  velocity. 

Designs  with  other  skin  materials,  such  as  L72  aluminum  alloy,  nickel 
and  titanium  alloys,  a CFRP-L73  alloy  laminate,  and  bonded  construction, 
were  also  investigated.  For  high  strength  titanium  and  nickel  alloy  leading 
edge  skins  on  adequately  strengthened  support  structures,  the  penetration 
velocity  is  given  by  factor  0.8  x tensile  strength  of  alloy/tensile  strength 
of  L73.  A 1.1  mm  thick  reinforcement  of  the  carbon  composite  inside  the 
leading  edge  resulted  in  the  laminated  skin  being  much  less  resistant  to  bird 
impact  than  a light  alloy  skin  of  equivalent  weight.  Bonded  nose  construction 
using  the  cold  curing  epoxy  Araldite  (USA  process  specification  DHA455) 
increases  the  impact  resistance. 
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APPENDIX  II 

SURVEY  OF  EXPERIMENTAL  DATA 
1.  BIRD  STRIKE  EXPERIMENTS 

A wealth  of  data  from  experiments  with  bird-impacted,  mounted  window 
panes  exists.  For  example,  Table  XVIII  summarizes  data  for  stretched 
acrylic  composite  windows.  These  data,  from  different  sources,  may  be 
interpolated  for  comparison  purposes  to  account  for  single  parametrical 
changes  such  as  larger  or  thicker  panes,  different  attachment  conditions, 
higher  impact  velocities,  and  heavier  birds. 

The  available  experimental  data  is  discussed  in  the  following  subsec- 
tions. 

a.  BIRD  G’JNS,  TEST  FRAMES,  INSTRUMENTATION 

Performance  capability  of  a Canadian  compressed  air  bird  gun  with  a 
6 in.  diameter  barrel  is  250  to  700  ft/sec,  with  the  ability  to  extend  the 
velocity  to  1000  ft/sec  using  4 to  8 pound  birds  (Reference  4) . A consid- 
erable increase  in  the  speed  potential  of  the  Royal  Aeronautical  Establish- 
ment Farnborough  (R.A.E.)  gun  has  been  gained  by  sealing  the  muzzle  of  the 
barrel  with  a frangible  ’Melinex'  diaphragm  and  evacuating  the  air  (Reference 
12) . Tests  referred  to  in  Reference  19  were  performed  with  a bird  gun  at 
the  Centre  d'Essais  Aeronautique  de  Toulouse  (C.E.A.T.).  R.A.E.  systematically 
tested  25"  x 19"  panels  with  rounded  corners  (Reference  12) , while  the 
French  Air  Ministry  used  a typical  frame  of  52C  mm  by  820  nun  (Reference  19) . 

Every  actual  windscreen  installation  for  bird  impact  resistance  should 
be  validated  by  tests  under  the  required  operational  conditions  (Reference 
54).  Force-time  relations,  as  they  occur  during  in-flight  bird  strike,  can 
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be  simulated  using  air  or  powder  guns  instead  of  more  expensive  rocket  sled 
tests.  Each  source  of  impact,  however,  requires  calibration  (Reference  74). 
Experience  has  shown  that  the  projectile  speed  achieved  is  within  3%  of  that 
predicted  from  the  calibration  graphs  (Reference  12) . The  means  of  visual- 
izing what  is  taking  place  during  an  impact  is  a very  important  area  of 
experimental  investigation.  Interferometric  techniques  using  photomechanics 
and  pulse  and  continuous  wave  lasers  are  used  to  determine  impact  resistant 
structural  design  features, 
b.  MOTION  PICTURE  OBSERVATIONS 

A list  of  phenomena  was  prepared  after  a review  of  Douglas  DC8/DC9/DC10 
airplanes  bird  impact  testing  movies,  photographs,  and  test  results  (Refer- 
ence 75) . The  films  reviewed  showed  the  impact  of  4 pound  birds  at  300-450 
kt.  The  bird  specimen  ruptured  immediately  upon  impact.  The  bird  either 
bounced  or  bagged,  if  there  was  no  penetration.  Blood  traces  radiated 
from  the  point  of  impact.  In  oblique  impact,  the  trace  of  splash  radiated 
from  the  point  of  impact  in  all  directions.  The  intensity  of  the  trace  was 
heavy  along  tangential  velocity  direction  and  reduced  gradually  in  the  oppo- 
site direction.  The  trace  for  the  point  of  impact  was  smaller  than  the 
packaged  cylinder  diameter  of  bird. 

A radial  crack  pattern  was  observed  for  glass  laminates,  while 
the  pattern  for  stretched  acrylic  and  chemically  strengthened  glass  wind- 
shields radiated  in  all  directions.  The  crack  pattern  in  the  tangential 
direction  was  centered  at  the  point  of  impact  and  extended  concentri- 
cally to  a finite  radius.  Inside  this  radius  there  was  permanent  defor- 
mation. The  crack  pattern  was  similar  on  both  inner  and  outer  glass 
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face  sheets.  The  cracks  radiating  approximately  normal  to  the  permanently 
deformed  circular  shapes  were  apparently  not  affected  by  the  windshield 
edges. 

From  these  observations  the  following  hypotheses  were  developed  (Reference 
75) . Bird  impact  at  certain  speeds  is  isothermal  with  no  entropy  change 
during  impact.  Total  energy  generated  from  the  impact  dissipates  in  splash- 
ing the  bird  and  converts  into  internal  strain  energy  of  the  windshield. 

The  material  conforms  according  to  its  nature:  linear  elastic-.-perfect 

plastic,  or  work  hardening.  Momentum  is  transferred  into  the  windshield 
instantaneously.  The  windshield  subjected  to  an  initial  velocity  after 
impact  has  no  instantaneous  displacement  response.  After  that  it  is  free 
from  surface  traction, 
c.  EMPIRICAL  STRENGTH  DESIGN  FORMULAE 

Two  types  of  design  equations  have  been  derived  from  bird  strike  tests. 

One  determines  structural  performance  and  the  other  estimates  impact  forces 
(Appendix  Il-ld) . 

The  empirical  equations  for  structural  performance  (References  10,  12 
and  19)  give  the  relationships  between  the  thickness  of  monolithic  glass 
or  plastic  panel,  the  bird  weight,  the  aircraft  speed,  and  panel  mounting 
angle  with  the  trajectory.  Mott  (Reference  12)  in  his  equations  further 
considers  clamping  width,  main  ply  thickness,  main  ply  composition,  temp- 
erature and  impact  position.  Poullain  and  Clamagirand  (Reference  19) 
investigated  shape,  size,  and  contour  characteristics  of  the  windows,  and 
found  that  a reduction  in  the  window  radius  would  appreciably  increase  the 


170 


AFFDL-TR-74-155 


critical  impact  speed  of  the  curved  panel. 

Similar  equations  were  developed  for  the  thickness  of  the  vinyl  inter- 
layer of  a laminated  panel  as  well  as  for  newer  materials.  The  variation 
of  impact  damage  (assuming  failure  speed  proportional  to  the  reciprocal  of 
the  cosine  of  the  sweep  back  angle  (Reference  12)),  deviated  from  Newton's 

O 

law  (m  • V ) for  test  data  recorded  from  thermally  toughened  glass  panels 
of  .8"  total  glass  thickness  according  to  m . V3  (Reference  18).  The  formulae 
for  the  rectangular  shape  of  Reference  19,  are  different  for  different  mate- 
rials and  use  and  (weight  of  bird)^/3  for  any  glass  and  acrylic  material. 

As  the  failure  analysis  of  bird-impacted  windshield  was  mostly  experi- 
mental, empirical  design  formulas  were  established  based  on  the  test  data. 

One  of  the  early  formulas  for  toughened  glass  is 

t = .136  (1  - .348  cosoO  e(V> cos  (*/873)  (75) 

where  t,  V,  and  o<  are  the  thickness  in  inches,  penetration  speed  in  mph 
for  a 4 pound  bird,  and  windscreen  angle  between  direction  of  impact  and 
normal  to  windshield,  respectively.  This  formula  is  valid  for  clamp- 
mounted  panes  ranging  in  thickness  from  1/4  to  1 1/4  inch,  and  c<  =40°  to 
60°  (Reference  18).  At  an  angle  of  50°,  500  (600)  mph  would  require  a 
thickness  of  4 (8)  inches  of  toughened  glass  to  avoid  penetration  of  a 4 
pound  bird. 

Tests  on  transparencies  (Reference  18)  have  indicated  that  penetration 
velocity  is  proportional  to  one  over  the  cube  root  of  the  bird  mass,  M.  The 
proportionality  constant  for  a vinyl  sandwich  panel  specific  thickness  under 
identical  conditions  is 

MV3=5.1  • 107  (76) 
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and  for  stretched  acrylic  it  is  6.75  • 10°  with  M varying  from  1 to  8 pounds 
and  V in  mph.  The  design  formula  for  20%  plasticized  PVB  for  the  total  PVB 
thickness  is 

t = .05  • e<v-cosK/180)  (77) 

Vcoso<  varies  only  5%  at  a range  of  <X=41°  to  46°.  Penetration  resistance 
decreases  substantially,  as  the  temperature  departs  from  the  optimum.  Scatter 
of  data  occurs  due  to  mounting  conditions  (Reference  18).  In  one  case, 
identical  screens  failed  around  the  edge  attachment  at  376  mph  in  December, 
and  at  460  mph  at  63°F  outside  temperature  in  July.  The  main  vinyl  layer  was 
at  40°C  in  both  tests.  The  design  formula  obviously  does  not  include  the 
temperature  of  the  edge  attachment  and  the  quality  of  mounting  design. 

A laminated  stretched  acrylic  panel  with  PVB  interlayers  of  .875"  total 
thickness,  inclined  at  45°,  defeated  a 4 pound  bird  at  320  mph  (Reference 
18) . This  speed  is  some  20-40%  higher  than  that  for  toughened  glass  at  the 
theoretical  glass  speed  of  265  mph.  The  high  penetration  speeds  of  the 
laminated  stretched  acrylic/polyvinyl  composite  panel  reflect  optimized 
mounting  attachments. 

rt  has  been  observed  over  the  years  (Reference  54)  that  the  impact  energy 
of  the  bird  to  produce  similar  damage  is  a constant  value,  and  that  the 
failing  speed  is  proportional  to  the  reciprocal  of  the  cosine  of  the  sweep- 
back  angle.  For  thicker  windscreens,  the  above  formulas  predict  failure 

J 

speeds  below  the  speeds  achieved  in  clamped  edge  panel  tests.  If  compared 
on  the  basis  of  weight  and  not  thickness,  the  stretched  acrylic  windscreens 
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would  be  twice  as  resistant  as  those  of  thermally  toughened  glass,  and  1 1/2 
times  as  resistant  as  Triplex's  "10-20"  high  strength  glass. 

The  following  formulas  are  for  thermally  toughened  glass,  stretched 
acrylic,  and  "10-20"  glass,  respectively  (Reference  54). 

V(mph)  = 400  • t ( in) /(-/w (pound )'*coso<)  (78) 


V(mph)  - 470  in )/(Yw(pound)'* cosoO  (79) 

V(mph)  = 520  • t(in)/(yw(pound)‘*cos«.)  (80) 

(W  = bird  weight,  t=total  thickness  of  specified  material) 


These  equations  do  not  reflect  support  lengths.  Panels  of  25"  x 19"  and 
32"  x 20"  were  tested.  No  information  is  presently  available  that  reflects 
the  effects  of  panel  size  and  shape  in  such  an  empirical  equation.  Similar 
equations  were  derived  independently  (Reference  19)  and  are  reported  for 
comparison: 

(i)  Laminated  Glass  of  several  plies  with  thin  layers  of  PVB  in  between 
(1  to  3 mm),  tested  at  ambient  temperature: 


e3/2c_v2  p2/3  ginrt . \ 

If  all  the  plies  do  not  have  the  same  tempering: 


e3/2( 


elcl+e2c2+’”  9 

) = vz 

e 


p3/3  sin0(  . \ 


(81) 


(82) 


where  e=total  thickness  of  glass  in  mm;  e1=thickness  of  ply  with  coefficient 
cl>  e2= thickness  of  ply  with  coefficient  C2  (thickness  of  PVB  not  taken  into 
account);  C=strength  coefficient  for  bird  impact  complying  with  specification 
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FAR  25  (no  penetration  of  bird  or  shattering). 

(ii)  Stretched  Acrylic  windows; 

e2  C = v2  P2/3  sin of • A (83) 

where  e ■ total  thickness  without  interlayer;  A = 1 for  civil  aircraft;  .5 
for  side  windshield  of  fighters;  .35  for  wrap  around  fighter  windshield. 

(iii)  As-cast  Acrylic  flat  windows  of  civil  aircraft  and  front  windshield 

of  a fighter:  The  equation  is  the  same  as  for  stretched  acrylic,  but  with 

C = 28,  A =1.  Curved  side  windows  of  fighter  aircraft: 

48  e2  = v2  P4/3  sin«  (84) 

for  25"  to  35"  curvature. 

The  formulae  fit  520  mm  x 820  mm  test  frame  size.  At  edge  rebates  not 
too  deep,  e is  the  pane  thickness.  Otherwise,  the  speed  is  very  similar  to 
that  for  the  thickness  in  the  rebate  area. 

A typical  requirement  for  helicopter  windshields  is  to  defeat  bird  impact 
in  the  velocity  range  of  200  to  250  mph  (Reference  76).  Specimens  manufac- 
tured and  tested  by  Pittsburgh  Plate  Glass  of  .25"  stretched  acrylic  and  the 
light  weight  glass  configuration  failed  at  impact  speeds  as  low  as  100  mph, 
when  impacted  with  a 4 pound  bird  at  room  temperature.  In  comparison,  the 
penetration  velocity  for  .25"  monolithic  polycarbonate  is  approximately  250 
mph.  It  is  the  same  for  composite  designs  of  glass  and  polycarbonate  weighing 
4.5  pound/ft2.  Table  XVII  provides  test  data  from  different  sources  for 
stretched  acrylic  composite  aircraft  windshields. 

The  resistance  to  bird  impact  of  the  modified  2'  x 2'  size  windshield 
of  the  Canadian  Air  Force  CL-41  Tutor  aircraft  has  been  investigated  (Ref- 
erence 77)  utilizing  the  high  speed  rocket  sled  facility  at  Holloman  AFB,  NM. 


174 


AFFDL-TR-74-155 


Two  (2)  layers  of  vinyl  (25  parts  plasticizer),  .18"  and  .12"  thick  were 
added  to  the  original  .45"  laminated  glass  windshield  (inner  and  outer  layers 
semi-tempered,  center  tempered  glass)  and  attached  to  the  frame.  Prior  to 
modification,  the  main  glass  yielded  a penetration  velocity  of  only  127  mph 
for  collision  with  a 4 pound  bird.  At  the  indicated  air  speed  of  298  mph 
(IAS),  a 3.5  pound  chicken  did  not  penetrate  the  right-hand  panel,  which 
was  at  a temperature  of  70°F.  Penetration  did  occur  at  this  speed  into  the 
left  hand  panel,  which  was  cooled  to  a temperature  30°F  below  that  for 
maximum  strength  of  the  vinyl.  Using  the  M.E.  106  formula  for  vinyl,  the 
theoretical  penetration  speeds  for  .12",  .18"  and  .3"  thicknesses  are  210, 

290  and  415  mph,  respectively.  This  reveals  that  the  add-on  design  modifi- 
cation partly  utilized  the  material  strengths  of  the  PVB. 

The  following  example  indicates  how  misleading  windshield  designs  for 
high  impact  energy  can  be  arrived  at  by  extrapolating  empirical  formulae 
such  as  the  one  of  Reference  84  for  medium  impact  speeds.  It  is  estimated 
that  a windshield  of  laminated  vinyl-glass  required  to  defeat  the  impact  of 
a four  (4)  pound  bird  at  1.2  Mach  and  500  ft.  altitude,  has  to  be  2.75"  thick. 
At  thi3  material  thickness,  luminous  transmission  is  an  estimated  50%,  which 
is  far  below  the  acceptable  value. 

Flat  25"  x 19"  panels  of  thermally  strengthened  glass  with  rounded  edges 
and  clamped  edge  conditions  have  been  tested  in  connection  with  the  develop- 
ment o'  the  Concorde  windscreen  (Reference  54).  I .at  of  a four  pound  bird 
at  a speed  of  512  mph  is  the  operational  requirement.  Experimentally  observed 
penetration  speeds  were  considerably  higher  than  those  predicted  using  the 
formula  of  R.A.E.  Tech.  Mem.  106  (1952)  while  not  being  as  high  as  those 
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predicted  by  the  newer  formula  of  Note  373  of  1967.  In  the  experiments, 
the  edge  attachment  featured  a relatively  stiff  light  alloy  cast  mounting 
frame,  two  alternative  widths  of  packing  pieces  and  clamping  rings  (.84" 
and  .6  " encastration) , rubber  gaskets  and  no  glasst-'to-metal  contact.  No 
appreciable  differences  in  impact  behavior  between  the  two  encastration 
widths  were  observed.  The  forward  facing  abrasion  screen  was  reduced  in 
size  to  form  an  edge  rebate  to  permit  clamping  of  the  main  load  bearing  plies 
only.  Deflections  were  determined  by  measuring  the  distance,  between  the  glass 
face  and  the  displaced  plunger  rod.  Accuracy  of  the  measurements  was  pro- 
bably better  than  10%  measuring  deflections  of  .08"  to  .28".  Comparable 
deflections  for  polycarbonate  windshields  are  several  inches.  Data  for  a 
different  Concorde  test  panel  (Reference  54)  of  the  thermally  strengthened 
glass  "10-20"  are  as  follows:  three  (3)  main  plies  laminated  together  with 
PVB,  1.15"  total  thickness,  size  820  mm  x 520  mm  (32  1/4"  x 20  1/2"),  rounded 
corners,  tested  at  7°C  ambient  temperature,  53°  declined.  This  panel 
satisfied  the  operational  requirements . 
d,  EMPIRICAL  IMPACT  FORCE  EQUATIONS 

Impact  forces  due  to  birds  hitting  a non-~denting  plane  at  relative  normal 
velocities,  V,  are  tabulated  in  Reference  42,  The  minimum  possible  impact 
force, 

F = MV2/2s  (85) 

caused  by  the  change  of  velocity  to  zero  in  the  time  available,  was  obtained 
ay  assuming  constant  deceleration  rate.  As  a baseline,  this  force  was  deterr- 
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mined  for  a 2 pound  herring  gull  of  4"  effective  bird  dimension,  s.  The 
bird  dimension  was  scaled  for  larger  and  smaller  birds  as  the  cube  root  of 
bird  weight.  If  a sinusoidal  acceleration  were  assumed,  the  peak  acceleration 
would  be  11/2  times  mean  acceleration.  Forces  would  be  greater  if  the  bird 
carcass  did  not  spread,  as  it  would  result  in  the  acceleration  taking  place 
in  a shorter  distance  than  "s".  In  the  case  of  an  oblique  impact  at  an  angle 
of  incidence, oC , the  relative  normal  velocitv  is  V’coscrt  instead  of  V. 

A bird  acts  as  a semi-rigid  object  forming  a three-dimensional  mound 
(Reference  41).  This  mound  is  assumed  to  be  1/4  of  the  effective  bird  dimen- 
sion in  height.  With  no  structural  deformation,  its  acceleration  occurs 
during  the  linear  motion  of  the  center  of  gravity  of  the  bird  through  a 
distance  of  3/8  s.  The  formula  for  sinusoidal  peak  force  for  normal  impact 
(Reference  41) , 

F = .702  • W2/3  • V2  (86) 

(W  in  pound , V in  kt) 

was  experimentally  verified  in  connection  with  work  on  engine  intake  bird 
deflectors  from  readings  of  pressure  transducers  fitted  to  inclined  rigid 
grills  supporting  the  deflector  plate  and  measuring  both  the  vertical  and 
horizontal  forces  (Reference  72).  The  formula 

P = .685  • W2/3  • V2  • cos*  (87) 

was  found  to  be  valid  in  the  range  of  W - .1  thru  2.25  kg,  V = 105  thru 
320  kt  ando<=  45  thru  75°. 

e.  BIRD  IMPACT  ON  LEADING  EDGES,  FAN  BLADES 

Equations  86  and  87  were  also  applied  to  other  impacted  structures,  such 
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as  wing  leading  edges  (Reference  72) . The  response  of  composite  material 
cantilever  beams  to  soft  body  impact  representing  bird  impact  on  jet  engine 
fan  blades  was  studied  in  Reference  78.  Agreement  between  the  experimental 
results  and  one  degree  of  freedom  - classical  beam  theory  for  free-clamped, 
free-hinged  and  free-free  inodes  of  attachment  was  found  to  be  good.  Extension 
of  this  work  to  a two-degree  of  freedom  lumped  parameter  model  has  provided 
excellent  predictions  of  measured  torsional  vibrations  induced  by  impact.  The 
second  part  of  these  studies  deals  with  the  transmission  of  planar  forces, 
created  by  a pressure  profile,  from  low  modulus  transparencies,  such  as  homo- 
geneous, laminated  and  fiber  reinforced  materials,  into  different  types  of 
adjoining  fuselage  structure. 

2.  WINDSHIELD  EXPERIMENTS  OTHER  THAN  BIRD  STRIKE 

An  analysis  method  is  mentioned  (Reference  79)  that  predicts  the  temp- 
eratures of  and  heat  transfer  within  window  materials  by  both  conduction  and 
radiation.  The  heat  balance  equations  in  the  computer  routine  account  for 
emission,  attenuation,  and  absorption  cf  radiant  energy  within  the  glaze. 

The  heat  transfer  test  results  provide  experimental  data  on  realis  tic  window 
components.  This  effort  could  be  considered  a prerequisite  for  conducting 
thermal  stress  and  failure  analysis  of  window  systems  in  thermal  flight  en- 
vironments . 

The  thermal  shock  profile  of  a windshield  is  the  temperature  variation 
thru  the  transparency  assembly  at  the  specific  flight  profile  that  produces 
the  maximum  temperature  differential.  The  thermal  shock  temperature-time 
history  of  the  original  glass  windshield  of  the  F-lll,  repeated  for  75  cycles 
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of  1600  sec  each,  is  typical  for  high  performance  aircraft  windshields.  This 
windshield  was  composed  of  .11”  Herculite  II,  .08"  interlayer,  .11"  Herculite 
II  (Reference  80) . The  inner  and  outer  initial  test  temperatures  were 
T]_  = -35°F  and  T2  = -18°F,  respectively.  The  temperature  differential  was 
high  during  half  of  the  flight  cycle,  Ti  peaking  at  350°F  and  T2  at  275°F. 

The  development  of  one  windshield  required  extensive  structural  testing 
of  the  materials  as  well  as  the  interlaminar  bonds  (Reference  9) . The  tests 
were  oriented  toward  pressurization  loading  in  one,  and  separate  impact 
performance  tests  in  another  program.  The  pressurization  loading  was  simula- 
ted in  a test  rig  using  a composite  windshield  of  1.7"  total  thickness 
and  subjecting  it  to  over  1,000  hrs  of  cyclic  loading  at  a range  of  temper- 
atures from  -65 °F  to  70°F.  Figure  40  shows  the  pressure  cycle  used  in  the 
investigation.  The  windshield  survived  without  damage  and  greatly 
increased  the  understanding  of  the  structural  capabilities  of  this  type 
of  windshield. 

3.  OTHER  IMPACT  ANALYSES 

The  following  research  investigations  were  found  to  be  valuable  back-up 
information. 

a.  TRANSIENT  BEAM  ANALYSIS 

Transient  analysis  was  conducted  to  investigate  the  impact  behavior  of  a 
free-free  beam  under  an  impulsive  force  applied  at  the  beam  center  in  Ref- 
erence 81.  Non-dimensional  deflections,  wpl/PL^,  were  computed  in  terms  of 
non-dimensional  beam  location  (actual  locatTon/tip  location)  for  different 
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non-dimensional  instances  of  time, 

T = t/tj  with: 

tj  - 7T/f0o  (88) 

(Oo-lST 

V ML**  (first  mode  of  free  vibration) 

where  Ma  mass  of  beam  per  unit  length  (pound-sec^ /in^) , and  P/L=  inertial  load 
on  the  beam  (pound/in).  Figure  41,1  shows  the  initial  deflectional  response 
up  toT=.25.  A "whip  lash"  motion  of  the  end  of  the  beam  may  be  noted  in 
the  spanwise  wave  travel.  The  time  histories  of  center  deflection  and  tip 
deflection  relative  to  the  center  of  mass  (the  reference)  are  presented  in 
Figure  41.2. 

b,  ELASTIC  PLASTIC  ANALYSIS  OF  PLANE  STRUCTURES 

The  following  is  a description  of  a very  elaborate  elastic-plastic  impact 
analysis  (Reference  48) : 

(1)  The  elastic-plastic  response  of  certain  plane  structures  under- 
going large  deformations  due  to  dynamic  loading  has  been  investigated. 

Various  methods  were  used  in  loading  the  specimens.  In  some  tests  heavy  lead 
rings  were  dropped  on  flat  and  pointed  anvils.  On  other  tests,  stationary 
lead,  copper  and  aluminum  rings  carrying  weights  around  the  periphery  and 
those  with  no  additional  attached  masses  were  subjected  to  contact  explosives 
directed  radially  inward  at  a single  point.  Aluminum  cantilevers,  rectangu- 
lar and  triangular  frames  and  semicircular  arches  were  loaded  by  magnetomotive 
and  explosive  impulses. 

(2)  The  aluminum  rings  with  5"  to  6.5"  outer  diameter,  1/8"  to  1/4" 
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[FROM  REFERENCE  81] 


Figure  41.  Beam  Response  Due  to  Impulsive  Load 
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radial  thickness  and  up  to  3/8"  wide,  were  machined  from  tubing  and  annealed 
at  350 °C  for  15  min.  The  aluminum  had  the  following  mechanical  properties: 
14.9  ksi  yield  stress,  8 to  10  x 10^  psi  elastic  modulus  and  .2  to  .82  x 10^ 
psi  plastic  modulus.  Cantilevers  and  frames  of  aluminum  were  made  from 
.032"  to  .036"  thick  aluminum  sheet  material,  10"  long  and  3"  wide.  The 
masses  attached  to  the  plane  structures  were  1/2  pound  steel  blocks, 

(3)  The  sharp  indenter  initiated  the  plastic  deformation  process  at 
a free  falling  speed  of  14  ft/sec.  Point  explosive  loading  by  No,  6 I.C.I. 
nylon  detonators  caused  1.5xl0-2  pound-sec  impulses  with  detonative  defor- 
mation times  in  the  order  of  ,01  sec.  0.37"  diameter  and  3/4"  long  mild 
steel  bullets  moving  at  695  ft/sec  were  used  in  the  high  speed  bullet  impact 
tests.  The  distributed  magnetomotive  loading  method  imposed  a prescribed 
uniform  radial  velocity  on  the  ring.  Magnetomotive  as  well  as  detonative 
loads  delivered  to  the  frames  were  calculated  from  the  swing  of  1.2  pound 
ballistic  pendulum  rigidly  fixed  to  the  frame.  Initial  structural  velocities 
were  determined  to  be  7800  to  9000  in/sec,  requiring  high  speed  photographing 
at  15,000  frames  per  second. 

(4)  The  impact  behavior  of  the  ring  subjected  to  impulsive  point 

loads  was  simulated  by  large  deflection-analysis  assuming  elastic-linear 
strain  hardening  or  rigid-plastic  material  behavior  and  the  formation  of  four 
plastic  hinges:  two  stationary  hinges  at  the  point  of  impact  and  a point 

diametrically  opposite,  and  two  moving  hinges.  Rigid-plastic  analyses  were 
found  to  be  valid  only  in  cases  for  which  the  initial  kinetic  energy  was 

an  order  of  magnitude  greater  than  the  elastic  strain  energy  of  the  struc- 
ture at  the  onset  of  plastic  deformation. 
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(5)  Two  analysis  methods  were  used  to  investigate  the  loading  of 
rings,  and  particularly  arches,  cantilevers  and  frame,  One  method  assumed 
a simple  approximate  model  of  the  ring,  consisting  of  concentrated  masses 
around  the  circumference  connected  to  each  other  by  weightless  links  possess- 
ing the  strength  properties  of  the  material,  The  dynamic  equation  for  each 
mass  of  the  ring  were  solved.  The  second  method  utilizes  Lagrange's  energy 
equations.  Both  methods  were  in  close  agreement. 

(6)  Insignificant  difference  between  the  calculated  and  experiment- 
ally determined  profile  of  the  deformed  structure  were  observed  as  long  an 
the  ratio  of  structural  length  to  the  link  length  remained  more  than  ten: 
As/s >10.  Prediction  improved  with  Increasing  number  of  nodal  points.  Time 
interval, At,  cannot  be  chosen  arbitrarily  and  depends  upon  the  wave  speed 
thru  the  structural  link  length, As.  It  was  found  that  divergence  occurs 
when 

(At/As)(E/j)1/2>l  (89) 

Another  analytical  effort  (Reference  82)  demonstrated  the  injury  reduc- 
tion capability  of  energy  absorbers.  Existing  statistical  data  on  crash 
accelerations  and  man-weight  distributions  were  used  as  inputs  to  a MIMIC 
computer  program.  Dynamic  response  indices  and  stroke  lengths  were  computed 
for  each  acceleration-man-weight  combination.  This  program  can  be  used  as 
a guideline  in  selecting  candidate  energy  absorbers. 

c.  EXPLOSIVELY  LOADED  ALUMINUM  SHELLS 

The  computer  code  PETR0S  (Plastic-Elastic  Transient  Response  of  Cylindri- 
cal Shells  with  Large  Deflections)  used  finite-difference  method  to  calculate 
dynamically  induced  deformations  of  general,  thin  shells  (Reference  83). 
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PETR0S  handles  symmetrical  loading  conditions  only.  No  structural  reinforcing 
supports' were  taken  into  consideration.  Classical  boundary  conditions  were 
employed.  This  program  was  reviewed  to  determine  its  ability  to  simulate 
the  actual  F-lll  windshield  bird  strikes  test  results.  The  boundary  conditions 
that  can  be  handled  and  symmetry  restrictions  make  the  program  inappropriate 
for  simulating  the.  F-lll  windshield  test  results. 


d.  HAILSTONE  ANALYSIS 

The  "DEPR0SS  3"  program  may  be  applicable  to  the  bird  strike  impact  pro- 
blem once  the  failure  criteria  for  transparencies  material  are  clarified 
(Reference  75).  This  program  analyzes  the  dynamic  elastic-plastic  axisym- 
metric  responses  of  Impulsively  loaded  plates  and  shells  of  revolution. 

Thomson  and  Hayduk  adapted  "DEPR0SS  3"  with  minor  modifications  in  the  ini- 
tial velocity  input  suction  and  the  mass  ratio  of  hailstone/plate,  improving 
their  previous  work  by  including  the  effects  of  combined  bending  and  membrane 
action.  They  succeeded  in  predicting  the  size  of  dents  caused  by  hail  impact. 


e.  OTHER  COMPUTER  PROGRAMS 

The  following  programs  can  be  used  to  simulate  experimental  bird  strike 
test  results.  SABOR  5/DRASTIC  V is  applicable  only  for  impact  analysis  of 
shells  of  revolution  without  stiffening  members.  MAGIC  and  NASTRAN  have  the 
required  stiffened  shell  segment  capability.  In  MAGIC,  however,  the 
transient  force  (or  velocity)  input  response  capability  is  lacking^  while 
NASTRAN  possesses  it. 
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APPENDIX  III 

DESCRIPTION  OF  FINITE  DIFFERENCE  COMPUTER  ROUTINES 

User  oriented  descriptions  of  the  four  FORTRAN  program  versions 
developed:  FLAT  PLATE,  PLATE  WT-18,  CONICAL  SHELL  SEGMENT,  and  F-lll  W/S 

(FM-1)  are  presented  in  this  appendix.  These  computer  programs  are  based  on 
the  finite  difference  methods  of  Section  III  and  were  applied  in  the  analy- 
ses of  Section  IV.  Listings  of  the  programs  are  provided  in  Appendix  IV. 

For  familiarizing  the  user  with  the  programs,  the  appropriate  card  listings 
accounted  for  by  numbering  from  the  start  of  the  individual  program  may  be 
looked  t'p  as  soon  as  they  are  referred  to  in  the  following  guidelines. 

All  functions  of  executing  the  computations,  such  as  assigning  the  time 
varying  variables,  zeroing  variables  for  initial  conditions,  controlling 
variables,  as  of  stiffness  and  support,  and  finally  solving  the  structural 
finite  difference  equations,  are  contained  in  the  four  different,  separate 
computer  routines  "MAIN".  The  only  type  of  subroutine  that  has  been  es- 
tablished, provides  graphical  support  to  process  output  data  by  means  of 
the  CALCOMP  plotter.  This  subroutine,  "GRAPH’,  is  listed  only  in  connection 
with  the  F-lll  W/S  (FM-1)  program.  Similar  subroutines  may  be  devised  and 
attached  to  the  remaining  programs. 

1.  FLAT  PLATE  PROGRAM 

a.  Introduction 

Based  on  a mathematical  model  for  impulse  loading,  a computer  pro- 
gram for  calculation  of  time  histories  using  finite  difference  equations  of 
the  types  of  Equations  50  has  been  developed  in  the  Analog  and  Hybrid 
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Computation  Branch,  Computer  Center,  WPAFB,  Ohio.  This  program  (Reference 
43)  was  written  for  a flat  plate  of  6 x 6 elements  to  consider  flexural 
bending  of  the  plate,  rotational  clamping  stiffness  of  the  attachment  frame, 
and  viscous  damping  and  spring-back  behavior  of  the  mounting  system.  This 
required  a 7 x 7 mesh  that  also  needed  half-way  points.  Variables  were 
dimensioned  in  13  x 13  arrays.  "D0"  loops  were  used  to  define  problem 
variables  within  the  mesh  and  on  the  boundary.  Numerical  integration  util- 
ized the  point  slope  formula  (Euler's  Method)  which  proved  to  be  adequate 
for  time  steps  smaller  than  .05  msec. 

The  sequence  of  operating  in  a marching  procedure  on  the  representative 
straight  beam  equations  in  cross-beam  arrangement  to  calculate  the  instan- 
taneous variables  is  shown  in  Figure  42.  The  "x"  and  "y"  markings  in  the 
patterns  represent  the  generation  of  numerical  values  of  x-  and  y-  variables, 
respectively.  Initial  conditions  consider  all  variables  (internal  shear 
loads  (VX,  VY) , bending  moments  (MX,  MY),  slopes  of  deflection  (OX,  OY), 
and  displacements  (W))  zero.  Figure  43  reflects  the  optional  distribution 
of  constants  in  the  finite  difference  scheme:  distribution  of  bending  stiff- 

nesses (EIX,  EIY,  EIF),  mounting  properties  (KW,  KWD,  KOX,  KOY,  KOXD,  KOYD) 
along  the  edges,  time  dependent  pyramidal  impact  loads  (Q1>Q2>Q3)  and  mass 
properties  (DMM,  DMF) . 

For  reason  of  simplicity,  the  basic  flat  plate  program  is  presented  in 
its  shortest  form  (Appendix  IV-1).  Bending  stiffnesses  are:  EIX-EIY^EIF^EI . 

Mounting  properties  along  the  edges,  KW,KWD  and  K0X=K0Y=K0,  are  constant. 
Rotational  dampening  is  avoided:  KOXD=KOYD»0.  The  components  of  load  can 

be  specified  at  any  grid  point.  Distributed  loads  must  be  treated  as  con- 
centrated at  the  nodes.  The  pyramidal  load  distribution  (Figure  43.3)  is 


187 


AFFDL-TR-74-155 


TYPICAL  REPRESENTATIVE  EQUATION 


GENERATED  PATTERN 

(FOR  GRID  SIZE  NX*NY=13*13) 


0X( I , J ) »(W( 1+ 1 , J)-W(I^ 1 , J ) )/DS 
(AND  FOR  0Y,  INSIDE  AND  BOUNDARY) 
(LISTED  AT  LINES  67  AND  68  OF  FLAT 
PLATE  COMPUTER  ROUTINE  USED) 


It  II 


1,1  r?TxTxT-x-,-  x-i  1,15 

y+y+y-f-y-hy-j-y-j-y 

f~  x x -f—  x -4—  x — j—  x-4-  x -j 

y + y y-+- y-f- y y -f- y 

f-j-  x-+-  x-j 

y+y+y-ry-f-y-j-y-H-y 

13,!  l-x-l-x_X-x_J_x_J_x_J_x_li5ti5 


0X ( 1 , J ) =0X( 2 . J )/( 1 . +DS • K0X/( 2 . • EIC ) ) 
(AND  FOR  0Y , BOUNDARY) 

(LINES  72  THRU  75) 
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13,13 


MX(  II , J )=-K0X*0X( 1 1, J) • SIGN  "2" 

(AND  FOR  MY,  BOUNDARY)  (LINES  80  AND  82) 

VX(IT , J)  =(KW(  I)  • W(  II , J ) + KY/JD ( IK)  "2K 

•WD(II,J))«SIGK 

(AND  FOR  VY,  BOUNDARY)  (LINES  92  AND  102) 


Figure  42.  Sequence  in  Cyclic  Loop  for  Calculating 
Finite  Differences 
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MX(1I,J)=(0X(II+1  , j)-0X(  1 1-1  , J ) ) • IilC/DS 
(AND  FOR  MY,  INSIDE  AND  BOUNDARY) 
(LINES  109  AND  112) 


1,1  1,1'} 


J-  J— If V_i_v !_i_] i_J  .1 


13,i 


LJ 


- l-X-l-X—i. 


13,13 


VX(II,J)=(MX(II+1 ,J)-MX(II-1 ,J))/DS  "1" 

(AND  FOR  VY,  INSIDE  AND  BOUNDARY) 

(LINES  117  AND  118) 


VX0( I)=2.  • VX( 1 , 1 )-VX(2 , 1) 

(AND  FOR  VYO,  ETC.,  OUTSIDE  AND  BOUNDARY) 
(LINES  122  THRU  125) 


WDD(I,J)=(Q(I,J)  -VXP+ VXK-VYP+- V YM  )/DM 
(AND  FOR  WD  AND  W) 

(LINES  141,  142  AND  144) 


1" t "i _r  t"  f “i-  f”i 

- f -i- -i-  — ! — i — l j J— j 
lltiJL.t  * iiifi 
P— 2 z z~.— z — !“Z — — z 


f 

i 

[_ J-Lfj  f L 1 j 

13,1  ^ Z 2 -i 


Figure  42.  (Concl.) 
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presented  as  an  example  of  external  transient  load  application  only.  In 
the  basic  flat  plate  program,  loads  of  .1  pound  as  "quasi-zero"  loads  remain 
constant  at  all  interior  nodal  points  and  along  the  edges  (Q1=Q2=Q3=Q) . 

All  nodal  masses  are  constant:  DMM=DMF=DM. 

b.  Basic  Program  Listing 

Program  FLAT  PLATE  has  been  listed  in  two  versions  (Appendix  IV-1)  to 
demonstrate  ways  of  simple  modification  of  the  same  basic  routine  for  versa- 
tile utilization.  The  unmodified  listing  is  described  first,  excluding  the 
three  modifications  each  in  the  options  "ONE"  and  "TWO"  by  computer  card 
additions.  For  this  purpose,  the  comment  card  statements  are  repeated  with 
additional  comments. 

Following  the  title,  FLAT  PLATE  PROGRAM  FOR  FIXED  GRID  SIZE  (13  x 13), 
PARAMETRICAL  RUN  *****  ONE  *****  (or  PARAMETRICAL  RUN  *****  TWO 
*****),  the  second  comment  is:  ASSIGNING  THE  TIME  VARYING  VARIABLES. 

This  is  followed  by  dimensioning  the  variables  W,  WD,  WDD,  OX,  OY,  VX,  VY, 

Q,  MX  and  MY  and  the  single  array  variables  VXO,  VX14,  VYO  and  VY14,  To 
change  the  applicable  problem  size  of  the  program  to  a grid  size  other  than 
13  x 13,  the  dimensions  of  these  parameters  as  well  as  dimensional  notations 
in  the  equations  of  the  routine  must  be  changed. 

The  comment  INPUT  DATA  REQUIRED  (NUMERICAL  EXAMPLE  IS  ARBITRARY)  is 
followed  by  six  necessary  plate  constants.  The  geometry  of  the  structure 
is  determined  by  the  grid  parameter,  DS,  that  reflects  the  side  of  the 
square  panel  in  inches.  The  bending  stiffness,  El  (pound-inch^) , of  the 
pane  and  the  frame  boundary  depends  upon  Youngs  modulus  (psi) , grid  size 
and  pane  thickness  (inch).  DM  (pound-sec^/in)  is  the  lumped,  constant  mass 
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at  the  grid  point  location.  Boundary  conditions  are  presented  at  the 
nodal  edges:  spring  constant,  KW  (pound/in),  rotational  spring,  KO  (pound- 

in/rd),  and  viscous  damping,  KWD  (pound-sec/in).  In  addition,  there  are 
three  time  variables  controlling  the  program.  As  a numerical  example; 

DT  = .00001  (time  step,  sec)  and  NEND  = 1000  (number  of  last  time  step) 
yield  plate  vibration  of  1/100  sec  real  time.  NP  determines  the  print 
sequence,  to  suppress  print-out  at  time  steps  not  desired. 

To  SET  INITIAL  CONDITIONS,  all  distributed  plate  variables  (13  x 13) 
are  zeroed  at  time  step,  NNNN  ■ 1,  using  a double  DO-loop,  These  are: 


0X(I,J),  0Y(I,J) 

Slopes 

(rd) 

KX(I,J),  MY(I,J) 

Bending  moments 

(pound-in) 

VX(I,J),  VY (I, J) 

Shear  loads 

(pound) 

WDD(I,J) 

Normal  accelerations 

(in/sec^) 

WD(I,J) 

Velocities  of  displacement 

(in/sec) 

W(I,J) 

Deflections 

(in) 

The  external  loads  Q(I,J)=.l  are  "quasi-zero"  loads  of  .1  pound  (or  another 
small  value)  each  and  are  assigned  to  stabilize  the  equations  when  the 
variables  above  eventually  turn  to  very  small  values.  Much  higher  loads, 
Q(I,J),  replace  this  .1  pound  loading  in  case  of  an  impact  load  analysis 
(such  as  in  run  "TWO").  For  another  type  of  transient  analysis,  the  DO-loop 
may  be  used  to  initially  set  WD(I,J)  non-zero  (run  "ONE"). 

The  MAIN  D0-L00P  OF  COMPUTATION  PER  TIME  STEP  NNNN  is  to  DEFINE  0-,  M- 
AND  V-VARIABLES  IN  13  x 13  ARRAYS.  This  is  the  executing  part  of  the  routine 
with  multi-nested  DO-loops  recalculating  OX(I,J),  OY(I,J),  MX(I,J)  etc. 


» 


A 
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In  the  sequence  of  determining  the  finite  differences,  the  typical  repre- 
sentative equations  with  line  reference  numbers  (line  numbers  of  run  "ONE") 
are  repeated  from  Figure  42:  OX,  OY  INSIDE/BOUNDARY,  (a)  OX(I,J)= 

(W(I+1) ,J)-W(I-1,J))/DS  (and  for  OY,  inside  and  boundary)  (listed  at  lines 
66  and  67),  (b)  0X(1,J)  = OX(2,J)/l.+DS-KO/(2. *EI))  (and  for  OY,  boundary) 

(lines  68  thru  71);  MX,  MY,  VX,  VY  ON  BOUNDARY,  (c)  MX(II,J)=-KO-OX(II,J) • 
SIGN  (and  for  MY,  boundary)  (lines  78  and  79),  (d)  VX(II,J)=(KW(I)-W(II,J)+ 
KWD*WD(II,J))'SIGN  (and  for  VY,  boundary)  (lines  80  and  81);  MX,  MY  INSIDE/ 
BOUNDARY,  (e)  MX(II ,J)=(0X(II+1,J)-0X(II-1,J)) ‘EI/DS  (and  for  MY,  inside 
and  boundary)  (lines  86  and  87);  VX,  VY  INS IDE /BOUNDARY,  (f)  VX(II,J)= 
(MX(II+1,J)-MX(II-1, J))/DS  (and  for  VY,  inside  and  boundary)  (lines  92  and  93); 
VX,VY  OUTSIDE  BOUNDARY,  (g)  VXO (I)“2. *VX(1,I) -VX (2 , I)  (and  for  VYO,  etc., 
outside  and  boundary)  (lines  99  thru  102) ; CALCULATE  WDD  AND  INTEGRATE, 

(h)  WDD ( I , J) = (Q ( I , J ) -VXP+VXM- VYP+VYM ) / DM  (and  for  WD  and  W)  (lines  114, 

115  and  116).  Finally,  the  IDENTIFICATION  OF  VARIABLES  FOR  PRINT-OUT  and 
WRITING  OUTPUT  VARIABLES  yields  the  required  print-out  at  any  time  step 
NNNN.  This  is  the  END  OF  MAIN  D0-L00P. 

c.  Different  Versions  "ONE"  and  "TWO" 

Example  "ONE"  is  the  first  listing  in  Appendix  IV-1.  Comment  cards 
inserted  at  three  locations : BEGIN  OF  RUN  "ONE"  MODIFICATION  and  END  OF  RUN 

"ONE"  MODIFICATION  indicate  the  changes  made  to  the  basic  program.  Following 
SET  INITIAL  CONDITIONS  (line  32),  all  time  varying  variables  are  zeroed, 
except  for  the  quasi-zero  load  Q(I,J)=,01  (pound)  and  the  initial  velocity 
WD(I,J)=100  (inch/sec).  The  main  D0-loop  of  computation  per  time  step, 

NNNN  (lines  56  thru  157) , is  activated  two  times  (lines  25  thru  159) . 
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The  first  activation  generates  all  time  varying  variables  of  five  consecutive 
steps  and  print  all  these  variables;  the  second  one  calculates  50  steps, 
printing  every  10th  step  only.  The  reactivation  of  the  DO-loop  a second 
time  is  not  necessarily  required  to  get  the  results.  This  was  done  for  ease 
of  programming.  It  can  be  done  with  one  use  of  the  DO-loop. 

For  structural  reasons  as  well  as  loading  symmetry  and  symmetry  of 
initial  conditions,  only  parts  of  the  generated  patterns  of  numerical  values 
are  printed-out  in  such  a way  that  no  value  in  a pattern  is  repeated  unne- 
cessarily. This  print-out  is  on  the  first  and  second  page  following  the 
listing  "ONE".  It  can  be  seen  that  zero  values  half  way  between  the 
nodal  points  of  grid  generation  are  omitted.  The  WRITE  and  FORMAT  state- 
ments under  WRITING  OUTPUT  VARIABLES  reflect  the  tabulation  per  page.  Use 
of  this  basic  13  x 13  plate  program  for  nonsymmetrical  cases  suggests  that 
these  statements  are  modified  such  that  the  entire  pattern  of  numbers  is 
printed. 

The  outer  DO-loop,  DO  3 NN=1,  3 (line  24),  conducts  two  additional  runs 
resulting  in  four  more  pages  of  print-out.  The  print-out  on  the  third  and 
fourth  page  has  the  same  velocity,  WD(I,J)*100,  initially  applied.  The  quasi- 
zero load,  Q(I,J),  was  set  at  .01  pound  instead  of  ,1  pound  to  show  the  dif- 
ference. The  generated  numbers  on  first  and  third  (and  second  and  fourth) 
pages  should  be  compared  for  similarity  $f  results.  The  third  and  last  run 
with  reversed  input  signs,  WD(I,J)**-100  and  Q(I,J)*-.01,  yields  the  output 
with  the  same  numbers  as  the  second  run  but  with  reversed  signs. 

Example  "TWO"  is  the  second  listing  in  Appendix  IV-1.  Comment  cards 
inserted  at  three  locations  show  where  the  mbdifications  to  the  basic  program 
were  made.  The  initial  condition  is  a linearly  increasing  load  at  the  center 
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of  the  plate,  Q(7,7),  rather  than  an  initial  non-zero  velocity  condition. 

The  linear  load  increase  is  changed  to  a constant  load  level  at  .1  sec. 
Computation  is  terminated  at  .25  sec. 

Two  runs  where  conducted  for  comparison  of  results,  using  the  loop,  DO 
» NN=1,  2 (line  26).  The  first  run  required  1250  steps  to  yield  .25  sec 
of  time  event;  the  second  required  20  times  as  many.  Displacements  can  be 
seen  to  increase  almost  linearly  thru  .1  sec  and  stay  fairly  constant  from 
there  thru  .25  sec  in  both  cases.  The  accelerations,  WDD,  reflect  the 
highest  deviations  when  comparing  corresponding  numbers  of  the  two  runs. 

2.  PLATE  PROGRAM  FOR  CASE  WT-18 

This  is  an  expansion  of  the  basic  flat  plate  program  to  simulate  the 
deflection  data  of  impact  case  WT-18.  The  comment  cards  of  the  listing 
(Appendix  1V-2)  are  repeated  with  an  explanation  of  the  differences  between 
this  and  the  basic  plate  program.  The  self-explanatory  title  is  FLAT  PLATE 
PROGRAM  FOR  VARIABLE  GRID  SIZE  (UP  TO  21  x 21)  SET  UP  FOR  CASE  WT-18. 

All  dimensioned  variables  are  sized  such  that  21  is  the  maximum  allowable 
grid  size  (ASSIGNING  THE  TIME-VARYING  VARIABLES) . The  expanded  variables 
as  compared  with  those  of  the  basic  program  are:  displacement,  W,  is 
replaced  by  W,  WMDT  (deflection  at  previous  time  step);  velocity,  WD,  and 
acceleration  WDD  remain  the  same;  OX,  OXMDT  (slope  at  previous  time  step), 

OXD  (rate  of  slope)  for  OX;  OY,  OYMDT,  OYD  for  OY;  VX,  VY,  Q,  VXO,  VYO  remain 
the  same;  VXP1  and  VYP1  for  VX14  and  VY14,  respectively;  MX,  MXEL  (elastic 
component  of  bending  moment) , MXVC  (viscous  component  of  bending  moment)  for 
MX;  MY,  MYEL,  MYVC  for  MY.  The  constant  boundary  conditions  are  replaced  by 
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these  distributed  ones:  KW,  KWD,  KOX,  KOY,  KOXD  and  KOYD. 

These  are  the  pre-calculated  parameters  for  grid  sizes  (VARIABLES  FOR 
DIFFERENT  GRID  SIZES):  L(l)-5,  L(2)=9,  L(3)=13,  L(4)=17,  L(5)=21,  L(7)=29 

and  L(8)=33.  One  set  of  the  variables,  L,  A and  B may  be  chosen  for  grid 
size  determination.  A and  B depend  on  grid  size  L. 

Input  data  (INPUT  DATA  REQUIRED  (CASE  WT-18))  are  expanded  to  side  lengths 
X and  v instead  of  DS  of  the  basic  plate  program.  Any  choice  of  Nl=l  thru  8 
selects  the  parameters  E,  A and  B.  The  parameters  EMOD,  THCK,  EIXF  and  EIYF 
stand  for  El,  and  MPANE  and  MFRAME  instead  of  DM.  The  basic  KW,  KWD  and  KO 
are  replaced  by  KWDN  (KW  in  downward  mode):,  KWUP  (upward) , KWDC,  KOXC,  KOYC 
(all  for  local  distribution) . Added  parameters  KOXDC  and  KOYDC  are  rotational 
edge  dampening;  directional  viscous  in-parie  and  in^frame  properties  are 
KVSX,  KVSY,  KVSXF  and  KVSYF . QTOT  is  the  total  impact  load  that  is  pyramidally 
distributed.  TQ  is  the  impact  time.  DT,  NEND  and  NP  remain  as  in  the  basic 
program. 

The  reader  is  referred  to  lines  55  thru  73  for  CALCULATION  OF  DEPENDENT 
VARIABLES.  These  are  dependent  grid  (NG  thru  DY) , stiffness  (EIX  and  EIY) , 
mass  (DMM  and  DMF)  and  impact  endurance  constants  (TQP,  TQM)  to  set  up  the 
finite  difference  model.  The  DISTRIBUTION  OF  INITIAL  CONDITIONS  AND  BOUNDARY 
CONSTANTS  are  lines  75  to  131,  where  the  main  DO-loop  of  computation  begins . 

All  (I,J)  variables  are  zeroed  for  the  first  step  NNNN=1  (lines  81  thru 
98);.  Uniform  distribution  of  boundary  constants  KW,  KWD,  KOX,  KOY,  KOXD, 

KOYD  is  accomplished  (DOWNWARD  AND  UPWARD  MODE  OF  KW).  For  SPECIAL  CLAMPING 
FOR  CASE  WT-18,  some  KW  values  In  the  uniformly  distributed  pattern  are 
selectively  superposed  by  new  values.  The  operations  under  DISTRIBUTE: 
STABILIZING  PERMANENT  LOAD  (lines  121,  122,  123):,  and  SUPERPOSE  PYRAMIDAL 
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IMPACT  LOAD  WT-18  (lines  127  thru  131)  are  self-explanatory. 

All  the  following  comments  concern  operations  within  the  main  DO-loop 
starting  at  line  135  (MAIN  DO-LOOP  OF  COMPUTATION  PER  TIME  STEP  NNNN) . 

IMPACT  TIMING  WT-18  (line  140)  replaces  the  continuous  quasi-zero  load  with 
the  time  dependent  impact  load  distribution  and  reestablishes  the  continuous 

i 

load  after  impact  (line  36).  The  remaining  comments  are  identical  with  those  j 

of  the  basic  plate  program:  DEFINE  0-,  M-  AND  V-VALUES  IN  NG  x NG  ARRAYS;  | 

j 

OX,  OY  INSIDE/BOUNDARY;  ON  BOUNDARY;  MX,  MY,  VX,  VY  ON  BOUNDARY;  MX,  MY  INSIDE/  j 

BOUNDARY;  VX,  VY  INS IDE/ BOUNDARY;  VX,  VY  OUTSIDE/BOUNDARY ; CALCULATE  WDD  AND  I 

INTEGRATE;  IDENTIFICATION  OF  VARIABLES  FOR  PRINT-OUT;  WRITING  OF  OUTPUT  ! 

DEFLECTIONS  FOR  WT-18  and  END  OF  MAIN  DO-LOOE.  t 

i 

The  reader  may  compare  the  computational  expansions  of  this  program  with 
the  basic  one.  The  sample  of  results  printed  out  following  the  listing  repre-  J 

t 

sents  Figure  32  ana  is  self-explanatory.  ' 

3.  CONICAL  SHELL  SEGMENT  PROGRAM 
a.  Introduction 

The  concept  of  finite  difference  techniques  in  the  computer  routine 

"flat  plate"  (Appendix  IV-1)  was  applied  in  the  development  of  the  program  i 

t 

for  dynamic  stress  analysis  of  conical  shell  type  windscreens  and  canopies.  f 

Such  a program,  which  later  could  be  used  for  typical  systems  hardware  ; 

design,  was  not  readily  available  at  the  time  development.  This  monolithic 

conical  shell  routine  of  equal  grid  spacing  and  consistent  properties 

(Appendix  IV-3)  was  developed  as  a prerequisite  for  the  F-lll  windshield 

routine  (Appendix  IV-4) . j 

The  user  should  have  no  difficulties  applying  this  basic  code  and  gaining  j 

! 
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further  understanding  of  it  by  studying  the  F-lll  applied  code  that  is 
described  in  detail  (Appendix  III-4) . 

Serving  as  a proof  for  operational  readiness  of  this  dynamic  shell  pro- 
gram, symmetrical  17  x 17  stress  and  deflection  matrices  were  observed  as 
a result  of  a symmetrical  impact  loading  history.  The  current  listing  is 
set  up  with  initial  velocity  WD(I,J)*50  (in/sec)  and  provides  symmetrical 
results.  The  example  print-out  at  the  end  of  the  listing  includes  variables 
in  reduced  17  x 17  matrix  tabulation  at  time  steps  NNNN=1  and  4000,  equiva- 
lent to  .0001  and  .4  seconds  (At*. 0001  sec)  after  the  start  of  time  history. 
The  attached  two  by  three,  self-explanatory  print-out  tables  were  generated 
by  the  WRITE  statements  of  lines  197  thru  230  (Appendix  IV-3) . 
b.  Discussion  of  Program  Listing 

Additional  comments  concerning  the  modifications  necessary  for  a switch 
from  the  previously  described  flat  plate  geometry  are  provided.  The  reader 
is  referred  to  Appendix  III-l  for  the  basic  comments  still  valid  here.  The 
title  is  self-explanatory : DYNAMIC  CONICAL  SHELL  PROGRAM  (17  x 17) , ANY 

LOADING  HISTORY,  FIXED  TANGENTIAL  SUPPORT,  X-VARIABLES : MX,  UN,  FT  AND  OX, 
FN,  UT  (6  FIRST  ORDER  DIFFERENCE  EQUATIONS),  Y-DIRECTION:  MY,  UN  and  OY,  V 
(4  EQUATIONS  COUPLED  WITH  X-VARIABLES) . 

ASSIGNING  THE  TIME-VARYING  VARIABLES  (line  9)  incorporates  the  changes. 
The  previously  used  plate  variables  W,  WD,  WDD  are  renamed  by  displacements, 
velocities  and  accelerations  in  normal  directions:  UN,  UND,  UNDD.  The 
tangential  displacement,  UT,  is  newly  introduced.  Slopes  OX  and  0Y  remain. 
Internal  load,  V,  is  replaced  by  FN,  FT,  V.  External  loads  are  QN  and  QT 
(normal  and  tangential)  instead  of  Q.  Masses  DMM  are  calculated  instead  of 
using  a constant  distributed  mass,  DM.  Additions  are  one-dimensional 
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variables:  FNO,  FN1,  VO,  VI,  R (radii),  DX,  DS,  DTHX  and  THX  (radial  constants), 

and  EIY.  MX  and  MY  remain  the  same. 

INPUT  DATA  REQUIRED  (NUMERICAL  EXAMPLE  IS  ARBITRARY)  are  geometrical 
dimensions  SIDE,  RMAX,  UMIN,  THETA  that  replace  DS  of  the  flat  plate  program. 

A switch  was  made  from  parameter  El  to  THICK  and  E (thickness,  Young’s 
modulus),  and  specific  density,  SPDEN,  instead  of  DM.  Two-directional  edge 
stiffnesses  (KW)  are  KUN  and  KY;  rates  are:  KUND  and  KYD.  Rotational  springs 
KOX  and  KOY  stand  for  KO  of  the  basic  plate  program  (Appendix  IV-1).  DT, 

NEND  and  NP  stay  identical.  The  CALCULATION  OF  THE  REMAINING  INPUT  is 
FINITE  DIFFERENCES  OF  CONE  SEGMENT,  AVERAGE  DENSITY  OF  ELEMENT,  APPROXIMATE 
TANGENTIAL  STIFFNESS  INCLUDING  EFFECT  OF  CURVATURE,  and  frame  stiffness  using 
the  ASSUMPTION:  FRAME  = 5X  SHELL  STIFFNESS.  These  dependent  constants 

are  self-explanatory  in  lines  41  thru  72.  For  details,  refer  to  the  shell 
geometry  explained  in  the  context. 

SET  INITIAL  CONDITION  is  required  for  variables  UN,  UND,  UNDD,  UT, 

OX,  OY,  MX,  MY,  FN,  FT,  V,  QN,  QT  (lines  76  thru  94)  prior  to  initiation 
of  the  MAIN  DO-LOOP  OF  COMPUTATION  PER  TIME  STEP  NNNN  to  DEFINE  0-,  M-, 
and  V-  VALUES,  and  UN  AND  UT  IN  17  x 17  ARRAYS.  The  computational  steps 
in  the  DO-loops  are:  OX,  OY  INSIDE  BY  ALL  UN  AND  UT  INSIDE ; OX,  OY  AT 
BOUNDARY  BY  INSIDE  OX,  OY,  BOUNDARY  CONDITIONS  MX  AND  MY,  AND  APPROXIMATED 
OX,  OY  OUTSIDE  THE  BOUNDARIES;  MX,  MY  AND  FN,  V AT  BOUNDARY;  MX,  MY 
INSIDE  AND  AT  EDGE  BY  OX,  OY  INSIDE;  FN,  V INSIDE  AND  AT  EDGE  BY  MX  AND 
MY;  FN,  V OUTSIDE  THE  BOUNDARY.  An  addition  to  the  program  is:  FORWARD/ 

BACKWARD  CALCULATION  OF  FT,  an  alternate  approach  from  the  direct  one  to 
calculate  FT.  The  equations  are  physically  correct  assuming  a realistic 
distribution  of  FT.  Use  of  Equation  40,1  instead,  has  resulted  in  "blow- 
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up"  in  the  calculation  procedure,  CALCULATE  UNDD  AND  INTEGRATE  remains  an 
unchanged  procedure.  As  quoted  for  FT,  the  FORWARD/BACKWARD  CALCULATION  OF 
UT  replaces  Equation  40.2  for  reasons  of  "blow-up".  WRITING  OF  VARIABLES 
AND  END  OF  MAIN  DO-LOOP  completes  the  algorithm. 

4 . F-lll  WINDSHIELD  AND  PLOTTING  PROGRAM 
a.  Introduction 

This  computer  code  is  based  on  the  basic  shell  program  described  and 
listed  in  Appendices  III-3  and  IV-3.  The  finite  difference  grid,  dimen- 
sioned 17  by  17  in  the  present  case  (Figure  44) , correspondingly  requires 
as  many  17  by  17  arrays  of  computer  storage  as  variables  employed  at  grid 
locations  (Figure  45) , About  one  quarter  of  these  locations  is  actually 
utilized.  This  way,  simplification  of  index  notations  in  the  operational 
equations  is  achieved.  These  equations  could  be  rewritten  in  a compact 
form  to  eliminate  this  expenditure  of  storage.  A more  refined  analysis 
could  be  conducted  utilizing  the  mesh  completely;  however,  the  number 
of  calculative  steps  would  be  several  times  larger. 

At  a grid  size  of  17  x 17,  the  running  time  of  the  current  computer 
routine  was  less  than  30  seconds  per  set  of  input  data  for  the  selected 
parameter  printout  of  a 2500  step  history,  totaling  .25  seconds  of  real 
time*  Computer  time  is  reduced  when  printr-out  is  kept  at  a minimum  and 
not  required  at  every  time  interval.  In  addition,  non-sequential  print-out 
could  be  triggered  by  desired  threshold  values  of  any  timely  variable. 

b.  Time  Dependent  Variables 

The  variables  dimensioned  1=1  thru  17  in  radial  (X)  direction  and  J=1 
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Figure  45.  Typical  Locations  of  Variables  in  Grid 
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thru  17  in  longitudinal  (Y)  direction  (Figure  45)  are; 

QN(I,J)  normal  component  of  impulse  force  (pound) 

QT(I,J)  tangential  component  (pound) 

UN(I,J)  normal  displacement  (inch) 

UNMDT(I,J)  normal  displacement  of  previous  time  step  (inch) 

(not  used  in  current  option,  established  for 
control  of  attachment^conditions) 

UND(I,J)  normal  velocity  (in/sec) 

UNDD(I, J)  normal  acceleration  (in/sec^) 

UT(I,J)  tangential  displacement  (inch) 

0X(I,J)  radial  slope  of  deflection  (rad) 

0XMDT(I,J)  radial  slope  of  previous  time  step  (rad) 

0XD(I,J)  rate  of  slope  of  radial  deflection  (rad/sec) 

0Y ( I , J ) longitudinal  slope  of  deflection  (rad) 

0YMDT(I,J)  slope  of  previous  time  step  (rad) 

0YD(I,J)  rate  of  slope  of  deflection  (rad/sec) 

MX(I,J)  radial  bending  moment  (pound-in) 

MXEL(I,J)  elastic  component  of  bending  moment  (pound-in) 

MXVC(I,J)  viscous  component  (pound-in) 

MY(I,J)  longitudinal  bending  moment  (pound-in) 

MYEL(I,J)  elastic  component  of  beading  moment  (pound-in) 

MYVC(I,J)  viscous  component  (pound-ln) 

FN(I,J)  internal  normal  load  in  radial  direction  (pound) 

FT(I,J)  radial  tangential  load  (pound) 

V(I,J)  shear  load  in  longitudinal  direction  (pound) 
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Tlie  histories  of  the  normal  and  tangential  load  components,  QN 
and  QT,  at  time  step  NNNN  are: 

QT(I-1,J)  = SIN(THX(I-l,J)*Q(I,J)/2.  | 

I 

QN (I , J)  = C0S(THX(I,J))*Q(I,J)  I 

QT(I+1,J)  = SIN(THX(I+l,J))*Q(I,J)/2. 

The  present  computer  routine-  contains  the  Q(I,J)  data  of 
bird  impact  case  FM-1  as  a built-in,  but  flexible,  input  (see  l 

i 

Section  1V-4) . 

One  dimensional  arrays  of  order  17  for  definition  of 
normal  loads  outside  the  grid,  FWO,  FN18,  VO,  V18  are  defined  j 

T 

from  the  assumption  of  approximate  boundary  conditions  of  related  1 

variables  to  satisfy  the  number  of  equations  required  for  solution 
of  unknown  variables . 

c.  Fixed  Variables  I 

j 

Two-dimensional  arrays  of  order  17x17  are:  j 


DM(I, J) 

mass  distribution  (pound-sec^/in) 

DX(T,  J) 

finite  differences  of  radial  secant  (inch) 

DS (1, J) 

finite  radial  arch  difference  (inch) 

DY(I,J) 

longitudinal  finite  difference  distance  (inch) 

THX(I,J) 

radial  angular  spacing  (rad) 

DTH(I,J) 

related  angular  expression  (rad) 

SNDTH(I, J) 

angular  expression 

CSTHTH(I, J) 

angular  expression  (inch) 
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EIX(I, J) 
EIY (I,  J) 
KUN(I, J) 


radial  inplane  bending  stiffness  (pound-in^) 

o 

longitudinal  bending  stiffness  (pound-in  ) 


translational  spring  mounting  (pound/in) 

(only  utilized  at  the  edges) 

K0X(I, J)  rotational  sill  stiffnesses  (pound-in/rad) 

K0Y(I,J)  rotational  bow  or  arch  stiffnesses  (pound-in/rad) 

KUND(I,J)  viscous  translational  dampening  (pound-sec/in) 

K0XD(X, J)  viscous  rotational  sill  dampening  (pound-in-sec/rad) 

K0YD(I,J)  rotational  bow  or  arch  dampening  (pound-in-sec/rad) 

One-dimensional  arrays  of  order  17: 

RD(J)  radii,  variable  in  longitudinal  direction  (inch) 

CIRC(J)  arch  lengths  (inch) 

TH(J)  central  angles  (rad) 

SD(X)  variable  longitudinal  lengths  (inch) 


Determination  of  finite  differences  in  the  x-directlon  is  as  follows: 
RD(1)  and  RD(17)  are  minimum  and  maximum  radii,  RMIN  and  RMAX  respectively, 
assuming  a grid  size  of  17  (Figure  44).  Radii  of  bow,  RD(2),  and  arch, 
RD(16),  are  based  on  given  bow  and  arch  widths  in  proportion  to  the  size 
of  grid  in  the  y-direction.  The  radii,  RD(2)  thru  RD(16)  vary  linearly. 

All  the  other  radial  variables  are  as  indicated  in  Figure  46. 


T 
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I-0RD1NATE 
(ALONG  X-AXIS) 


DS=CIRCULAR  FINITE  DIFFERENCE  1 

DX=SECANT  FINITE  DIFFERENCE  J DS~DX  THX=ANGULAR  GRID  DISTANCE  FROM  £ (RD) 
DTH=ANGULAR  FINITE  DIFFERENCE  SNDTH=2 • SIN (DTH/ 2) 

CSTHTH-CURV*THICK*(1  - C0S(DTH/2))Z 


Figure  46.  Geometrical  Arrangement  of  Radial  Dimensions 
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The  types  of  mass  elements,  DM(1,J) , developed  at  the  corners 
(1,1),  (17,1),  (1,17)  and  (17,17)  are: 

DM(1,1)=DX(1,1) *MB0WtDY(1,1) *MS1L 

For  MBOW  and  some  c '.her  constants  in  the  following  equations,  see  the 
input  data 
Along  the  frame: 

DM(2,1)=DX(2,1)-MB0W 

Mass  properties  of  the  windshield  (window  pane)  calculated  from  windshield 
dimensions: 


DM(T,J)=SPDEN-DS(I,J) -DY(I,J) -THICK 

In-plane  pane  stiffnesses: 

EIX(I , J)=EP ANE • THICK3/12 • DY (I , J) 

EIY (I , J)=EPANE -THICK3/12 • DS (I, J) 

+CURV  • EPANE  -THICK-(1-C0S  (DTH(I , J)  /2)  ) ' RD  (J)  2 • DS  (I,  J) 

The  stiffnesses  along,  and  transverse  to,  the  frame  (sill,  head  beam,  arch, 
bow)  are  calculated  using  corresponding  eight  boundary  stiffnesses,  such 
as  EIARC,  EITAC . 

When  some  of  the  attachment  restraints  or  stiffnesses  are  not  effec- 
tive (equal  to  zero) , very  small  values  were  actually  used  to  stabilize 
the  calculation  procedure. 
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i 


d.  Input  Data 

An  example  input  for  operation  of  the  computer  routine  is  provided 
in  Table  XVII  showing  the  time  increment  of  each  step,  DT,  the  number  of 
calculative  steps  after  which  the  program  stops,  NEND,  and  the  frequency 
of  print-out,  NP.  CURV  is  a constant  to  account  for  additional  bending 
stiffness  due  to  curvature.  The  index,  ISHEAR,  set  equal  to  the  location 
of  centerline  of  impact,  arranges  the  string  lengths  of  "forward  and  back- 
ward"  solutions  of  tangential  equations.  C0R  considers  the  increased  stiffs 
nesses  in  the  corners. 

A description  in  sequence  of  the  remaining  variables  in  Table  XVII: 
THICK  thickness  of  window  pane  (inch) 

LHEA  length  of  head  beam  (inch) 

WHEA  width  of  heat  beam  (inch) 

LSIL  length  of  sill  (inch) 

WSIL  width  of  sill  (inch) 

RMIN  actual  radius  of  bow  (inch) 

THMN  minimum  central  angle  of  cone  segment  (rad) 

WB0W  width  of  bow  (inch) 

RMAX  actual  radius  of  arch  (inch) 

THMX  maximum  central  angle  of  cone  segment  (rad) 

WARC  width  of  arch  (inch) 

SPDEN  specific  density  of  window  pane  (pound-sec^/in^) 

MHEA  mass  of  head  beam  per  running  inch  (pound-sec^/in^) 

MSIL  mass  of  sill  (pound-sec^/in^) 

O n 

MB0W  mass  of  bow  (pound-sec  /in*) 


208 


AFFDL-TR-74-1 55 


MARC 

EPAME 

vise 

EIARC 

EITAC 

KUNAC 

K0ARC 

KUNDAC 

K0DARC 

EIHEA 

EITHD 

EISIL 

EITSL 

KUNSL 

K0SIL 

KUNDSL 

K0DSIL 

EIB0W 

EITBW 

KUNBW 

K0B0W 

KUNDBW 

K.0DB0W 


mass  of  arch  (pound-sec2/in2) 

Young's  modulus  (pound/in2) 

pane  viscous  constant  (pound-in2 /rad) 

bending  stiffness  of  arch  (pound-in2) 

transverse  stiffness  of  arch  (pound-in2/ln) 

spring  mounting  of  arch  (pound/in) 

rotational  arch  stiffness  (pound-in/rad) 

viscous  translational  arch  dampening  (pound-sec/in) 

viscous  rotational  arch  dampening  (pound-in-sec/rad) 

bending  stiffness  of  head  beam  (pound-in2) 

transverse  stiffness  of  head  beam  (pound-in) 

bending  stiffness  of  sill  (pound-in2) 

transverse  stiffness  of  sill  (pound-in) 

spring  mounting  of  sill  (pound/in) 

rotational  sill  stiffness  (pound^in/rad) 

viscous  translational  sill  dampening  (pound-sec/in) 

viscous  rotational  sill  dampening  (pound-ln-sec/rad) 

bending  stiffness  of  bow  (pound-in2) 

transverse  stiffness  of  bow  (pound-in) 

spring  mounting  of  bow  (pound/in) 

rotational  bow  stiffness  (pound-in/rad) 

viscous  translational  bow  dampening  (pound-sec/in) 

viscous  rotational  bow  dampening  (pound-in-sec/rad) 
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e.  Calculation  Procedure 

The  calculations  start  by  zeroing  the  variables  for  initial  condi- 
tions (1=1,17  and  J=l,17):  UN,  UND,  UNDD,  UT,  OX,  OY,  MX,  MY,  FN,  FT,  V. 

Load  inputs  to  initially  stabilize  equations  are  quasi-zero  loads  of  ' 

.01  pounds  everywhere  on  the  grid  for  QN,  and  QT  on  the  right  side 
(1=1,8),  QT=0  at  the  center  line  (1=9),  and  QT=  -.01  on  the  left  side 

(1=10,17).  The  calculation  continues  similar  to  the  "flat  plate"  j 

I 

routine  of  Appendix  IV-1.  Plotted  results  for  this  example  may  be  found 
at  the  end  of  the  listing  of  this  program  (Appendix  IV-4) : the  three 

page  17  x 17  matrix  print-out  per  time  step,  NNNN,  at  .025  sec  and  .05  1 

sec  (DT  = .0001,  NEND  = 500,  NP  = 250).  The  print-out  of  the  selected 
strain  and  deflection  histories  (of  case  "1",  see  page  272)  is  self- 

i 

explanatory.  Here,  NP  = 5 yields  100  values  per  column.  For  further  j 

details,  the  user  is  referred  to  the  comments  in  the  program  listing 

which  are  practically  identical  with  the  basic  conical  shell  program  that  j 

was  adapted  to  this  bird  strike  case,  EM-1.  One  then  simply  follows  the  j 

F0RTRAN  statements  in  both  the  listings  (Appendix  IV-3  and  4).  j 


i 
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APPENDIX  IV 

COMPUTER  PROGRAM  LISTINGS 

The  attached  F0RTRAN  IV  source  programs: 

1.  FLAT  PLATE,  one  starting  at  page  222,  the  other  at  231, 

2.  PLATE  WT-18  (page  236), 

3.  CONICAL  SHELL  SEGMENT  (page  245),  and 

4.  F-lll  W/S  (FM-1)  on  page  256, 

contain  built-in  input  parameters  and  flexible  input  data.  A modest  amount 
of  printed  output  has  been  generated  at  selected  time  steps. 

The  source  decks  of  the  computer  programs  described  herein  can  be 
obtained  by  contacting  AFFDL  (FBS/JANSENX,  W-PAFB,  OH,  45433,  (513)255-5548. 
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